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SSPARAMA:  A NONLINEAR,  WAVE  OPTICS  MULTIPULSE 
(AND  CW)  STEADY-STATE  PROPAGATION  CODE  WITH 
ADAPTIVE  COORDINATES 


INTRODUCTION 

Several  methods  of  propagating  CW  high-energy  laser  beams  through  the  atmosphoiv- 
have  been  reported  previously  [1,2].  This  report  will  describe  a method  for  propagating 
multiply  pulsed  laser  boams  in  a nonlinear  atmosphere  by  adapting  the  coordinate  system 
to  the  amount  of  thermal  blooming.  This  technique  increases  the  accuracy  of  thermal 
blooming  calculations  and  extends  the  capability  of  the  code  in  the  case  of  extreme  hcuu- 
distortion. 

The  computer  code  SSPARAMA  calculates  the  steady-state  intensity  pattern  of  n 
train  of  high-energy  laser  pulses  propagating  through  the  atmosphere  in  the  presence  >>r 
thermal  blooming.  Steady  state  is  achieved  when  enough  equally  spaced,  equal-energy 
pulses  have  been  propagated  for  transients  in  air  heating  to  have  died  out.  In  the  stoutly 
state  a single  pulse  will  propagate  in  an  atmosphere  that  has  been  heated  by  many 
preceding  pulses  which  have  the  same  energy  distribution  as  the  pulse  one  is  calculating. 
The  pulse  widths  are  assumed  to  be  short  compared  to  the  sound  transit  time  across  the. 
face  of  the  beam,  so  that  self-blooming  will  not  take  place.  Blooming  occurs  only  as  a 
result  of  air  heating  by  preceding  pulses.  However,  to  avoid  problems  of  plasma  forma- 
tion, the  pulse  width  must  be  sufficiently  long  that  the  critical  intensity  for  air  break- 
down is  not  exceeded.  Finally,  us  the  pulse  is  propagated  from  one  coordinate  plane  to 
another,  coordinate  transformations  are  performed  to  insure  that  the  transverse  scale 
lengths  are  adapted  to  the  amount  of  thermal  blooming  induced  on  the  pulse  train  by  the 
negative  lensing  Influence  of  the  heated  atmosphere. 

Another  requirement  for  steady-state  propagation  is  that  a cooling  mechanism  exist 
for  removing  heated  air  from  the  path  of  the  beam.  In  SSPARAMA,  cooling  is  provided 
either  by  a wind  moving  perpendicular  to  the  propagation  direction  or  by  beam  sluing 
about  an  axis  in  the  aperture  plane  perpendicular  to  both  the  wind  and  the  propuguiirm 
directions.  Tire  steady-state  density  changes  Ap  Introduced  in  the  path  of  a given  pulse 
by  energy  absorption  from  all  preceding  pulses  can  then  be  expressed  as  [3] 

_ i 0.  2 

Ap  = — aEpe-a*  ]T  \<p(x-  nbl,(vo  + Or),  y,*)|  , (1) 

c»“  M- 1 
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whore 

x ■ the  distance  in  tho  propagation  direction  measured  from  the  aperture  plane, 

x " the  distance  In  the  wind  direction  measured  from  beam  maximum  intensity 
in  the  aperture  plane, 

7 ■ the  ratio  of  atmospheric  specific  heats  (*  1.4), 
cg  ■ the  speed  of  sound  in  air  {*340  m/s), 
a ■ the  absorption  coefficient  for  the  laser  radiation, 

Atj  * tho  pulse  spacing, 

Ep  ■ the  energy  of  each  laser  pulse, 

vq  * the  wind  speed  along  the  x direction  perpendicular  to  the  direction  of 
propagation,  and 

12  - the  angular  sluing  rate  of  the  beam  about  the  y axis. 

Finally  0 is  the  normalized  steady-state  energy  distribution  of  each  pulse  at  the  x plane: 

J \<t>(x,  y,  x)\2  dxdy  ■ 1.  (2) 


This  density  reduction  Ap  changes  the  index  of  refraction  from  its  ambient  value  no, 
where  no  * 1,  to 


na  * n0a  + 3 NAp, 

where  N is  the  molecular  refractivity  of  air  (*0,164  cm3/g).  The  distribution  0 must 
then  be  calculated  self-conslstently  from  the  propagation  equation: 


o, 


(3) 


where  It  - 2ir/\  is  the  wavenumber  of  the  laser  radiation.  It  is  assumed  in  SSPARAMA 
that  at  x ■ 0 the  pulse  train  has  a spherical  phase  front  and  a truncated  intensity  profile. 
For  example,  when  truncated  Gaussian  pulses  are  propagated 

0(.t,.V,O)  - Nu<pK(x,y),  xa  + .v2  < 2oa, 

■ 0,  xa  + ,v2  > 2a2,  (4) 


where 


0( t(x,y) 


_L~  l l(.va*lv1,)/n»  1/2 
Oy''? 


(6) 
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and  Nu  Is  a normalization  constant  insuring  that  Eq.  (2)  is  satisfied  at  z * 0,  Two  scale 
lengths,  a and  f,  are  defined  in  Eq,  (5).  The  scale  length  f , the  initial  curvature  of  the 
phase  front,  defines  the  distance  from  the  aperture  to  the  focal  plane.  At  a distance  a 
from  the  aperture  center  the  beam  intensity  falls  to  1/e  of  its  maximum  value,  and  the 
beam  is  truncated  at  1/e3  of  maximum  intensity. 

Altogether  eight  variable  physical  quantities,  a,  f , h,  a,  Ep,  At,,  vq,  and  *1  appear  in 
Eqs.  (1)  through  (5).  All  variations  will  not  however  lead  to  a mathemuttcidly  distinct 
problem.  In  SSPARAMA  Eqs.  (1)  through  (5)  are  scaled  so  that  distinct  propagation 
problems  are  defined  in  terms  of  five  dimensionless  parameters.  The  progre  . io  designed 
to  accept  either  the  set  of  data  with  dimensions  or  the  dimensionless  set,  end  both  sets 
are  printed  out. 

The  scaling  of  Eqs.  (1)  through  (5)  is  carried  out  via  the  coordinate  transformations 

* ■ J.  Sf  85  “•  * s * (6) 

and  the  variable  transformation 

$>(2.?.  *)  H a0(#,y,z).  (7) 

By  multiplying  Eq,  (3)  through  by  a3,  one  can  write  the  propagation  equation  in  a form 
which  identifies  the  five  dimensionless  parameters  characterizing  propagation  in 
SSPARAMA: 


0-0. 


(8) 


and 


, N„ , and  N„  are  defined  as 

N*  - he2  If, 

(9) 

„ 3/Vfc(7  - DotfEp 

Nc  c,2a2 

(10) 

Na  “ af, 

(ID 

_ 2a 
0 ('oAt,’ 

(12) 

n,  - nf/VQ, 

(13) 

Nk  1s  the  Fresnel  number  of  the  f roe-propagation  problem,  and  Ne,  Na,  N„,  and  N,  are 
coupling  strength,  absorption,  overlap,  and  sluing  parameters  respectively.  was  intro- 
duced by  Wallace  and  Lilly  [4]  and  called  the  pulses-per-flow-time  parameter.  It 
measures  tho  number  of  preceding  pulses  which  have  heated  the  air  across  the  beam 
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aperture  as  the  pulse  under  study  begins  to  propagate,  The  solution  to  Eq.  (8)  is  ob- 
tained subject  to  the  energy  normalization 


Il$(xty,0)ls  dxdy  “ 1 (14) 

and  the  initial  condition 

0(S,y,O)  - |0|  (IB) 


where  |$|  ■ 0 for  52  + ?2  > 2. 

Equations  (8),  (14),  and  (15)  are  numerically  solved  in  SSPARAMA  on  a 64-by-64 
grid  in  the  xy  plane.  Since  one  would  like  to  use  as  much  of  the  computational  grid  as 
possible  to  describe  the  variations  in  beam  intensity,  a scheme  for  adapting  the  coordi- 
nate grid  to  the  propagation  must  be  used.  For  example,  as  the  beam  propagates,  the 
initial  focusing  causes  the  beam  intensity  pattern  to  decrease  in  size  until  the  negative 
lenslng  effects  of  the  heated  utmosphere  accumulate  to  thermally  defocus  it.  Moreover, 
since  the  wind  removes  heated  air  from  the  path  of  the  beam  from  left  to  right,  a 
thermal  gradient  is  established  that  deflects  the  beam  from  right  to  left.  If  the  computa- 
tional grid  were  not  moved  or  changed  In  size  as  the  beam  intensity  was  calculated  from 
aperture  to  focal  plane,  the  intensity  pattern  would  either  be  poorly  sampled  as  it  de- 
creased in  size  or  it  would  expand  or  deflect  to  reach  the  boundary  of  the  grid  and 
Invalidate  the  calculation. 

A technique  for  adapting  the  computational  grid  to  local  changes  in  the  size  or 
location  of  the  beam  intensity  pattern  has  been  developed  by  Herrmann  and  Bradley  [5], 
A slightly  modified  form  of  their  technique  has  been  incorporated  into  SSPARAMA  and 
will  be  described  in  the  next  section  of  this  report.  In  the  third  section  the  numerical 
procedures  used  in  SSPARAMA  will  be  described,  and  in  the  fourth  section  the  code 
usage  will  be  explained. 


COORDINATE-SYSTEM  ADAPTION 

The  dimensionless  form  of  the  propagation  equation  cun  be  rewritten  more  com- 
pactly as 


+ a*  + a?  +/,'V(«2  i)]5  - o,  US) 

where  n2  - 1,  the  nonlinear  index  of  refraction,  depends  on  $ as  given  by  Eq.  (8).  The 
$$£  coordinate  system  is  normalized  to  the  constant  lengths  n and  f,  and  is  fixed  in  space. 
In  this  system  therefore  the  beam  will  lie  symmetrically  ubout  the  origin  of  the  .v.v  plane 
only  at  i m 0 with  an  extent  of  order  1 (see,  for  example,  Kq.  (15)).  When  z 4=  0,  a new 
set  of  xy  coordinates  is  needed  to  maintain  the  two  properties  that  the  beam  be  centered 
about  the  xy  coordinate  origin  and  be  of  order  1 in  extent.  In  general,  one  can  relate 
the  xy  and  xy  coordinates  by  a set  of  scale  parameters  D\  and  Da  and  a deflection 
parameter  X , which  ore  functions  of  z,  Since  one  would  like  to  solve  Eq.  (16)  in  a set  of 
coordinates  that  adapt  to  changes  in  beam  size  and  direction,  the  coordinate  transformation 


NHL.  REPORT  8074 


must  be  related  to  these  beam  changes  as  determined  by  the  linear  and  quadratic  terms 
of  the  phase  front.  By  analogy  therefore  with  the  transformation  to  dimensionless 
parameters,  one  must  perform  simultaneous  coordinate  and  variable  transformations. 
The  form  of  these  transformations  is  suggested  by  linear  propagation  theory: 


x-X 

(17) 

y 

(18) 

~N~k  ’ 

(19) 

and 


S . — .va t +•?!» ) (2 

The  constant  scale  change  from  * to  z la  done  for  convenience  to  eliminate  Nu  from  the 
4-derivative  term  In  Eq.  (18): 


2 iNkbi  -*  219*. 

The  factor  \j^/D\D%  is  removed  from  $ to  insure  the  form  invariance  of  the  energy 
normalisation: 


1. 


(21) 


When  Eqs.  (17)  through  (20)  are  substituted  into  Eq.  (16)  and  when  the  nonlinear 
term  is  of  negligible  size  and  the  beam  has  a Gaussian  profile,  D i,  Dg,  X,  di , dg,  j3,  71 , 
and  7a  os  functions  of * can  be  analytically  determined  for  all  z.  However,  when  the 
nonlinear  term  is  important  or  when  a non-Gaussian  beam  is  propagated,  the  d’s  and  (J , 
which  represent  the  effective  quadratic  and  linear  phase  changes  throughout  the  xy  plane, 
can  no  longer  be  so  determined.  One  must  adopt  a more  limited  strategy  for  the  em- 
ployment of  Eqs.  (17)  through  (20). 


Consider,  for  example,  that  the  quantities  0i,  Dg,  X,  di,  d2,  jj,  7i,  and  72  are 
known  at  z ■ r0  and  that  their  dependence  on  t is  to  be  analytically  determined  as  one 
propagates  to  a neighboring  xy  plane  at  + Az.  Since 


_L 

D\ 


9*. 


(22) 


± 

o2 


9 


2 

.v  ’ 


(211) 
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mid 


one  finds  that 


9; ' * [s' - (f s' ln 01  * ~§i)' 9*  ■ s' 9* ln °2 1 *4 


[2 tNnl)£  + 3 2 + 3|  + *V(»2  - 1)]---^— 

V^1^2 


| .{* +fl*.v  * +/J.v+f  i ♦fa ) 


(24) 


. 2i  fbt  - | 3,  In  D,  3,  - Wv  - | 3,  In  D2  3vj 


i(3,lftO,  + 3,ln»2)  - 23,(71+^)  + ~ 3.v  ~ [23,^#  + *)  + 0]2 

4 /Jl 


+ (2 «t  (v/?J,  * ♦ tf)  + (J]3V  + 2(S,  + -f  3*  - 4fi22D2y2 

V^I  iJ2  ' 


+ 4 i&n.vdy  + 2f«2  + fe2(j2(/i2 


^ - 0. 


(25) 


For  vunishingly  srnull  /ia  - 1 und  for  tt  reul  Gaussian  profile  ^(.v,  y,  z0)  one  would  deter- 
mine  Du  Da.  X,  Si,  da.  ft,  1u  and  72  from  the  requirement  that  Eq.  (25)  be  cupuble  of 
being  put  in  the  form 


2(3,  + -f  (32+l 
*'1 


.v2)  + 


— ( i)y  + 1 ~ y2)  + Ataua(#i2  - 1) 
lh2  ' 


\p  “ 0. 


(20) 


Then,  os  wus  propagated  to  to  + Az,  it  would  acquire  no  z dependence  und  would 
remain  real  and  Gaussian;  thut  is,  all  of  the  z dependence  of  0 would  have  been  accounted 
for  in  D\,  • ■ 7a- 


For  the  imaginary  terms  of  Eq.  (25)  other  than  21 3,  to  vunish,  the  quantities  D\, 
D a.  and  A\  which  determine  the  scale  and  location  of  the  xyz  coordinate  system,  must 
satisfy  the  equations 


3,  In  D 1 

" 4d) , 

(27) 

3^  In  D2 

- 4 da, 

(28) 

'0tX  n 2di,V  + ft. 

(29) 

On  the  other  hand,  for  the  real  terms  involving  0.v  und  3V  to  vanish  und  for  the  scale 
functions  /)|  and  Da  to  be  factorable  from  the  remaining  ,v  und  y terms  respectively,  the 
phase  functions  at,  da.  (5,  ylt  und  72  must  sutisfy  the  set  of  equations 


fl 
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2/}1d,a1  + 45 ,2/),  - , 

I}\ 

(30) 

2/)23j52  + 432,"fJ2  “ 1 

D 2 

(31) 

3,0  + 2*3,5,  + 25,(25,*  + $)  - 0, 

(32) 

23,7!  + 2*z3,5,  + 2*3,$  + (25, *+$)z  - - — , 

D\ 

(33) 

^■'37 

(34) 

Thus  Eqs.  (27)  through  (34)  will  determine  oil  of  the  z dependence  of  $ when  \p(x,  y,  z0) 
Is  real  and  a Gaussian  function  of  x and  y and  there  Is  no  lensing  effect  cuused  by  heating 
of  the  atmosphere;  that  is,  Eqs.  (27)  through  (34)  will  describe  beam  focusing  in  the 
absence  of  diffraction  and  nonlinear  media  phenomena.  They  are  of  more  limited  utility 
when  such  phenomena  are  present.  In  this  case,  during  the  displacement  of  <t>  from  z0  to 
z0  + Az,  linear  and  quadratic  phase  changes  will  arise  from  two  sources.  As  a result  of 
focusing  ut  z • z0,  the  Initial  phases  5i(zo),  62{zo),  and  (J(*o)  will  become  5](zo  + Az), 
Sa(zo  + Az),  and  /3(z o + Az)  through  the  solution  to  Eqs.  (27)  through  (34).  In  addition 
however  \p  at  zo  + Az  will  acquire  linear  and  quadratic  phases,  A/?,  A5,,  and  Aa2  respec- 
tively, as  a result  of  diffraction  and  thermal  blooming.  Thus  at  zo  + Az  a new  factoriza- 
tion of  $ must  be  made,  namely, 


3>(X  V i0+Asl=  ^ (*t.V,iO_+Az)_  e(|5i(<0+AzXva*dii(Zo+Az).vi<+^(*o+A*).v+7i+7i| 

<i/D[U0*Az)D2(z0*Az)  ^ 

if  i//,  which  is  to  be  propagated  from  zo  + Az  to  zo  + Az  + Az', 
quadratic  or  linear  phase  terms,  After  each  step  in  propagation 
must  bo  redefined  as 

is  not  to  initially  have 
therefore  5,,  a2,  and  $ 

5',(z0  + Az)  - 5j(z0  + Az)  + AS,, 

(30) 

«!j(zo  + Az)  - S2(z0  + Az)  + Aa2, 

(37) 

and 

$'(*o  + A z)  - j}(zQ  + Az)  + A$ 

(38) 

In  order  to  adapt  the  coordinate-system  determination  from  Eqs.  (27)  through  (29)  to 
changes  in  phase  that  result  from  focusing,  diffraction,  and  thermal  blooming. 
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In  SSPARAMA,  ^ is  propagated  from  one  z plane  to  another  by  finite-differencing 
a phase-transformed  version  of  Eq.  (28).  Then  Aai,  Aa2»  and  A/)  are  found  in  the  xyi 
coordinate  system  using  the  method  of  phase  minimization  discussed  by  Herrmann  and 
Bradley  [5] . One  requires  that 

f |i//|2[V(A<r1*2  + Aa2y2  + Aj3x  - 7)]2dxdy  - minimum,  (39) 
*«*o+A  * 

where  \j/(x,  y,  zo  + Az)  s It  follows  that 


and 


where 


and 


Aon 


A(3 


Act 


D\E  ~ B1C1 
2(Atf?-Bi2)  ’ 

(40) 

A,Cj  - BXDX 
AXE  - Bx 

(41) 

'*  " 2Aa  ’ 

(42) 

Ai  a J'#2 10|a  dxdy,  A2  » Jy2li^la  dxdy,  (43) 

Bx  S jxm2  dxdy,  (44) 

Ci  s Im  Jij/*dx\p  dxdy,  (45) 

7?i  h Im  |x^/*  9*1// dxdy,  D 2 s Imjy\p*dy[p  dxdy,  (40) 

F,  & jty\2  dxdy  - l.  (47) 


The  factorization 


i//(x,  y,  z0  + Az)  ==  0'ei(Aat*a+dtt3/',+A|3*)  (48) 

will  then  define  \p'  at  zo  + Az  as  a wave  function  of  minimum  quadratic  and  linear  phase. 
In  particular,  if  \p  is  exactly  0 Gaussian  beam,  \I/‘  will  be  real. 
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The  relationship  between  (Acq,  Ac^, 
tuting  Eqs.  (17)  and  (18)  into  Eq.  (48): 

AS] 

Aa2 


A(3}  and  {ASi,  A&2,  A/J)  is  found  by  substi- 


Ac*] 
~z  » 

(49) 

0 1 

Attn 

(50) 

d2 

and 


A0 


Vd~i 


2A  otiX 
Ox 


(51) 


A similar  set  of  equations  will  hold  between  {5i,  Sa,  (J } and  (ai,  c*2.  M < which  are 
computed  directly  in  the  xyz  coordinate  system.  When  reexpressed  in  terms  of  ai,  a2 
and  |3.  Eqs.  (27)  through  (29)  become 


BtDt  ■ 4oq, 
btD2  “ 4oi2> 


and 


dsX 


L_ 

Vox' 


and  Eqs.  (30)  through  (32)  transform  into 


■ — (1+40:!2), 


" ^,-(1+4CV22)' 


and 


5,(3  “ 


2o<l0 

D, 


(62) 

(63) 

(54) 

(55) 

(56) 

(67) 


Eqs.  (52)  through  (57)  must  be  solved  in  terms  of  Initial  values  at  *q.  The  solutions  arc 


0 1,2(*)  “ ^i,2(*o) 


^i,a(*o) 


2 2 ■> 
+ J-f!L  1 
.•®i.a(*o)J  J 


(58) 
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«t,2(*)  “ “i,2(*q)  + -J-U  + l2ali2U0)]2} 


P(z)  » /3(20)- 


1 + 


2ft  i (*q  ) 
®1  (*q) 


ll2  + 

~z  - 2q" 

1 

bZ>i(«o). 

and 


(3(*o) 

X(z)  - X(z0)  + -7====  (*  - *0)- 

Finally  the  procedure  for  solving  Bq,  (26)  In  SSPARAMA  is  similar  to  the  one 
described  in  an  earlier  report  [2].  A phase  transformation  on  ^ is  made: 


(59) 

(60) 

(61) 


where 


«!»(*,  y,  2)  sb  4/(x,  y,  z)  1 , 


g(x,y,z)  m (l-*a)  + ~(1  - y2)  + *2aV  - 1). 


The  equation  for  <1>  follows  from  Eq.  (26): 

[2/9,  + H(x,  y,  *)]<!>  * 0, 

where 


(62) 


(63) 


(64) 


(66) 


By  picking  z b to  lie  between  zq  and  2o  + Ar,  one  can  propagate  lt>  from  zq  to  z o + A20, 
with  first-order  accuracy,  by  solving  the  equation 

[2i9,  + H(x,  y,  2q )] - (2it,  + ~ 9 *2  + “9 , ,2)*  - 0.  (86) 

Equation  (68)  Is  solved  by  Fourier  transforming  <1>  [6] , 

<$(/«!,  ft2,20)  « J*«*i*+*»y)*(*,y1a0)djCdy,  (67) 

and  propagating  $ to  zQ  + Az: 


I 


M 


MMkUMlIUk 
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The  inverse  transformation  to  Eq..(67)  then  yields  <1*,  and  Kq,  (62)  yields  \p(x,  y, 
zq  + A z), 


NUMERICAL  PROCEDURES 


The  phase  function  g(x,  y,  z)  of  Eq.  (63)  can  be  written  more  usefully  in  the  form 


Jl(*)  + St(y) >»  *) 

Di(«)  C2(*)  V^](*)D2(z)  ' 

where 

*i(*)  si-*2, 

£a(y)  s 1 - y2, 

and 


(89) 


(70) 

(71) 


*a(*.  y.  *)  5 NkNce~NaN>l* 


?>[*- 


2n 


N0y/D&) 


(l+N,Nkz),y,z 


(72) 


This  expression  for  #3  is  found  by  substituting  the  new  variables  x,  y,  z,  and  4>  into 
Eq,  (8).  The  phase  integral 

A 0 a I g{x,y,z')dz' 

**o 

appearing  in  Eq.  (62)  can  now  be  partially  evaluated  and  expressed  in  the  form 


AO  - gx(x)AZx  + g2(y)AZ2  - f (73) 

J,’  ^Dx(z)D2(z) 

where 


A^ia 


dz' 


tan/i'1  (d  + £2of!  2(-*o)l 2 ) + 2“1 


,2(*0  )j 


* ■* 


• (74) 


* **o 


The  differential  quantities  A Zx  and  A Z2  are  similarly  named  as  the  coordinate  differ- 
ential A Z that  was  used  in  earlier  code  calculations  which  involved  only  a single  scaling 
function  D(z). 

To  complete  the  evaluation  of  AO,  one  must  know  the  z dependence  of  £3,  that  is, 
the  z dependence  of  |<I>|2.  Two  options  are  provided  in  SSPARAMA,  for  evaluating  AO, 
depending  on  whether  one  has  determined  l«I>la  at  one  or  both  of  the  integration 
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endpoints.  The  procedures  work  as  follows:  Suppose  first  that  the  solution  for  ip(x,y,zo) 
has  been  obtained.  Then  one  can  compute  g(x,  y,  zo),  since  lfl>(;e,  y,  zo)l2  ■ IM#»  y,  *o)l2* 
To  find  <!>(*,  y,  »o  ),  however,  one  must  evaluate 

r *q 

AO'  h g(x,  y,  z‘)  dz',  (75) 

where  zo  lies  between  zq  and  the  plane  zo  + Ax  to  which  one  would  like  to  propagate  0. 

If  ^ is  known  only  at  zp,  the  zoroth-order  approximation 

A O'  * gxMAZ'i  + g2(y)AZ'2  ~ gi(x,y,z0)AZ,li  (70) 


must  be  made,  where 


AZ'X2  3 


r*°  - dz' 

\ Vd&WW] j ' 


(77) 


Equation  (66)  can  now  be  solved  for  <!>(*,  y,  zq  + At)  by  the  use  of  Fourier  trans- 
formations. Finally  on  performance  of  the  phase  integral 


r* o+A* 

AO"  a I ( g(x,  y,  z')  dz'  (78) 

^(x,  y,  zo  + Az)  can  be  obtained  from  «!>(*,  y,  zo  + Az).  In  keeping  with  the  accuracy 
with  which  A O'  was  approximated,  AO"  can  be  approximately  evaluated  as 

AO"  m gx(x)AZ'i  + g3(y)AZ'i  - g2(x,  y,  z0  + A z)AZ'[2.  (79) 

The  differentials  AZ",  AZ2,  and  AZja  are  defined  by  the  integrals  of  Eqs.  (74)  and  (77) 
with  the  integration  limits  as  specified  in  Eq.  (78). 

Suppose  however  that  initially  both  \p{x,  y,  zo)  and  4/(x,  y,  zo)  ore  known  and  that 
the  values  of  \p  at  zo  are  to  be  propagated  to  the  plane  at  zo  + Az.  In  this  cuse  tho 
phase  integrals  defined  in  Eqs.  (75)  and  (78)  can  be  approximated  using  the  integration 
formula 


c 


■*0  + A* 


f(x)g(x)dx  * u>if(x0)  + w2f{x o+A*), 


(80) 


which  has  first-order  instead  of  zeroth-order  accuracy.  The  weights  u>i  and  w2  are  thus 
determined  such  that  equality  will  hold  In  Eq.  (80)  whenever  f is  a linear  function  of  x\ 


vt 


1 + 


5)f 


x0+Ax 


g(x)  dx 


2_ 

Ax 


*n  + 


xg(x)  dx 


(81) 
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and 


W 2 


- JL 

Ax 


r*( 


*0  + Ax 


xg(x)dx 


2x0  /'*o+^ 


_5  | 

Ax  J 


#(x)  dx. 


(82) 


Than,  for  example,  in  place  of  Eq,  (78)  one  would  have  that 

A O'  * gi(x)AZ\  + f2(y)A#a  - #3(*.  y,  *0)^3  ” #3(*>  y>  *o)AZ\,  (83) 

where  AZ'3  and  AZ\  are  related  through  Eqi.  (81)  and  (82)  to  AZ\^  and  an  integration 
over  the  function  *VZ)i(«yJDa(*) : 


and 


A z'3 


*0  * * 0 A„,  _ 2 f t'  dt' 

*0  ” *0  18  *0  “ *0  *t0  \/f3i(*')Oa(2’) 


A Z\ 


2 f f *°  a'd*' 

*b  - *0  IA0  v'i>i(i')X>aUr) 


- XqAZ\2 


(84) 


(86) 


Although  integration!  over  Df1  and  Dj1  can  be  carriod  out  analytically  in  terms  of 
inverie  hyperbolic  tangent!  (ai  in  Eq.  (74)),  integral!  over  lA/ETT^a  produce  elliptic 
function!.  Both  set!  of  Integrations  are  handled  in  SSPARAMA  numerically,  with  third- 
order  accuracy,  using  a second  Integration  formula: 


, *o+Ax 


J f(x)dx  * j[A*o  +Ax,)  + f(x0  + Axa)), 


(88) 


where  Axi  s (1  - lfs/TS)Ax/2  and  Axa  B (1  + l\/5)Ax/2.  Again,  as  an  example,  con- 
sider Eqs.  (84)  and  (85)  and  define 


f\  55 


(87) 


and 


/ sx * 

v'^,i(*a)°a(2a) 


(88) 


where  *1  s «o  + (1  _ l/\/3)t(*b  _ ■*o)/2]  and  »%  = *0  + (1  + l/v/3)t(«o  - «o)/2].  One  can 
complete  the  numerical  evaluation  of  AZ']  and  A Z\  bv  rewriting  Eqs,  (84)  and  (85)  with 
the  use  of  Eq.  (86),  in  terms  of  f\  and 
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A*S»  “ ~h]  (89) 

and 

- (*0-«0)(^7^-)  - A*3-  (90) 

The  procedure  by  which  Eqs.  (80)  through  (90)  are  employed  requires  that  two  sets 
of  values  of  \ p be  stored  at  any  time  by  SSPARAMA.  At  the  beginning  of  the  propaga- 
tion step  described  above,  the  two  arrays  contain  the  values  of  i/»(x,  y,  zq)  and  i//(x,  y, 
z'o),  where  *o  < *o  < *o  + Ax.  At  the  end  of  the  propagation  step  the  values  of 
0(x,  y,  z o)  have  been  roplaced  by  <//(«,  y,  z o + Ax),  These  new  values  con  then  be  used 
to  propagate  \p(x,  y,  xq)  to  ip(x,  y,  xb  + Ax'),  where  now  *o  < *o  + A*  < xb  + Ax’.  The 
process  of  alternatively  propagating  one  and  then  the  other  of  the  two  arrays  is  repeated 
until  the  focal  plune,  defined  by  the  Initial  beum  curvature,  is  reached. 

Since  both  arrays  are  initially  assigned  the  values  i//(x,  y,  0),  the  process  of 
propagating  one  array  past  the  other  cannot  begin  until  after  the  first  propagation  step. 
The  first  z step  is  therefore  taken  using  Eqs.  (76)  and  (79)  to  determine  AO'  and  AO”.  In 
general  the  incremental  steps  Ax  are  selected  in  SSPARAMA  according  to  a criterion  that 
the  phase  changes  induced  by  g3  os  computed  from  Eq.  (70)  be  no  larger  than  some  pre- 
assigned value  of  order  1 for  all  x and  y.  However,  to  carry  out  the  first  advancement  of 
^ nt  xo  * 0,  half  of  the  initially  computed  Ax  value  is  used.  This  leapfrog  procedure  is 
summarized  for  the  first  few  z steps  in  Fig.  1. 

The  advantage  conveyed  by  using  Eqs.  (70)  and  (79)  to  evaluate  the  phase  integrals 
A O'  and  A 0"  is  that  only  one  <p  urray  is  needed  in  carrying  out  the  calculation.  Because 
of  the  reduced  accuracy  In  computing  A O'  and  AO”,  however,  smaller  x steps  are  in 
principle  required  to  obtain  the  same  results  as  when  two  urrays  at  different  x planes  are 
used.  To  allow  a quantitative  comparison  of  these  two  procedures,  both  options  for 
propagating  ^ were  installed  in  SSPARAMA  and  can  be  selected  according  to  the  value 
of  one  of  the  input  parameters  to  the  code.  For  the  same  reason,  another  input  param- 
eter is  also  available  that  allows  one  to  adapt  or  not  adapt  the  coordinate  system  to  the 
amount  of  diffraction  or  thermal  blooming  occurring  during  beam  propagation. 


PROGRAM  OPERATION 

This  section  will  describe  the  input  parameters  required  to  run  SSPARAMA  and 
explain  the  data  included  In  the  output.  A complete  listing  of  SSPARAMA  is  included  in 
Appendix  A. 

To  use  program  SSPARAMA,  two  input  cards  are  required.  The  first  specifies  cer- 
tain numerical  parameters  und  selects  vnrious  program  options,  and  the  second  defines  the 
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*T1P»0  < 


DWlNI‘f|«<y 
NMI  CHANQI 
IIT  «i  - »|  • WNfc, 

AMD  »ITfl*  A*,  >i4,  • A1-G. 


■IIP  1 < 


•TIP  i d 


•TIP  i d 


1-0 


t-ki 0 


AdVBfMt 


NOW  PIND  a,,  •},  AND  | 
ANO  U\,  Ut  AMD  ip, 

Ul«  TO  apply  PHAII 
C HANOI  l^ll-Ol-Ajd-Ok 


1*0 


»-*A> 


1-0 


Afraof  A|jli  ii|  ♦ lui 

#‘14*1  PINO  AND /I  ANO 
/U,,  4«)  AND  if  AOAlN 
AND  till  TO  APPLY 
PHAII C HANOI 

Y,IAl*J-PlUlcl. 

<-A4q  A - Alp  « Alt  f-A*»Al,  « Al| 

*'IU»I  *»IM#  • A*ll  UNO.,,.,,  ANOfl 

AND  itt|,  A-.,  AND  M 
AMlY  l-HAtl  CMANQI 


Pig.  1- Leapfrog  procedure  for  sdvsnclng  the  w«ve  function  <l> 


particular  phyiical  situation,  litis  second  card  can  contain  the  actual  physical  parameters 
or  a set  of  dimensionless  parameters. 


First  Input  Card 

The  parameters  read  from  the  first  card  are  listed  in  Table  1.  A description  of  each 
of  these  parameters  is  as  follows: 


Table  1 —Parameters  Specified  by  the 
First  Input  Card 


Columns 

Name 

Format 

Columns 

Name 

Format 

1-5 

PHIMXX 

F5.0 

38-40 

NPM 

15 

8-10 

ROCULT 

F5.0 

41-45 

NBM 

15 

11-15 

HXY 

F5.0 

48-50 

NPLOT 

15 

16-20 

NXY 

15 

51-55 

NCT 

15 

21-25 

NCW 

15 

56-60 

NRS 

15 

28-30 

NAD 

15 

81-65 

NPUNCH 

15 

31-35 

NM8 

15 

75-80 

NID 

A0 

15 


D-.Li.T/t:.  Dttl.-l. 


irinmai 


lagamtaami 
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PHIMXX.  Thii  Is  the  maximum  allowed  phase  change  in  radian*  for  any  point  in  the 
computational  grid  at  each  i *tep.  It  i*  uaed  to  define  the  newly  computed  » increment* 
HZN  at  each  itep,  where 


HZN  ■ PHIMXX, 

VDiDj 

in  which  fj(x,  y,  «)n nu  i*  the  maximum  value  in  the  computational  grid  of  g3,  given  by 
Eq.  (78).  PHIMXX  i*  nominally  entered  a*  1.0,  If  more  t itep*  are  required,  PHIMXX 
can  be  decreaaed.  In  this  caie  the  a increment  is  tied  to  the  amount  of  heating  in  the 
atmosphere,  becoming  smaller  automatically  as  large  density  changes  take  place  or 
becoming  large  and  efficient  when  near-vacuumlike  propagation  occurs.  If  HZN  exceeds 
0.1  of  the  total  propagation  distance,  the  smaller  of  these  two  * increments  is  used.  If 
HZN  at  any  time  is  less  than  10  “7  times  the  distance  to  be  propagated,  the  program  exits 
and  an  error  message  will  be  printed. 

ROCULT.  This  is  used  when  propagating  uniform  circular  beamshapes  with  an  obscuring 
disk  or  a uniform  rectangular  beamshape.  In  the  former  cose  ROCULT  is  the  ratio  of  the 
occulting  radius  to  the  total  radius.  For  a rectangle,  it  is  the  ratio  of  the  y to  the  x 
dimension.  ROCULT  is  used  only  when  NBM  equals  4 or  6. 

HXY.  This  parameter  defines  the  slae  of  the  computational  grid  relative  to  the  nperturo 
radius  by 


Ax  - Ay  ■ HXY 

where  Ax  and  Ay  are  the  sizes  of  Individual  computational  cells,  which  start  out  squoro. 
Depending  on  the  beamshape,  values  between  0.1  and  0.3  are  typical. 


NXY.  This  is  the  number  of  individual  computational  cells  along  the  edge  of  the  entire 
computational  grid.  The  FFT  routine  is  more  efficient  when  NXY  is  a power  of  2,  and 
NXY  is  normally  entered  as  64. 

NCW.  This  parameter  permits  CW  propagation  to  be  included  by  allowing  the  summation 
in  Eq.  (72)  to  be  replaced  by  an  Integral  [7].  Before  the  summation  is  replaced,  Eq. 

(72)  can  be  written  in  terms  of  physical  parameters  as 


my  l)*aagj>g'a< 


][]  l‘i'U  - n(i>o  + , y , *1  lu . 

n-l 


This  summation  is  performed  when  NCW  - 0.  When  NCW  - 1,  the  program  is  in  the  CW 
mode,  and  Eq.  (72)  is  replaced  by 


3N(?-l)*aa/>e-««v'B7 

c,u(to 


|4»(*  + x\  y,  *)|a  dx', 
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whore  P is  the  average  power  o f a CW  laser  (P  m Ep/At).  The  integration  is  performed 
using  a simple  trapezoid  rule. 

NAD.  When  NAD  ■ 0,  the  coordinate  system  adaption  is  not  included.  When  NAD  ■ 1, 
it  is  included. 

NMS.  When  NMS  *»  0,  the  midplane  integrations  are  not  used.  When  NMS  ■ 1,  they  are 
used. 

NPM.  When  NPM  ■ -1,  the  second  data  card  contains  physical  parameters.  When  NPM  * 
+1,  the  second  card  contains  dimensionless  parameters. 

NBM.  This  parameter  selects  one  of  the  five  beamshapes  available  within  the  program: 

NBM  ■ 0- infinite  Gaussian,  with  WIDTH  (a  parameter  read  from  the  second  input 
card)  being  the  Intensity  radius; 

NBM  ■ 1 - Truncated  Gaussian,  with  WIDTH  being  the  e“*  intensity  radius,  trun- 
cated at  vflf  X WlDT^i  c'r  c‘a  intensity  radius; 

NBM  * 2 -Uniform  circular  ape 'tor, , with  WIDTH  being  the  actual  aperture  radius; 

NBM  ■ 3 - Uniform  square  aperture,  with  WIDTH  being  the  dimension  from  the 
center  of  the  square  to  the  edge  (half-side  dimension)  in  the  x or  y 
direction; 

NBM  ■ 4 - Uniform  circular  aperture  and  an  occulting  disk,  with  WIDTH  being  the 
total  aperture  radius  and,  as  stated  previously,  with  ROCIJLT  being  the 
ratio  giving  the  occulting  disk  radius; 

NBM  ■ 6 - Uniform  rectangular  aperture,  with  WIDTH  being  the  half-side  x dimen- 
sion and  ROCULT  being  the  ratio  giving  the  y dimension. 

NPLOT.  This  determines  the  type  and  the  number  of  plots  given  in  the  output: 

NPLOT  - 0-No  plots; 

NPLOT  ■ 1 - Final  contour  plot  only; 

NPLOT  ■ 2 -Final  contour  plot  plus  a plot  of  average  intensity  and  peak  intensity 
versus  s; 

NPLOT  ■ 3 - Preceding  plots  plus  a plot  of  flux  und  arou  versus  Imullanco; 

NPLOT  ■ 4 - Preceding  plots  plus  u contour  plot  of  aperture  intensity; 

NPLOT  ■ B - Preceding  plots  plus  Fourier-trensform  contour  plots  of  aperture  and 
final  Intensity  distributions. 

NCT.  This  determines  the  contour  levels  used  in  the  contour  plots: 

NOT  * 0 — Contour  plots  use  contour  levels  with  10%  Increments; 

NCT  ■ 1 — Contour  plots  use  3-dB  contours  (0.B'\  n • 1,  2,  ....  10). 

NR8,  When  NR8  ■ 1,  the  final  contour  plot  is  corrected  and  standardized  according  to 
an  internal  criterion,  to  remove  the  effects  of  different  amounts  of  coordinate  system 
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adaption  in  the  x and  y directions.  When  NRS  * 0,  this  plot  can  appear  with  nonuni- 
form axes. 

NPUNCH,  This  determines  whether  there  is  a punchod-curd  output: 

NPUNCH  ■ 0 — No  punched-card  output; 

NPUNCH  - 1— Punched-card  output  for  later  data  processing, 

N1D.  Up  to  six  characters  can  be  used  to  identify  a run  or  a series  of  runs  on  both  the 
printed  and  punched  output. 


Second  Input  Card 

Tire  data  contained  on  the  second  input  card  depend  on  the  value  of  NPM.  If  NPM 
■ -I,  the  physical  parameters  listed  in  Table  2 will  be  road.  A description  of  each  of 
these  parameters  is  as  follows; 

OM.  The  slew  rate  in  radians  per  second. 

HT.  The  interval  between  pulses  in  seconds,  or 
the  reciprocal  of  the  pulse  repetition  frequency 
(PRF).  For  CW  propagation  this  should  be  set  to 
1 second. 

ALPHA.  The  absorption  coefficient  a in  km"1 . 

ALPHAS.  The  scattering  coefficient  in  km  "l . 

ALPHAS  is  used  to  compute  the  total  extinction 
but  is  not  included  in  the  absorption  that  pro- 
duces atmospheric  heating. 

WIDTH.  The  aperture  radius  a in  centimeters. 

1110  particular  definition  is  given  In  the  preceding 
subsection  for  each  value  of  NBM. 

WN.  The  wavenumber  ft  ■ 2it/\  or  2ir/0\,  whore  0 is  the  beam  quality  and  \ is  the  bourn 
wavelength  In  centimeters. 

VO.  The  wind  velocity  vq  in  metors  per  second. 

ENERGY.  Tho  Individual  pulse  energy  E()  in  joules.  For  CW  propagation  ENERGY  is 
the  average  power  in  watts. 

F.  The  focal  length  in  kilometers. 

54P.  The  distance  at  which  tho  calculation  is  to  be  stopped  in  kilometers. 

As  already  shown,  the  propagation  is  u function  of  five  dimensionless  parameters. 
Different  combinations  of  the  eight  physical  parameters,  which  urc  required  to  define 

IN 


Table  2 -Parameters  Specified 
by  the  Second  Input  Card 
When  NPM  - -1 


Columns 

Name 

Format 

1-6 

OM 

F5.0 

6-10 

HT 

F6.0 

11-15 

ALPHA 

F5.0 

16-20 

ALPHAS 

FB.O 

21-30 

WIDTH 

E10.0 

31-40 

WN 

E10.0 

41-50 

VO 

E10.0 

51-60 

ENERGY 

E10.0 

81-70 

F 

E10.0 

71-80 

55  F 

E10.0 
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these  dimensionless  parameters  and  which  lead  to  the  same  values  of  the  dimensionless 
parameters,  will  produce  identical  results.  In  order  that  a unique  physical  situation  be 
specified,  some  physical  quantities  are  also  road  from  the  second  data  card  when  NPM  ■ 
4-1  (Table  3).  They  are  not  used  to  define  the  physical  situation  but  rather  to  assign 
units  to  the  derived  quantities  at  the  end  of  the  calculations.  The  quantities  read  when 
NPM  * 4-1  are: 


Table  3— Quantities  Specified 
by  the  Second  Input  Card 
When  NPM  “ +1 


Columns 

Name 

Format 

1-8 

F 

F5.0 

6-10 

HT 

FB.O 

11-20 

PNA 

B10.0 

21-30 

PNALF 

E10.0 

31-40 

PNK 

G10.0 

41-60 

PNO 

E10.0 

61-60 

PNS 

E10.0 

81-70 

PND 

E10.0 

71-80 

PNZ 

El  0.0 

F.  Focal  length  in  kilometers. 

HT.  Pulse  interval  At  in  seconds  (-1  second  for 
CW). 

PNA.  The  f number  - WIDTH/F. 

PNALF.  Absorption  number,  ALPHA/F. 

PNK.  Fresnel  number,  WN*WlDTHa/F. 

PNO.  Overlap  number,  2v^2’WIDTH/(VO*HT) 
for  an  infinite  and  truncated  Gaussian  beam  and 
2 • WIDTH/(  VO  * HT)  for  all  other  beam  shapes. 

PNS.  Slew  number,  OM  • F/VO. 


PND.  Distortion  number,  3Nk(y-  l)afEp/c?auoAt. 


PNZ.  The  ratio  of  the  distance  at  which  the  calculation  is  to  be  stopped  to  the  focal 
length,  ZF/F. 


Examples  of  Output 

A series  of  multipulse  runs  was  made  varying  the  pulse  spacing  and  energy  so  that 
the  average  power  remained  constant  and  using  a number  of  average  powers.  The  results 
of  these  runs  are  shown  in  Fig.  2 in  the  form  of  power  optimization  curves.  The  CW 
curve  is  included  so  that  the  convergence  of  the  multipulse  curves  to  the  CW  curve,  as 
the  limiting  case  when  pulse  Interval  Is  decreased,  can  be  readily  observed, 

To  test  the  SSPARAMA  code  in  the  CW  mode,  some  comparison  runs  were  made  to 
check  against  some  results  obtained  from  Jan  Herrmann  of  Lincoln  Laboratory,  who 
studied  the  propagation  of  a CW  infinite  Gaussian  with  a e"a  diameter  of  70  cm.  Tho 
absorption  coefficient  was  0.07  km'1,  with  no  scattering.  The  laser  was  twico-dif fraction- 
limited  DF  with  a wavenumber  of  8.B  X 103  cm'1.  Two  cases  were  considered  at  focal 
longths  of  2,  B,  and  10  km.  The  first  case  had  a power  of  10  MW,  a wind  speed  of  250 
m/s,  and  no  slewing.  'Hie  second  case  had  2 MW  power,  a 2-m/s  wind,  and  a 0,02-s'1 
slow.  The  results,  consisting  of  the  area  containing  83%  of  the  focal- plane  power  and  of 
the  peak  intensity  aro  summarized  in  Table  4,  Aro|  and  Lh  compare  these  quantities 
with  those  that  would  hovo  been  obtained  If  thore  were  no  thermal  blooming.  The  results 
for  these  highly  bloomed  eoBes  ngroe  within  about  B%  with  those  of  Herrmann. 


19 


WHITNEY,  MADER,  AND  ULRICH 


Table  4— SSPARAMA  Results  for  the  Propagation  of  a CW  Infinite  Gaussian 
With  a Wavenumber  of  8500  cm"1,  an  e"2  Diameter  of  70  cm,  an 
Abiorption  Coefficient  of  0.07  km"1,  and  No  Scattering 


Focal 

Length 

F 

(km) 

Area  A 
Containing 
63%  of  the 
Focal-Plane 
Power  (cm2) 

Relative  Area 
Arel 

Relative  To 
No  Thermal 
Blooming 

Peak 

Intensity 

Iponk 

(kw/cm2) 

Relative  Peak 
Intensity  lr,.| 
Rotative  To 
No  Thermal 
Blooming 

First  Case:  10  MW  Power,  260-m/s  Wind,  and  No  Slew 

2 

57.6 

20.3 

147 

0.0464 

6 

658 

37.0 

10,3 

0.0261 

10 

3543 

4B.8 

1,33 

0.0184 

Second  Case:  2 MW  Power,  2-m/s  Wind,  and  0,02-s'1  Slew 

2 

64.8 

22.8 

26.8 

0.0422 

5 

474 

26.6 

2.96 

0,0359 

10 

— 

2018 

28.4 

0.496 

0.0341 

Another  exumple  of  S8PARAMA  output  is  illustrated  in  Pig.  3,  numely,  the  final 
contour  plot  for  the  6-km  run  from  the  first  case  with  10%  contour  levels,  The  com- 
plete printed  output  from  S8PARAMA  is  included  in  Figs.  4a  through  4c. 
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Fig,  3 -Contour  plot  with  10%  contour  levels  for  tho 
0-km  run  from  tho  flr«t  ease  In  Table  <1  (PNAi.F  * 0,350, 

PNK  - 10,400,  PNO  - 0.003,  PN8  » 0.000,  PND  - 
80.000) 

Figure  4a,  the  first  page  of  printed  output,  is  almost  self-explanatory.  Both  dimen- 
sionless and  physical  parameters  are  listed!  one  is  computed  from  the  other,  depending  on 
which  was  entered.  The  program  options  indicate  the  mode,  either  CW  or  MP  and  the 
beamshape  etc.  The  results  summary  in  Fig.  4a  includes  the  final  value  of  the  energy 
conservation  integral,  Eq.  (2).  This  quantity,  which  is  ideally  equal  to  1,  gives  a quick 
check  on  the  validity  of  the  numerical  calculations.  One  factor  that  limits  the  accuracy 
is  the  use  of  a finite  mesh  size.  As  this  mesh  is  made  finer,  the  intensity  distribution  gets 
closer  to  tho  mesh  boundaries,  and  numerical  errors  may  enter  through  diffraction  and 
the  use  of  a discrete  Fourler-transform  routine  as  energy  is  reflected  off  the  boundary, 

To  avoid  this  reflection,  tho  outermost  boundary  of  tho  computational  grid  is  set  to  zero 
and  the  next  outermost  boundary  is  set  to  one  half  its  value  at  each  t step.  Thus  the 
sum  over  normalized  intensity  gives  an  Indication  of  how  much  energy  was  lost  due  to 
boundary-value  problems. 

The  area  that  is  given  in  Fig.  4a  is  the  area  containing  exactly  0.63  of  the  total  flux 
obtained  by  linear  Interpolation  between  adjacent  flux  fractional  areas.  This  area  will 
include  contributions  from  several  peaks  as  the  intensity  pattern  breaks  up  under  severe 
blooming  conditions,  so  its  meaning  may  ulso  require  a suitable  interpretation  of  the 
Intensity  contour  map.  In  addition  the  relative  area  and  maximum  intensity  are  calcu- 
lated relative  to  the  focal  area  and  Intensity  of  a vacuum-propnguted  infinite  Gaussian 
whose  e'1  diameter  is  equal  to  the  value  of  WIDTH  regardless  of  the  beamshape  being 
propagated. 
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Figure  4b,  the  page  containing  numerical  data,  begimt  with  a list  of  Internally  com- 
puted quantities  that  relate  to  the  problems  of  air  breakdown  and  t-cubed  self-blooming. 
They  are  printed  only  for  possible  future  data  analysis.  Assuming  the  breakdown 
Intensity  at  10.6  /urn  is  3 X 10®  W/cma  and  that  this  is  inversely  proportional  to  wave- 
length squared,  the  following  quantities  are  computed  as  a function  of  range:  the  mini- 
mum area  required  for  breakdown,  the  ratio  of  this  minimum  area  to  the  vacuum  area, 
the  maximum  pulselength  before  breakdown  occurs,  the  critical  power,  the  saturation 
time,  the  Intensity  produced  by  the  critical  power  propagating  in  a vacuum,  und  factors 
accounting  for  turbulence  with  values  of  of  10  15  and  10"11.  This  is  followed  by  an 
x und  y slice  through  the  aperture  to  check  the  initial  beumshape. 

The  quantities,  Including  the  values  of  MSN  in  «//<a2  units,  relating  to  the  coordi- 
nate system  adaption  are  printed  at  each  z step.  The  headings  D,  Dl,  D2,  ALPHA1, 
ALPHA2,  BETA1,  DALPH1,  DALPH2,  DBET1,  and  XCEN  correspond  to  P,  P\,  D2,  «t, 
tvg,  p,  Aoq,  Ac»2,  A/3,  und  X used  in  the  second  section  of  this  report.  Also  Included  is 
EPSMX,  the  maximum  value  of  the  summation  given  in  Eq.  (72);  PHIMX,  the  maximum 
value  of  the  positive  phase  change  applied  to  i p to  obtain  >l>;  and  PARM,  the  number  of 
pulses,  for  the  MP  mode,  that  occur  in  a computational  cell. 

Figure  4c,  the  output  data,  lists  In  the  top  portion  the  area,  flux,  the  area  fruction, 
und  flux  fruction  contained  within  each  contour  level.  From  these  data  the  03%  urea  is 
interpolated.  Tills  is  followed  In  the  middle  portion  by  the  z locations  of  the  maximum 
of  the  uverugo  and  peak  intensities,  the  minimum  03%  area,  and  the  minimum  z step  that 
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Fig.  4b— Second  page  of  the  output,  containing  numerical  data 


have  occurred  during  beam  propagation.  Then  in  the  lower  portion  the  peak  and  average 
intensities,  the  03%  area  and  the  location  of  the  peak  intensity  in  centimeters  are  listed 
at  each  2 step. 


Summary  of  Program  Structure 

When  the  half-step  integrations  are  used,  the  solution  is  advanced  twice  before  the 
information  at  each  2 step  is  3tored.  This  can  be  seen  from  the  flow  chart  of 
SSPARAMA  (Fig.  5).  Thus,  when  NMS  ■ 1,  the  program  actually  used  twice  the  number 
of  2 steps  that  are  printed  and  Included  points  approximately  midway  between  those 
listed. 

The  structure  of  the  code  SSPARAMA  is  explained  below  and  summarized  in  the 
flow  chart  in  Fig.  6. 

• The  call  to  subroutine  START  cuusos  the  Input  data  to  be  read,  The  real  purt  of 
the  64-by-64  army  ip  is  defined  uecordlng  to  the  beamshape  specified.  Initially  the  phase 
of  this  array  is  zero. 
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Pig.  4o -Third  page  of  ths  output,  containing  the 
remaining  numerical  data 


• The  initialization  procedure  continues  with  the  call  to  INTENS,  where  the 
aperture  intensity  is  computed  at  each  mesh  point. 

• The  call  to  DENS  computes  the  quantity  g(x,  y,z)  given  in  Eq.  (63)  and  then 
applies  the  phase  change  given  by  Eq.  (62)  which  converts  p to  4',  The  first  z increment 
is  also  computed. 

• The  main  program  loop  begins  here  with  a call  to  OUTPUT  to  store  various 
values  until  the  calculations  are  completed. 

• The  call  to  ADVANCE  applies  the  Fourier  transform  of  Eq.  (67)  and  then  the 
phase  change  of  Eq.  (68).  The  array  is  Fourier-transformed  back  to  yield  4>(z  + Az), 

• The  intensity  is  computed  with  the  call  to  INTENS,  and  the  boundary  values  of 
the  array  are  tapered  to  zero. 

• The  call  to  DENS  now  Includes  a call  to  VTEANS,  by  which  the  phase  change 
of  Eq.  (62)  is  reversed,  converting  <f>  back  to  i p.  The  quantities  {aj,  ag,  0 } and 
{Acq,  Aag,  A0}  are  found  in  VTRANS,  and  the  values  of  D\  and  Di  are  updated. 

After  the  return  to  DENS,  Eq.  (63)  is  solved  and  the  phase  change  of  Eq.  (62)  is  reap- 
plied, converting  p back  to  <1>  in  preparation  for  the  next  call  to  ADVANCE. 
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Flu.  5— Summary  of  the  code  SBPARAMA 
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• Now  that  one  cycle  df  propagating  the  solution  is  completed,  the  code  checks  if 
z final  has  been  reached  and  if  the  half-step  integrations  are  to  be  performed  as  outlined 
in  the  section  titled  Numerical  Procedures. 

• When  z final  has  been  reached  or  the  time  limit  of  execution  is  near,  the  last 
call  to  OUTPUT  prints  the  results  and  ends  this  run. 

The  Appendix  contains  a complete  listing  of  the  code  with  copious  comments 
included. 
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Listing  of  Cods  and  Commsnta 


PROGRAM  SSPAfl AMA 


COMMON  /«T9V/  Al(6A*i.AI*  A2(6<,,6Al*  I Tt-.NS  I A A * AM  I 
COMMON  /AAA/  ( MSI 6M,(,M> • EPSO I 6 A ,6M ) *AOUT 111*991 lUUUTI 9* 1C  I . 

A INI  2 » AM  I i DAL  Pul  ( 21  *DALPM/(2  I *.»n.Tl  (2)  *ALPlUO  ( 2 1 • ALPnZO  ( 2 ) * 

Bt  T 1 i ( 2 I * 010121*  021(21*  MDl.it?!*  KU2TM2I*  SrtTO 101 2 I * XC.LNOI2I 
COMMON  /BHB/  TtNSlfaN.bM) * 0116MI*  «i2  ( 1 • PuASHIbM)*  PriAft  2 t 6M  I » 

CONMINI 10  1 *M2(31  *SV1 (6M  I *SV2<  6M I * PARK  l hi; ) 

COMMON  /SINC.LS/  K * PDA*  PMALT * PKN*  PMO*  PNK*  PNO*  PN2.  HX  • HV, 
HZ*  2.  2 2 • 2 2 T- » 2NM*  ZHNALi  XZtkO*  YZT.RO*  \li)tn,  «LPmA,  v.n, 
V.IUTt  OMI'T*  dT»  INf.Kt.'i*  ALPhA(  * CS*  .'UiFK/l*  OAMMA*  l,’.»  UK* 

L JT<  J * RHT  t BOUT*  OAI<rA.  2*  Ts,  TP'.ILM  • A‘,2  *PO<  ♦ S l « TCOM  1 * TC0R2  • 
21*  K 1 6 3 M X * 2Z«  R1MXMX*  / i * APRS  * /.I,,  h/.vii*  OKARL*-*  TLNSMX* 
t.X*  l>llt"X,  fcPft.MX*  1 WWNX  * Dl.MX « Ml*  t'0|*'AX«  VTERH*  PH  ll-'XX  *M2NMS  * 
PI*  I "A  X i Jf.AX*  NX*  MY*  NAD*  MX/  * NY, 2,  NX  Y • NXD 1 M * NY01M*  NPT  * 

I PLOT,  Ml  TIH.M>1UF*NX  "»NY."*MM:.*M-LAl1*l.«l.  1 *1)2  *P1  «P/  * SR  Til  l * SKTUZ  * 
SliHUl  2 *XCLN*lLAbT  *S'.'|YT  8 *PNl.' lu  CU'.O  *OC<'N  *1  i)  I *nCZ  10*MCZ/O  *IH.2  1\  t 
HLZ2N*HCZ 12.ALPH1  *ALPn/*Bl.T  1 .CONI  t CON/ *ii2m«H2.‘> . L XQ , tXN*  ■■  T 1 * ,•  T2 
COMMON  /OUTS/  MB'1!  * SC  l.T  AC  *Nl*S » TIP-’’  * NO-  *NI:  X I T * NPLO 1 .Nt’ONCH 

t l!«(l  »**,'*  !•  U K »«**««»»  B»  II  CH**<»*«*Ut* 

CUMPU.X  Al*  A2 
LOO  I CAL  LS 

DATA  Iff*  Mi-  Oil  « 1 ) * (HlH<AC.«<i*  1‘>MI  * I f.AM"A  * 1 . A I * ( r-  T J*  I t 0>, -7  ) « 

I U:Tk.*1  *•'  I'-S)  * I P I ■ 3 * 1 A ) BN?!, '•  I * (Bl)l»i,lf  I'l  * II  J T*  J*1  .A  I -X) 

DATA  I NXD  | M«  6M  ) • I NYU  I '•  * AM  ) * (l|/*/'U 
BANK  * I ■ ) */'>9SI/ 

INPUT  A Mi'  1 Ml  1 I AL  I / A T I <i*' 

T:>1AI,T*T  | Ml  II  T T ( i ii  ■ K V y | 
l '.».I  Al  Cl  « 

I f'N  T I Mill 

•I  Al.  I 1 .\|.  IIIM 

*-■*  X I T * : 

N I Tl  ,{«■ 

I i'l  M * 
t i 1 *i 
//,’  = ■* 
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IT  II,!"|  .ID.  I l/«l 


(.All  INTI  M'i  ( Al  t ,1  Al  'I  . I 

t Al  L IINMA1,  Al,  //l,  1,  1,  ,1,1!'.) 

H/r-.l'  ... 
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1 h roNTINui 
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( MORI  V Al. i 1 1 Iik  i.i.TI  ( i'kKIi’I'T 

c. 

’ C Al  l.  OI'TT'IIT  I • T Al  M , ) 

C H-  T !-■-«(  Ml  ’’A I N I Ml.  I!,  IIS*.  Thai,  .1  1 1 : | In.\i  Mm  li.l  i I 
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1 TNOw*T|MFi.t  F n I 
0T«TLAST-TNOW 
TLAST  *TNO 

|F(  3*UT.LF«TN0NI  GO  TO  H 
PRINT  2?*? 

22  FORMAT ( //23X25U*#*  T | Mr  AlUjRT  AT  / I K"  I » I ] ■ t (.  . 1 A 3il»  » » / / ) 

00  TO  13 
S CONTI  huh 

c 

C ADVANCE  FROM  ll  TO  U+K<lt  C al(  uL<.  T I -0  r,i.„  a l-Lllui'tS  In  a 

c 

AO  CALL  A0VNO  (Al,l) 

CALL  INTFNSIAl,  .FAur.l-  . ) 

CALL  OtNSIAl.  A?,  ZZl*  It  |?t  .TRiu-,1 

IF  U .GF.  ZF1NALI  00  TO  13  C 

c 

c REPEAT  IF  HALF-STiP  INTI. (THAT  |(  f>  I',  INU.o:  l.,- 

C 

IF  I NMS  .LU.  ■ I C»0  TO  Aft 
CALL  Ai>VNCt(A2t2) 

CALL  I NTtNft I Ai!  * .FALSE.  I 

CALL  0 K N S I A 2 1 Alt  Wit  it  It  .tRUl.l 

IF  U »GF  • 2HNALI  GO  TO  I ‘i 


A3  CONTINUE 

00  TO  1A 
C 

C # » t * » » « » » » » » # # «■ » * P » tt  * * » II  > » 11  d . * ll  II  ::  i i I1  <;  i:  . 1 r i,  a « it  >'  t l.  n i>  c <i  il  #i.  * It  * « » « » n « * » * * » « » » it  ♦ » « » 

c 

c SET  NEXll  rouAL  1 FOR  I'MlIATOi-r.  tXlTS 

C 

13  NtXIT'l 

l ft  CONTINUE 

C 

C EXECUTE  ALL  OUTPUT 
C 

CALL  OUTPUT  1 1 TWOF  • I 
PRINT  16 


1ft  FORMAT  I lHl I 
17  CALL  fiTOPPLOT 

C PRINT  RUN  t l"F  I CP  f I Vl  I t 

TWUN«(T.START-t  | : .ft.  L t.  F t U ’ ,r  ) l/ij.  i 
PRINT  1 H t TRI  If; 

1A  FORMAT  I //ilhl  111*1  •»  kl'N  II  l*»tl6tFt*  I N.  i T , » I 

STOP 
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C 

c 

c 

c 

c 

r 
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NPUNCH  * ,i  « NO  IH/NCHFD  OUTPUT 

Kill)  IS  *.N  ll>  NliMUhH  KlAl)  IN  »t>  FOHNAl 


POR  bAllSSIANS  - WIDTH  » 1 / K INTt.NMfY  HAnlllS 
poll  UNIFORM  AHI.WTURfcS  - V'ID'H  • ilAUluft  OK  HALF 


Mn,  dimnmoN 


HEAD  la  PHlN.XXaROCUt  T,HXY,NXYaNO.NAI>.F  l!- »M>. '•  *MK- T'lPLvlT  t 
1 NCT.NRS.NMuNCH.NIO 
\ POHMATI  )Fb.O*tv  I % » Vx  *Ah  ) 

HX'HXY 

HY«HXY 
NX'NXV 
NY  *NXY 

|F  iNX.OT«NXIi|MI  N X ■ N X I • I M 
IF  I NY  »fiT  *NY|>  I M | NY»NYU|M 
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C 

C. 
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c 
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c 

c 
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r 
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F «F»  1 «E*S 
*F-ZF*l.r*S 
ALPHA»ALPHA«l.t-5 
ALPMAS*AlPHA<,«l.t“Y 
Vr,DT*V0«Ht 
ON'.DT  ■OMtMf 
PNA»W  IfJTH/F 

P.NALF  ■ ( AlPHA+ALPHAS  | *F 
PNK«WN*Wl0TH*»i|/F 
PNCW.utWlDTM/V'HlI 

pn5*omI)T*f/v;>dt 

KNFFl.O«CS»CS*WinTH/(Hr  F MAC  »».?*(  <|A'-MA- UP)  I 

PNni'iENFRCiY  / (LhFFU)»winTh**2*t  tj) 

(tCON*>* > « *HhFHAC*  t i»Aa  vah  ,)*wN»»a*t Hkl<(iY« alpha*  1 it *7/1  CS»C.4) 

PNl)»  3 • *HF.KH AC  »WM»  ICiAI'..AA»I  I*  ALPHA»K»M.hlUiY*  1 .(.  + 7/ 

1 ( f.!i»C'i*i,.,tPTH»V  0T) 

PN2  "itF  / F 

2A  IFINCw.K).ll  HT»l. 

C 

£ . # . . »:  ##.  II  It  " . . (I.  1-  . it  A ll  |l  . t»  UK  It  hilt:  *'  . *1  -i  i:  ii  If  k « w i tl  t‘  « il  If  * tl  II  it  III  ll  If  I'  |i  . It  . tf  « « |!  II  . i|  . * . 

C. 

C.  PH  I NT  INPUT  PAWA'-'FTFH'. 

C 

7 0 co, 'it  WI- 

PE «em:r(‘iy/i..i  i, 

P V 1 *V  ■/  I'll'. 

PALI  S«ALPHA«i*l  il.t) 

pali»alpha*i  «i-.i 

P\"»W  1 1)  I P / I P,'  * 

PE.F/l.I.A 

r 

p M X T m»M.< 
c.  1 r i’iiJ  VA  f I 7"F  t AA  . / / | 

PM  I NT  LA. 

MM  Fi'IU'A  M 1 ' XF  7il*  * * "P’mKTS'  1 1 Put  , A t A •>*«/// 

I A x « P I .1  V,  |i4'j  I '..'i  i-A.fAr  I Tl  ....  ■ »hx  » >•  ft  IV  1 1 A|.  »•/  .. ..  .1.  T t •*;>*' 

JMiX'N  t,  :i  1 r A t I’A  ■'A"!  T i #:*.'/  I 

P w 1 1 1 T : * l". A (PS.MX  * 

1 P X'  A I I .iV-i.r  feny* 

,'P‘. ■;  liV»  < 

• P V .f.v  I 
AM',  « , ' I 
*■  i > ■ . r » i ■ i i 

i> p 1 . i . i‘i  • 

N'  T . 

Mi'Al  E ' 

1 i r ■ftr'TAt  t !•  tu\  "A  »F  I ■ .p « r x 1 1 1 1 <a  I ■ m ■ i ' i i ■ ■ > i « | x<»'iiu  ‘i  ■ . / i / 

1 m A(»i  ii|,'.  [ i-  1 1-  i . , i ) ax  | ft  1;/. l.  piii  I 1 / * : l m 1 1 • i,  I 1 .-xu  • ii  iV  « I 

1 xmp':k  «l  I ■ . / 1 1 x * . 1 1 ...  l | /<. i.  | «ii  . r'XAi'S.v  « i ■■  » / 

<1  XAiif.tl  «l  1 1 • i 1 J Xo"V  (’/  I >1  ..  .1  •»At'NY  "I’m/ 

ill  M.iif,:,  >l  l.,/*  x 1 An  1 ’l  iiA  c..u . / 1 i ( I >1  |:  Pi/ 

I > hi  if,’  ■ .1  1 1 • / • 'ik  1 ,u ,,i  „ .i  u j ) »i  1 1 ■ i / 

PI  X 1. 1 fN  / " i 1 . a . 1 A x 7 1 I II  i ' ' ) t|l  « I i / 

/■>ax  ii, i T l M i l "1  j ■ ,0  . / 

A 'M  x I lll'A)  PH  A S I I /K  'I  I "‘  I .P.//I 

P If  I'-ll  |M7 

I Ii  T P'A."A  T I A*  HPIHMiHA'1  ' l 1 T I1  'I  / I 

I I I In  ’ , r tin  I I'K  1 f, T | Ml 

II  I NO",  I lf.  II  PIIIM  |p  7 

I 66  F 1 W’.'A  I l ,» >, x x. ii  '"I”  ‘'X.'m:  p I 
I t,y  K-'W-  A | | p>  v 1 1"  1 id  1 ‘ i x 1 1 ■ c I 

li  i >ii  r . i i. , . i piii*  r iam 

1 1 I ".l  .'ill  i"< IM  I a ■' 

1 1 l ’;i  -P  * I ilt/  1 P’f  INI  I ’i 

1 1 I P 1 1 u,  I I P.<  1 A T 1 ■- 1 

II  I I.  I I . N I PK  IM  1 PM  tl"’l  I|  I 
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IF  (NAD.EO.Q)  PRINT  142 
IF  (NA|), 8(1.1)  PRINT  144 
IF  (NMS.EO.(')  PRINT  IS* 

IF  INMS.FQ.l)  PRINT  144 
IF  (NPUNCH, EQ.p)  PRINT  146 
IF  I NPUNCH.  F.Oil  ) PRINT  1ST 
PRINT  148)  NPIOT 
IF  (NCT.80.0l  PRINT  162 
IF  (NCT.Nf.1')  PRINT  161 
IF  (NRS.80.U)  PRINT  16* 

IF  tNKS.NE.tl)  PRINT  166 
U8  FORMAT  ( 20X9HHEAMSHAPE  4X«  INF  IN  I TF  CiAUS.SlAN*)  <">„ 

1*9  FORMAT  ( 20X9HI<lAMMIAPr4X*TRllNf  ATI  IT  CiAHSS  I AN*  ) 

140  FORMAT  I 2l',XRH|,FAMSMAPk‘ 4X*UNII'liHM  C|MCU»I 

141  FORMAT!  20XPH|)t'APSHAPtt)X»UNlFORM  SOUAPl.*) 

142  FORMAT  I ? 1XHMADAPT ION4X7HNO I 

ll'J  FORMAT  I 2 1 XRHAHAPT  I0N4X  JMYrS ) 

14*  FORMAT  I RX2l!H(ALF»STt  P I NT).  OR  AT  I oN4X?HNo I 

144  FORMA  T I H XFl  IlFIALF  “ AT!  P 1 NT  I'tikAT  I -N4X  '1MY I 5 I 

146  FORMAT  ( K'Xl  VMPUMCIII:  |i  LaWD  wTP0T4X?HM)  I 

147  FORMAT  I KXlVHPU.MClin-  CAPO  OPTPUTHxlFiYI  S ) 

148  FORMAT! 1 AX  1 4HNUl.'HF.  H OF  PL0TS4X I i I 

162  FORMAT  I 1 lXIHtIUlW  LtVKL  CONToOi<S6X.'HNO  I 

163  FORMAT!  1 1X1HMU1W  Lt.Vt.U  l ON  TuUh.VjX  JIIYI  >, ) 

16*  FORMAT  ( SX26Hl<l,SCALt  FINAL  CiA'TOllK  Pl.!T6X2M'.OI 

164  FORMAT  I ‘JX26HNI  SC.ALF  FINAL  ( IIMTOIIR  PI.0T4X1IIYI  6 ) 

168  FORMAT  ( FOX  RUM  AMSI'APt,  '.IX  * UNI  Fuki  i Cl<Cll  - l.CCULTl  P RADIUS  • "FP.il 
C 

(;  ft  ft  « ft  ft  ft  ft  ft  ft  ft  ft  ft  • # ft  4 u ft  i'  ■ ft  « II  kl  It  >1  k e V ftu  i>  ■:  <i  u II ..  t ti  / * <m>  « ti  rii«  Immi  n # k NX  t « no  * « k U ft  « U k t «i 

c 

C REF  I NF  AND  STOWt  COMPARATIVE  PHYSICAL  DATA 
C 

C 1 3 ■ 1 ,11/1,11 
R7«Y,llitFt«Wl!'TII 

'•'l  * t 7 . i *P  I /lr.  N I * I . r* 
op  I *:*!>  I *1  1.  .6/".'L)*»i> 

RMT  * I . .'/HT 

POl  I T * F N FKfi  Y*  Rll  T * F J T Y J 
DAR) A.HXXHV 

f At  a ',1)P  f I A|  PHA«F'»F  *1  F,l  I'l  Y“  • 4 / IS  , .’!  «f.  » I I'  | l I ' 1 M ' • I I i ? ) | 

C !P*.t»  r*  !“  SWH  M '•■«*/  I ALPHA* I »l  " I • I iH  Y 1 I 

C PCK"  I l.PI  | »P  >1  / I I , '•'  * l . I - Al  Hi  H <H  »l  II  l | T » IP -I  « I - | «,■ 

/I  | '|Al  . /!•  »C.  TM  , VMS'. Pis 
//F  »PN / /PMC 
TPIIU'.F  ■m 
M* 1 1 r 
C 

DO  4*  1*1 «NP 

,’NP*F|  uA  r ( l-I  ) /I  I.V . A T ([.!■-  1 I ♦ I -3 
/!  TA"/»  •VPF.F 
ft/b."*r  »fTF. 

D./l  T A i*  / 1 TA+  I 1,  «/\'l|«»7 
AV*VV*P 

I x>!  XPI  -P*  A|  F "/h"  I 
Mtl  X * SORT  I ! X I 
A :,*.<;  a /,*/''.’■*  sri  x»»'t 
AbF "SORT  I AS* 1 
AR.AS2/A V 

PC  P ■(  PCR»AS*«"<  1 1/1  X 
TP«C  TP«A.'.,V  I SRI  H*/f 
T Pin  M ■ A'1 1 N 1 ( tPOl  SI  • U'l 
TS»PCR»TP/F1OllT 
Si"1'*  TS»P(.M*t  X/  ( P I " A V I 
IF  ( / ,M  , ■ ... ) U>  TO  10 
TC0R1 « 1 ,0 
rcoR/"  i 
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00  TO  12 

DTuRUl«OT|jRBM  l.UE-lb) 
TCORI«0/OTURH1 
D TUKU2«0TURBF I l.OE-UI 
TCOR2«t)/OTURB2 
CONTINUE 
B0UTI1 • 1 l»Z 
BOUT  I 2 • I 1 ■ AS2 
ROUT ( 3 * I l«AR 
HOUTIA.I  )«TP 
BOUT  1 5 i 1 )«PCR 
ROUT  1 b > I ) *TS 
HOU  T ( 7 » 1 ).S! 

BOUT  (8.1 )«TC0R1 
BOUT ( V » t l«TC0R2 
CONT INUt 


ft.  f~ 

k )i  /!!.  .*  # 


in . 


•m  Oj  r 


i!  ■ . ii  ' ?> 

t 


> ' u . a <i  « u n . ■ • .1  k t.  • h » . <.  li  <>  m 


OK  F I Ml'  INITIAL  AMPLITUDES  AT  API HTUHF 

TRUNCATED  OR  INK  INI TF  GAUSSIAN 
l INI  FORK'  CIMCLK  OR  SWUARF 

TKUNCATtO  OAUBflAN  IS  THUNCA1i.!>  AT  1/1  INTENSITY  KAUllJS 
OR  R I TKtlh.  | ■ 1 .it  U*A 

X2ERO»-(hX«l>*liX/2. 

YZr,RO»-INY-ll*HV/*. 

N X M»  h X * 1 
A!YM.«flY-| 

DX.ARf  A*NX*iVY»I'ARI'A 
00  hh  J«1 *NY 
Y«  I J«1  I *MY  *Y2l.H0 
021 J)«1,0«Y*Y 

00  till  I * 1 *MX 

X « ( |-1  I «MX*XZI  Mil 
if  u «fu«  11  0|  i 1 1 > i .i  -x*x 

r,Bl.«X*X+yt;Y 

N F | M iiAff  SIX':  A VI  I TU"1 

1 1.  1 !;•  >>  ,r,i  , f | on  ro  < ■ • 

hi  Al.»|  XP(»..!,».SSUI 

If  I hi  M,|  I;.  1 . A .‘:1  ..M.u.nT  t.i  I lfl  •'  l " . • l1 
t,  1 T 0 3 V 

ol.F  I M UNI  I O'-  " l li<(  I I it  I’l  I III!  I 

IF  ( A-iH « 0 I i / I nr.  TO  *1.’ 

HI  Al. •!.- 

1 I-  I OAtJ.C.I  » I ..  ) HI  Al.  «.■•>■ 

00  III  3 'll' 

n I F I I F.fM'AHl  AI'|-RTl)W 
IF  INI.M.r.t.l)  on  TO  IA 

IF  ( Al'.B  I H I a 1. 1 • 1 • a AN.  ‘ a A I *»  ( Y I . l i a 1 a ('  • i'  i\l  • Nl 1 1 a I ,J  a i I -il  I'l-  ■ ' • I 

1 r I All'll  X I •111.  1 .,'.OR.Ai«,  | Y I . t . T • 1 .F.iUHY  N'lA'.l  i.  . ) | V A|  afH. 
00  TO  IX" 

oif  INI  ni.rui.Tl  o un  I Four  c I Kelt 
rial 

IF  l am;. iti.  1.1  pi  Ai  *•  .< 

jl  I M.HT  I N'llJI  .1  I .lOiCUI.  T ) lit  A|.“  t.i' 
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C 

c load  initial  array 
c 

J60  CONTlNUr 

Ain  *J)  *CMRl.X(REAL,O.0  I 
64  CONTINUl: 

C 

r I Nl'/  NORMALIZATION  FACTOR 


CALL  INTLNSIA1*  a I AL  M- « I 
RNORMiu  ,(i 
DO  66  J<  1 ,NY 
DO  66  1 - l . NX 
RNl  W ■ RN0RM4  TENS!  I « J I 
66  CONTINUL 

rrnokm*  1,'vsurti  i<Ntm»DARi:  ai 
C 



c 

C |N  I T I A|.  I «!  I'  CONST  AMT. r.  AND  PAKA.'il  TL  NS  USL  D LATFK  If!  MiXOC.RA" 

C 

NUL/F  »NXD  I Nff  Y 
NXY«  ?*NXDIf‘*NYU  I M 


ny;>«ny/i» 

F’Vt  I l*?,i  *P|»(  |-l  l/INXi'HXl 

I LO»nxF* 1 

DO  It'  | <|l.O, MX 

!>V(  I I » ? * i ' * !•  I * I l-t-NX  I / INy'O.X  I 
1r  l'V(  I 1 ■ MV ( I IXI’VI  I I 
DO  7?  J» | ,nY 
MUA'I.  Zl.ll*  • *»  *•  «■>  V < J ) 

DO  I,'  IM,M 
C 

'""'"AI  I/I  A"IM  | Tul'i'  |>.  (nl;.  |(H 

AIM  tJl  ■Rl-'f  ' I"  - I.|  11*01 
A.  I I ,.ll  «A|  ( | ,J) 

I I1  '»  I I f J I * « 

I I I J « l , > |I  or,,,:,  1 1 | m ••  \ i ; i 

f.t  (.ft,  T |\",n 

I I ' 7<i  J.  | » •;  V 

I A I I'M’  I | i .1 1 « i , 

( 

r T'Vvl  X AN1'  v *,|  |(|  1.  r r i'l  :<  J , 
r 

I l*X<»JX// 

I I I v « f i y / ;’ 

i);i  <>,'  I < l »Mi 
AIM  I , I I.AI  I IPX,  I I 
A | \ I , I ) «A  1 1 1 • I "Y  I 

A I I.  II  7r 

i<  1 * i x.  • ■ x < r,.- 
A I " 1 1 1 * 1 . 'I  ] , 
m / "\ < t • i r i ■ 
n/ «n,n 

t 111,  I B",  I 1 , - 1 r ','KTI  I,  ■ I I 

cum;’  - o,  h*  (1,  ♦ 1 i.nr  ( < , n 

A1;1  I I I "l.iKif  I I • "NX  I /Uu.l  ( a I *.I.6 

fV  I / I *- 1 o('l  ( l.-'NYl  /l  i'l, i I ? , I I ,1, 

" ;m  i i 

< Al  I :,l  Till'  I ,.\V  I 'V  1 « '.v,1  • till  I-..M 
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OtF  I NE  CONTOUR  LEVELS 

IP  INCT.NF.OI  00  TO  210 

CTLVL*.9 

no  2oo  1.1.9 

CONMIM I | "CTLVL 
CTLVL«CTLVL-.l 
C0NM1NI  tdUO.Oi 
r,0  TO  2 JO 
CONMINI 1 ) ■0.50 

no  221)  1.2.10 

CONMINI I I "CONMINt 1-11*0.9 

CONTINUE 

UASWlMELEFTlUl 

INITIALIZE  PROGRAM  PARAMETERS 

Z«1>.0 

ZZ.o.o 

HCZ1N.J.0 

IICZ2N»i.'.0 

2NM.li. 

n»i. 

n?«i.o 

r>ORTn«,V:MT(H.  i 

no  BO  U]  .2 

mm  1 1*  i 

02ni  1 1 «i .; 

NOl'M  Mail. 
pn?i'<  li.i.;' 

Al. Pill  01  I I ,t>*T*Mlt 
AlPTOM  t I • I'j.PNI', 

HKTli  I I )>  . 

SRTD1-  I 11*1.1 
XCI  V I II  •'  .' 

'.-AL.Pi  . I ( | l«P.  ■ 

DALPuZ  I I l«l  .0 
miETil  | 

CHNT  |M)|- 
a m t n i • i ... 
r. n t r. ;j « t ,e 
RI.HU12.  1.0 
NPl  A Ci a -i 


PUT  I A|  I T I A L U.1I  *;.M  TV  . | 1 IT  1 1 A I T I ; I:  Al  API  i.TlRi 


IT  tNPLOT.l.T.L)  TV'  TO  2 AO 
CALL.  NV! PITOLI1  ."*«•■  .^.lfc.JnZ  *.  .i.JI 
CALL  Nu.iPTRIi  .Hb.U.  . . ] A . t . : • 1 . A I 'T  H . 0 I 
CALI  f.Y'inf.LI  l.Vall..  » . 1 A » All  ►.  M ..,p,(i! 

CALI  LAPFIIC,  .a.oi 
C Al  L I ' L n T I 2 a 5 . .■  .-5) 

XCTNU  R ■ 5 . ** 

CALI  .'VMikil  ( XC  T NTT.  R . 5 . i; . J a 1 A . 1 - 1 I 

XI>'.lN*X/l.mT.p|UTH 

VN1M.V2I  ROai.'ll  Til 

CALI.  TOI’CHiRATI  TLflN./.'Xl  I .NYP  I r .MX  il.V  . 1 . , , ,(.11.1  < 
I .'MIN  .MX  .AHT'b,  1 ,AHX  CM.  + 'i.YA’IN.hY.hlll  (i  . 1 «/.|  IV  CM, LI 
Ri.Tunr, 
rvn 
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SUBROUTINE  ADVNCEl A»NS) 

THIS  SUBROUTINE  FOURIER  TRANSFORMS  PH|(X.Y,Z)  TO  PH  1 ( K 1 • K.2  « Z ) * 

AND  ADVANCES  THE  SOLUTION  t)Y  APPLYING  THE  PHASt  CHANGE 

PHI  ( K 1 ,K2»Z4HZI*PHl  < K 1 »K2,£  l»EXP('.;.Mt*<Kl*«2«HZ/Dl  + IS2«*2Ri.  i/021  I 

AND  THEN  TRANSFORMS  BACK  TO  HEAL  SPACE 


COMMON  /AAA/  tPS I 6A  • 6A 1 »EPSO t PA  ,64  I , AOUT  ( 1 1 • V*/ 1 .BOUT  ( V , 10  I * 

A IN ( 2 ,64  I tDALPMl I 2 I , DAL PH? I 2 ) .OBET1 1 2 I .ALPH10I2! ,ALPh20l2l » 

BET  10  < 2)  • 010(  21*  0201,21*  R010121*  H020(2l*  SRT010121*  XCENOI  2) 

COMMON  /BUB/  TENS164.64),  01(641*  02I6A1*  PHASE  1 (641*  PHARE2I64I* 

CONMINI 10) *M2 (SI *SVl 1641 *SV2 (64 j *PARM1 HOI 

COMMON  /SINGLS/  F*  PNA*  PNALF*  »NK * PNO*  PNS . PNU*  PNZ.  HX.  HY, 

HZ.  Z,  ZZ  »ZZF  * ZNM,  ZFINAL*  XZERO*  YZt'RO*  WIDTH.  ALPHA,  WN. 

VOOT,  OMDT,  HT*  ENERGY,  ALPHAC.  CS*  RbPKAC*  GAMMA  * fTJ,  C'TKt 
EJTKJ.  RHTi  POUT,  DAREA,  '.v2,  TS,  1PULSE,  AS2  »HCH , S 1 , TCUR  1 » K0H2 , 

Zl,  RI63MX,  22 , RtMXMX,  2 3 , APS'H,  /.4t  HZMN,  OKAHEA,  TENSMX, 

EX,  PH  I MX  , EP3MX,  EHRMx,  DGMX  * Rl,  ItOIMAX,  VTERY,  PH  I MXX  .HZNMS  , 14 

PI,  IMAX,  JMAX , NX.  NY*  HAD*  MX?  * NY2 * NXY*  NXDlK*  NYDlf't  NPT * 

I PLOT  ,N  I rEH*NbUE*NXV.*NYM*NMS*MFLA(,*0*in  ,D2*P1.P2  if  N TUl , SHTLV  * 

RSRD1?  iXCEN  *TL  AST  , SURT  8 iPNOf'  ,GCONO,(<CON  nil)  I »HCZ10»MCZ2U»HCZ IN, 

HCZ2N*HCZ  12  *ALPH1  * ALPH?  if- 1 T 1 iCOMl  *CUN2  ,mZi.»hZN *E XO* EXN*  ' T 1 *,■. T 2 

,#*,**»*»» 

COMPLEX  A ( 64 * 64 ) 

c 

C P I NE  PARAMETERS  FOR  PhAFE  TRANSFt'Hi  ATIi.’N 
C 

HZ  1 »C  ON 1 *HZN 
HZ?«CON?*HZN 
201 1 * I HZ  1 +HZO I *RD1 : INS) 

2012*1  HZ  2 + HZ 0 1 *R01 . IMS) 

ZD21*1H21+H?0)*R02'  (MSI 
Z 0 ? ? “(  HZ  ?*H? ()  | *Rl I? ,■(  N!i  I 

1 1 1 1 *0 1 P ( Nf  ) » (I  1 »A|.Pl  l]  :.(M.F  | « ?ull  I«»?*Z!'U*Z.M1  I 

012*01  ( MS  I » I ( 1 • l + ?*(:»ALl'iii  - I ) « ?IM  2 I ««?+/0J  2*  Zi>  1 ? 1 
021*0?.  (MSI*  (II,,  ♦2,.l»A(.Pii2*  IMS  I*  /n?i  I «»?+/.,)/ i « Z(;2  ] ) 

022*02  (NS  I » I ( 1 ,->  + ?,  i:»AL  IN  l?.<(  MM  » /!)??  I • «?  + 2i)Z2*/.0«'2  I 
ROll* l.n/01 1 
RPI2*1.  >/H12 

R:VI*l.o/o?i 

RD2  2*  1 ■ i'/02? 

HCZ1N*w,'i*HZM*IMI)1  1 + 1< t ) 1 ? I 
nCZ2N*'.  ,S*HZM*  | NH21  ♦!<’.? 2 ) 

PERU1. 1 uFOUT  I Roll  * H i ’ ? 1 1 
R SR0S2  ■ SCR T I Ml>  1 2 *R0??  | 

Ilf  Z 1 2* , ,6*H2N*  IRSRl)Sl  + RSR'i:,2  I 
C 

C i.'HFM  VMS  - 1-  * P2*‘. 

r 

Wt  1 *P? * I HCZ  12-11/  1*N<-,N|;M-.|/2«R},|<US?  ) 
w T 2«hCZ  1 2»WT  1 

0 

C RESORT  ARRAY  IF  MX  LT  64 
C 

IE  I NX  « E.  0 » NXO I M | (,()  TO  US'.) 

00  116  Jrl.MY 
On  lib  l«l,NX 

II''  A ( I + ( J»1  I »hX  I * A ( I * „ I 
1 1 !■  4 CONTINUE 


no 
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1 

)• 


r perform  fourier  tramsfok1-  - ru  v.-spact 

c A | L FAs1T0ur-ilAll.il.  V2  . -‘-■V  1 • SV2*  -It  1 F E k R ) 
IF  (NX  .LT.  NX'IlM)  CiO  TO  1 0 


C 

C 

C 

c 


1 0 
c 
c 
c 


1 ? 
c 
c 
c 


CLOT  FOURIER  TRANSFORM  OF  INTENSITY  OlsTRlBOT  IONS  AT 
APF  RTuRC  Ann  Z«F1(|AI. 


CONTINUE 

ap,fly  phase  r. hanoi:  to  auvanci  the.  l»lculat  ions 
00  12  J * 1 1 N Y 

no  12  i*i *nx 

PH  t *R 1* ( ( HC? 10+HC 2 IN  1 *PhASU 1 ( 1 1 * ( HC220+MCZ2N > *PHASt 2 ( J I I 
A( jT)«A(JTI*CKPLX!COSIPHI ) t-SlNIPHl ) ) 

PERFORM  FOURIER  TRANSFORM  “ 70  REAL-SPACE 

CALL  FASTFOUKIAl  Ml'  H2*  ■'•-Vl*  SV2.  )*  IFERRI 


pesort  array  if  nx  i.t  6a 

IF  (NX  ,FQ*  NXUlHI  GO  70  USS 
DO  110  J* 1 t NY 
DO  1 1 fa  1 * 1 tNX 
16  A( I*JI«AI !+(J-ll»NXI 
ll<-e  RETURN 
FNO 


non 
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SUBROUTINF  PFNSIA.  H,  2C00R0.  N01.  ND2.  LP) 

THIS  SUlROuTINt  APPLYS  THfc  pilASt  CHANCit 

I 1-X»»2I/P1+I  1-Y#»2  )/D2-3Nir.AMWA-l  I *K«  »2*E /CS«"7/5uRT  ( 01  "02  I * 

SUMIPHl IX-XP.Y.2  ) )**?. 

THIS  CONVERTS  PS  I (OR  Al  TO  PHI 

#»#  «»#»«»#•*»  »*»#««#»«•»»#»«»»  « * # * « . 
COMMON  /AAA/  EPS! 64,64 ) . E PRO  ( 64  , 64  1 » AOt  'T  II 1 »YV)  tlaOuT  I Vi^V  * 

. A IN  I 2 ,64)  .PALPIH  I 2 ) t|)AuPH2!2  I .IHiLTl  12)  • ALPH 10121  ,ALP)i2  Pi  * I * 

. PE T 1 1 ' I 2 ) * I'l'H  2 I » 020(21  < M01CI2).  HLi?l>  12).  SRTP10I2I.  XCENOI  ? ) 
COMMON  / BH  fl  / TENS  I 64 .64  I » 01  I 64  I • 02164).  PHASE  II 64).  PHASE 2 164). 
. CONMlNI  IP  ) »M2  I 3 I . Svl  I 66  ) »SV2  I 64  ) . PARMI  ti1' ) 

COMMON  /SlNCiLS/  F . PNA . PklALk  » PNK  . PRO*  PN5 . MNP.  PN2  . HX  * HY. 

• HZ.  2,  22.Z2F,  ZN.Mi  /FINAL.  X2LR0,  YZfcWO*  WIDTH.  ALPHA.  WN  t 

■ VOPT  . OMDT . HI*  FNEKOY,  ALPHAC  . CS.  KlFKAC.  MAMMA » ITJ,  CTK. 

. KJTXJ.  RHT  . POUT.  MAKE  A . .,2.  TS.  TPUlst*  ASZ  »PO<»  S I » TO-K 1 1 TC0II2 , 
. Zl.  HI63MX.  22.  KIMXMX.  23.  APMN.  24.  H2MN . OK  ARE  A . te.NI.MX, 

. EX.  PHIMX.  t'PSMX,  IRKMX.  CK.MX.  Ml.  IIP  IMAX  . V T l.  KM  • PH  IMXX  mZNMS . 

• PI.  [WAX,  JMAX , NX,  NY,  NAP,  NX2,  NY2 « NX  Y » NXPlMi  MYUlM,  NPT  « 

. I PLOT  ,Nl  TEM,NHIIF,NXM,NYM,NKSiNFLAC,*P.PI  .02  »P1  »P<  *SKTP1  tSMTUZ  , 

• R SUP  12.XCEN.TLAST.  SORTS  »PNI; .OCONO  .liCuN  »HJ  1 »IIC2  t0.K220.HCZ  IN  . 

. KZ2M.HC2  12.ALPii1  .ALPM2  .BI  T 1 ,L(;N  1 . C0N2  .HZU.HZN .1  XO.kXN.Wl  1 ,oT2 

COMMON  .'OPTS/  HUM » SCI.  F AC  . k|l<S . MV  »(IC!‘.'  * HI  X I T .MPI.0  T .NH 'NCH 
C «#»*#.*«*«  n »»»  * n <i  »««  nil  i:  ««  a L <i  ««  u li  » is  i*  * 
COMPLEX  A I 64. 64  I . 11164.64  I 

logical  i.p 

P I Ml'  NS  I ON  T ) N (l.'l,  EPSL'Il..).  ll/SAVM?) 

C 

C INITIALIZE  2-STI  M PM  FIRST  CALL 
C 

IF  ll.PI  HO  TO  61 
H/Hur.  . 

H/iVWr>».  ,ll 
1 1/  SAVI.  I I I 
1 1/ SAVE  I 21  ■!'. 

P 1«:  .6 
PZuO.P 
P .i  * 1 . : 
i:xm*u. 

r 

' 7 F " / p | n a |.  / x.  » 1 :i  l li » « / 

H/“X---w.  1 »/2l 

I k ( mm:>  . I'.).  . I hZ-'x*.  • s * m/.-ix 

H2"  I f:  I « t ,'•!  “4#Z/k 
li  I V>. fr'."  ) HZ"  I M ».  If.  I 

4 1 / Cl'ORIs  n / ( ' i"  ■ 1 I l 1 1/  5.A  V1  (M‘-l) 

r //  ■=  /(Cm)  / x » '.'iidth11" / 

r 

<•?  zz»zcooi<n 

c 

r 7 MM  ■ 2 / F 

r 

;,vm»//»p  fix 

r / « 2 I x M ) 
c 

/ n/M.M»k’*(  T I I | , 1-6 

(. 

P*/2#22  + I 1 > '-/4'l ) * «;’ 

/ni)l.H/#RO|  I Nil  1 ) 
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D 1 ■ 0 1 0 ( N51  )* { ( 1.0  + 2.U*AlPH10lNt)l  I*ZI>D1  )**2+ZDDl*2DDl  > 
SRTD1«S0RT(D1 I 

VTERM  ■ DISTANCE  BETWEEN  PULSES  / DISTANCE  BETWEEN  GRID  POINTS 

VTERM«2.U*  ( l.r-+|>N$»2NM)/!HX*SRTDl*PNOl 
IF(VTERM.LT.O)  GO  TO  A5 
EPSMX>0,0 

JC-IFIXtl.O/VTERMUl 
DC-FLOATl JC-2) 

IF  (NCW.EO.l)  GO  TO  200 
IF  ( DC ) AB,A9,5U 

EXIT  OPTIONS 

n punch* o 
NEXtT*l 
PRINT  H'U*  7 

FORMAT  I 772X«  * AT  /•*«  E6.A,  » KM,  A DEAD  ZONE  IS  PRESENT  IN  THE  C 
• ALCUI.AT  IOM») 

CALL  OUTPUTI .TRUE.  ) 

Stop 


A A PRINT  101.  2 

101  FORMAT  | II  2X  » AT  2*  » F ft.  A » K'-'t  I HERE  AMI.  MORE  I HAN  lu  PULSES 
. PER  CELL  PRESENT  IN  Tht.  CALCULATION*) 

STOP 

C 

at  MPUNCH*t> 

N f'  X I T c 1 
PRINT  IPS.  I 

103  FORMAT  I //  2X  * A T /«  * F7.A  * KM.  tMI:  CALCOLAT  LU  ncl  IS  SMALL  t M 1H 
.AN  THE  MINIMUM  ALLOWED  VALUL*  1 
CALL  OUTPUT ( .TRUE  « I 
5 I OP 

c: 

£ ft  . # <i  H « <;  ,Ml  * li  u , v it  . * ii  ii  i: .»  « i’  v i-  ..  : i;  n : i . , * . ..  u {.  .t  . » i.  v . , 

C 

C.  Cl  A'  Till  INTI  MS  I 1 Y Al  IvimE  Till  l-KlU  ME  . <UL  I I - •MIL  M 

C I N T t.  wiR  A T E TIM  1 N T I N S I T V , . c i<OR  N 1 > il  t*K  1 1 - Mm  tv/ 

c 

c 

C I ESS  THAN  ONE.  PuLAE  I'fP  till. 

C 

/,il  Il»VT|R" 

I 11  * I I ■*  1 
FI ■ I l 1 “VTERM 
E 2*1.0 "FI 
DO  A .1*1.  NY 
DO  A I » ? » NX 

I E ( I-l  1 1 .r'f  . 1 ) no  To  AA 
EP'Ol  I »jl»|-Pt  ( I »JI 
EE'S  I 1 ,J| 

(,0  T 0 A 

A A I.PMM  I ,J)*E  PS  I I »J) 

l-PG  I I . J)  *C  1*  l Tl  Mr.(  I-l  1 . .1 1+;  PA  I | - I 1 , J|  I » 

. F7*  ( Tf  N!.  ( 1-11 1 , JI+LPS  I 1-  I 1 1 • .J  I I 

IE  I 1 E N.‘,  f I . J ) .(iT,  .i  *,*11  NV-'x  | i.v'i -‘x*  Aim*  M i h,..|a  ,1  w:,  \ 1 . j n 

A CON  TIMM 

C.u  TO  A] 


ONI  To  TVcd  Pill  !-l  P|  R CELL 

llll  HM*?,vi»VTE  H" 

E 1 * . vl-1 1 T E MM 

12*1 .H-l 1 
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5 


DO  5 J«t»NY 
E P S 1 • 0 • 0 


TFM1»u.5»(TLNS(  l*J)+TtNS(2,J>  I 
t PSOI 2 * J)  ( a * Jl 
fc'Hsi2.j)*Fi#Tt:fa+fta*jt;NSi(  i *ji 
£PSWX«AMAX  1 ( tf’S  t 2 t J ) * U'St'AI 
DO  i I - 3 ,NX 

EPS  ! *f  I #1  TtNSI  1-1 1 J1+KI»S(  1-1 1 j|  )+(•-;>«(  I EMI  +EI-  si  ) 

TPN 1 *()•  8*  ( TEMS  ( I - 1 » J ) + TKN5  1 I t J ) I 
EPSOI  I t J)*KP'.(  I , J) 

J • TFNH-KPf.l  J+f  ?*(  TK  ’^c  !-l  , J)  *r  l-l.JH 

CONTINUP*1  l,JI  *0T'  l|,r.r.X«A.MAxlli-M!,yx,tP'; 

GO  TO  hi 


ItJH 


WRL-  THAN  TWO  PoLSUK  I'Ll'  CLU. 

IF  UC  ,GT.  10 1 00  TO  14  6 
uTEf»M«FL0AT  ( JC  I *VTKRM 
F 1 "2*0-0 TERM 
F2«1«u«f1 
no  6 J«1*NY 
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1-2  (8 

EPsom-o.  o 

EPSO (21-0*0  : 

TENQIl  l-TENSI I-1*J)  J 

T E NO < 2 ) ■(  TENS  ( I * Jl.+  I FLOAT!  JC  1-1*0 t»  I ENbl 1-1* J)  t/El-OAi  i JC  ) ;j 

DO  70  JJ-3*JC  a 

FJ-PLOATtJJ“l I /FLOAT ( JC  I f 

TENOI  JJ)-FJ*TENSI  I.JI  + I 1 *0-F  J I * I tNi*  1 1 -1  * J I i 

EPSOI JJ I -FI* ( TENui JJ-1 l+EPaut JJ-1 1 1+F2*  t i En« . JJ-2 I +EC jj-2  I I ^ 

TO  CONTINUE  | 

EPSOI 1 .JI-EPSI I . J)  i 

EPS  I I «J) -FI*  t TENOI  JC  t+EPSu( JC  1 1+F2* ( TENui JC-1 1 ♦EPSOI JC-1 1 I fl 

FPSMX-AMAXKEPSI  1 * Jf  * EPSMXI  1 

DO  (.  I-S.NX  I 

EPSOin-KPStl-ltJ)  I 

TEND  1 1 1 ■ TENS  I 1-1  *J ) J 

EPSOI  21  «F1*(  TENOI  ll  + EPSOI  1)  I + E'2»  (TUNO I JC  1 ♦EPS.)  I JC  ) 1 » 

T t'NO  I 2 ) • I TENS  1 1 i J I + 1 FLOAT  I JC  I - 1 *(J  I » I EN5»  1 1 -1  * J 1 1/HOAII  JC  1 j 

DO  71  JJ-S*JC  I 

FJ-FLOAT(JJ-11/FLOATIJC)  i 

TENOI JJ I -F J*TENS( l *Jl  + ( l»u-FJt » i tiNsi  1-1  u I 1 

EPSOI  JJl-Fl*  I TENOI  JJ-1  ) ♦EPSOI  JJ-1 1 I+E2*!  1 |'Nu  I JJ-2  I ♦ I'.PSUI  JJ-2  1 I | 

71  CONTINUF  ,| 

EPSOI 1 * J I -KPS  1 1 * J I | 

EPS  I I * J I -K 1*  I TEND  ( JC  I +E  PSiC  I JC  I I +l;2*  I I F NM  JC-1 ) ♦tP&y I Jt-1 I I jj 

IF  (Tl'NSIllJ)  *0T*  *l'5»  I EiSSNX  1 ■ Ah.A  A 1 * E * E-v'i  > I * J 1 I f 

6 C0NT1NUI 

no  to  si 

! 

COMPUTE  C.W  INTI.CiHAl.  i 

200  DO  11. 1 J«liNY 

EPSOI 1 t J ) -EPS  I 1 * J ) 

EPS  I 1 iJl  -u.!>*HX»TI:NM  1 * J 1 »;.•  I litrl/  I Vl'iJ  t » I 1 * + PN:>»C.NMI  I 
1 »S(IPTID1I 

DO  111)  I - 2 *NX 

F PSD  I I * J ) -EPS! I * J I i 

EPS!  1 * J ) « F P S I I - 1 i J I 4 *»-t»X»<  Tl  NS  I 1 tJI  + UNM  I - 1 * Jl  I *•*  lulu/  1 

1 I VOOT » I 1 . + PNS-2MM  I ) 

1 »M)P  TIPI)  * 

IE  I f I NS  I 1 * J ) *l>1  *.'  * . S * TEN. -fiX  I IP:.  * ■ A:SA  M 1 1 IS,  -'A  <1  l>  j I I * J 1 1 
IK'  CONTINUE 
C 

£ # 41  <t  # 4MMMH1  41  *»  Jt  # *1  Mil  4 « tt  i*  1M.  I.  IMF  If  « 4Ui  'I  - L It  «t  tl  C |»  ,1  D M !:  4 ♦»  IMMM  >;  i»  V 0 M tt  * # W f U M « * « # * * • »t  * li 
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C 

51  EXO*CXN 

EX«EXP(-PNALF*ZNiMI 

exN«ex 

OCON,tiCONO«EX 
IF  (ID)  CALL  VTRANslA.hOl  I 
HZO*HZN 
HCZ10*FICZ1N 
fICZ20*MCZ2N 
C 

C fALCULATt  THF  Z-INTEGHAL.5  wHtN  N.v.a  • u 
C 

IF  INMS.NF..O)  C.0  TO  52 
HZNMS»HZO 
HZI*C0N1»HZ0 
MZ?*CON?»HZn 

Znii»nZ|*Rr«lP(Mnn 

Zni?»HZ?#«rilMN«U) 

2D2 1 *HZ  1 •R02:.  I Nl'l  I 
ZD?  2*H2  < NO  1 1 

Dl  1 *CUM  VDI  l*(  I l*''*2*v  *ALI*tUD(NHl  I */|)  1 1 1 **2+/n  1 1*201 1 I 
01?*Dlt)INDl  )•!  I 1.0  + ?.v*ALPHl,)(NI'l  I *201 2 I **2+201  ?«.'»)  12  I 
02  I »D?Ut  NDl  1*1 1 1»*>  + 2«**ALPH2i.  IN01  l*20i!l  I <1*2+4021*402  1 > 
02  2*020  ( N01  )« I (1,0  + ?,.  < ALPH20  INtU 1*2022  I **2  + 402  2*41)2 1 I 
Rim*i.vi/M  i 

H01?«l»../r>l2 

MO? \ *1.0/021 

00??4 l,(.'/o?7 

HSH05 1 * SORT  I KOI  1 »R0? 1 1 

MSHDS?*  SWT  I KOI  2»«0, ■> 2 ) 

C 

<:  tor  i z-t‘:nr,i!ALf> 

/* 

KZlon  .Mll/0*(HUU  + Fnl2l 
1ICZ  20*  ) » A *M?.('*  I HD?  1 + RD2?  ) 

DC 4 1 ? *0  , h *lt/(i*  l KSHI’M  ♦ ROM. 2 | 
v;T?*llC/  1 2 

<• 

C C|>".|MTI  MW  / l‘:c  10  I i ,* . T I*'  //:■/«<•■?  I F | I 
r 

? t'/.R*)1 1 11  A"  | M I . , ■<(*'.'  1 , ■ .<ii|Vi  l’M|l'XA/ 

. ( in?*1  ft":*  (I  i""x  + ] 1 -•>.  ill 

II  (ll/X|. I'.T, M/OXI  H2',»M2f‘X 
IF  |n/0  * l.  T . IIZ'HM  I I ('  10  <1  V 
IF  I’l/N  .oT.  //I  -/./“ll/M'-i  ..  I H4!,«/4l  *4  4*n» 

m2  »ii2(>  + H/N 
H2 5 avi:  (fin?  i*fi/ 
r 

rOl.'POTI  tl'F  Tlllll. I 2 - 1 M II.  OllAl  :■ 
f. 

r 00  ONLY  Oil  F | list  CALL 
r 

|F  (l.’>)  FiO  TO  »A 
IF  I MM*.  §f  0»o  I I’l*'  •*> 

H/.  1 *f.i?*il  »HZ 
II4?*C()N?«M/ 

701  1*H/1*HIjU  (NO?) 

ZO)?*H/?  + »ni.  I NO  ? ) 

ZI>?l*Fi2]»l»l'?  -INIJ7I 
Z'>??*ll/?  + HO?i’IND?) 

0 1 1*1)1.' I NO?)  “111..  *2.  * Al  i’ll  1 . ( ;V7  I * ,•  IU  l I * * .'  + 4l 1 1 l *■  1 1 1!  I 
l'l?»l  l"IMI'?l*  I I 1 • . + 2 •..•ALF’ill’  li.l’i.’! ' M>12  I :’«?♦.  I ’l."’,  I’l,'  I 
n?  1 >02  ■ I NO?  )*l(l*  ■-  + /•>  * A I.  I’M,''  |M',’  I*  / 1 1?  I I » *7 + 4 IV  I *1  > , i I 
0 ? ? •D?o  I NO?  I*  I I 1 ♦ ?•>  » A I I’M?,  l I 2 ) , I’?,1  I « I I , l |VI 

poll. 1 ,0/011 
P 0 1 2 * 1 .0/01? 
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R021»l,0/tl21 

RD?2«i.i:/P21» 

RSR[)M»5UHMMUll»RD2l  I 
RSRDS2  »M)R  T(R012»RD??I 
MC2U;>>  . •»  »n/ « ( KU 1 1 1 ? t 

MC2;c;«.),!>»n2t(m)21+Riij'?  i 
HC2  I 2*i: « 4*H2  • I WSUINS 1 *WSKD&2  I 

\l  T 1 *v  ■ i 
'V  T 2 ■ tic  2 1 ? 

riiMPUTC  TIT  IWSI'  < I AKC.I  IN  Til  I L<  r |-  |t|  ,N  fr.K  ,1 

CON  1 1 MIC 
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PHIMX«0,U 

IF  (NMS  . EO*  0|  GO  TO  8U 
00  Sft  1 »NY 
Y«VZtRO+FLOATIJ-n*MY 
00  SS  I »1 *NX 
X«XZER04FLOAT  t l“l I *HX 

0N6w»GC0Nn* | wn«EXO*EP50 I I * J I 4Wf  2»fcXN»EPS I I *JI  I 
PHJ«o.5*Pl*(HCZ10*GH I (♦HCZZOtr^i Jl-GNEw) 
PHI«PH|-X*X#0ALPHl(N0?)-Y*»»nAUPH2iN0?t-A»r>flE.  1 ,ND2) 
B(  I *JI»M  I .J)»CMPLXlCOsiPHl  I i oINiPMl  I I 
IF  I TENS!  I » J)  • L I • u»u5* i I Gu  ,U  5S 
PH|MX«AMAX1(PH|MX* 

89  CONTINUE 
00  TO  AO 
C 

C COMPUTE  THr  PHASE  CHANGE  IN  THIS  LOUP  ,.HtN  NMS  » u 
C 

80  (10  89  J«J*NY 

Y«Y2tR0+FL0AT(J»ll*MY 
DO  HA  I ■ 1 *NX 
X"XZERO*FLOAT< t-lt»HX 

GNE'i‘.*(iCONv  * ( WT  1 *K  XO*  l:  I'SU  I 1 < J 1 4 T 2*t  XMEPM  I i Jl  I 
PH|»L.A»Pl*IHCZU)»Gl  I I |4HC<8o*<i2UI«UHt»'  I 
PH  | » PH | >X*X*DAUPHl  ( NI!J  I ■Y*Y*()ALPM2  I M>? ) »x "HHl  i HN02I 
At  I *JI«AI  I * JIKM'LX  I COS  U>0 1 I . s IN,, in  II  I 
IF  t TENS  ( I . J I .IT.  o»,)A*l  LMSI'X  I 00  i 0 S9 
PHI  MX" AMAX 1 1 PH  I *X . i . 9 *GNI V I 
HA  CONTINUF 

A<'  RETURN 
END 
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SUBROUTINE  INTCNRIA,  LM 
C 

C THIS  SUBROUTINE  TAPERS  THE  BOUNDARIES  OK  THE  COMPUTATIONAL 
C ARID  TO  ZFRO,  FINDS  TmE  iNItNalil  Ai  F.ACh  (mlO  culm 
C AND  PLOTS  THE  CONTOURS  OF  t HE  KOuRltR  i KANaTOMMfcD  ,.,AiMa* 

C 

ci  • I • t • I • • • • • • • • I • • « « > • I • I « < 4 « • I • « < I I 

COMMON  /AAA / L PS ( 64  *6A I i E PRO I A A • b A I » AOU I llltV91  ibOUItViluli 

• A INI  2 iA4 1 • DAL  PM  1 1 2 ) iUAlPHZ  I 2 I *DI«E  I 1 1 i I i AL  PH  1 y i i I ,ALPh2y i2)  • 

• BE  T l.i  1 2 1 1 l>lu(2l,  02*i  I 2 1 • KDlyt2>*  wiJ/yt/li  wK.UUi/l*  nLLIly  i 2 ) 
COMMON  /Bhll/  TENSI  6Ai6A  ) I <>)ibA)i  (i2  i 6A  I • PM/uH  i 6A  I *'  Pl-A  jl.  2 t fiA  I i 

, CONMlNI  |0|  |M2(SI  iSVl  IbA  ) ,SV2  i 6A  ) iPARMi  S ) 

COMMON  /SINGLE/  F I PMAi  PNALE  • PM/ 1 MNOi  PMSi  (>NDi  PNx , NX*  NY, 
i HZ  » /.*  Z2*2ZE*  £NM  i If  INAL*  XZtRU*  y At  NO  i viltllH,  mLP«A»  ivN* 

• VuOT*  O.VUTi  HTi  TNEKGY,  ALPHAC*  Cb*  KtKwAC.  GAi-iMA*  Ei.j  CiAi 

• EJTKJi  RHT  * POUT  * DAi'EAi  ,2*  is*  tPuLsf*  As2  ipCn  *j  I • iCunl  * i Lui*2 1 

• 21*  RI63HX*  22*  RMXMA,  a3i  AP*  A. * cu>  Ht."N*  ONARI.m*  ilNaM-a* 

« , Ex*  PhIMX*  EPSMX,  ERPMX,  Dow**  Ml*  IUjIMAa*  V'tKM,  PHIKaa  I Mi  N'lft  1 

• PI*  IMAX*  UMAX*  NX*  MY*  NAD*  NX?a  NT2  I NAY*  NXDIM,  NyDIM,  Nl*l* 

« I PLOT  *N  I Tfc'R.NItuF  *NXM  *NYI . ,NH.  A(. « I;  *01  *02  #l>  1 »P2  *il*  ID  1 INK  III/, 

• MSKD12*XCtN«TLASTiNOHI  A * Iji'Um  *OCONu  mCON  *I)U  1 iIm  1U  mC  l20*Ii(,  a IN  t 

• HC22N  *HC2 12  *ALPH1  iALPh  / il  l.  I 1 *CUN  1 iCurtii  »tuu  *IHi.  *L  AU  *t.  ai*  * ..  i 1 * i i 
COMMON  /OUTS/  ‘lU'USCLFAC  *NRs*I.P. . i.'lC*  *.*tA  1 1 *NPuu  i »NPwNCM 

C*  • I I 1 *»**»«*  f * II  « 1 « .1  K It  « II  II  l‘  II  It  II  II  II  . II  I ■ t I 

COMPLEX  A (64  •(■<*) 

LUO  I CAL  U 
C 

C TAPER  r.OUMiARY  VALUES  A.ji 
C 

IE  (Lll  00  TO  12 
DO  1,!..  |.l, MX 
At  I *1 

A ( I *NY  I »('»■ 

A I I *2  I « A | 1 ,21  » ■,* 

10'  A t 1 ,INY«|  M »A  I I , (NY-  1 | I • ; V' 

DO  lid  J a 1 , N V 
At  1 ,J|  •••,(! 

A(NX*JI»  . 

A t 2 , J I .A  ( ,<* 

II  At  (NX-1  I ,JI  • A ( (!,X-1  | ,.;l  n»'. 

1?>  CONTI  Not 

r 

( rOVPi.  Tt  Itlf  INTI  'AIIY  At  i Ait*  « .:•  1 . * I|||M  Ml  |n(.„|i  ln|.  ' i,x  | I'll'' 

Tl'K’iVX  ■ r •• 

DO  '!  J*  1 ,NY 
DO  ')  1*1  , NX 

TENS  t 1 , j I » A | | , J I *C()N  Jt,  ( A I 1 , J I l 
IE  (TENMItJl  tit,  T l ’:M  X I (if  TO  U 
t MAX • 1 
J''AX  • J 

TENVAaTfHRI  I ,J) 

'I  ruNtlNur 

c 

C pi  Town  if  mot  pi.uTT  Img  H'OKIi.n  i.*a,:.i*ih 

c 

IF  t .NOT.  I I ) On  10  l‘j 
C 

r III  MUD  AH  MAY  H’l:  I'll'TT  l.-.l,  I . ; . * | i.A  ,i 

r 

no  lb  J ■ t 1 1, y p 

HO  111  1 « 1 , t : X ? 

Mf  L. , J"  T t ■»'.*,(  I ,.J) 

It  N!  ( I i J I < II  NM  I +NX2  • J*(ii / I 
1 1 NM  I MMX/  t.l  + NY?  I ■ I IOL  ■ 

III  f.iNTlMn 
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DO  16  1»1*NX2 
HOLD*  TENS  I NNX2#JI 
TENS! 1+NX2»J|"T£NSI I *J+NY2> 

TEHS1 1 * J*NY2 j «H0LD 
16  CONTINUE 
C 

C PLOT  FOURIER  TRANSFORMS 
C 

IF  (NPLOT.LT.SI  QO  TO  IS 

CALL  SYMBOL  ( 0.  <8.0 1 , 14  *3H2  ■tC.tJI 

CALL  NUMBER  I . it , 8. 0 • • 14  *2  *U**4mFB.9> 

CALL  SYMBOL  1 1 .32 *8.0 1 . 14 , 4H  KM  *U**4) 

CALL. LABELI0.U.4.0I 
CALL  PLOT  (2.S*0,t-3l 
call  symbol  is.o*  s.v*  o.it*  j.  o«  -it 

CALL  TOPOARAFI TENS*  NXOIM*  N Y D I M * NX*  NY*  0*0*  0.0*  lit  1U*0* 

• 10.0*  FPSO  t X2ER0*  HX * 4HFt.lt  1 HA  * -1  * » aEHO *H » *4HF6. 1 * lh i * 1 I 

DO  74  J«1*NV 
74  EPSOI 1 • JlaOtU 

IS  RETURN 


END 
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SUDROUTINt  VTHANS I A .NS  I 

THIS  SUt'HOU  T I Nt  OlANfilS  PHI  HACK  li,'  Pal  (OH  A I «flli  AUrtP  i o 
THE  COOMjINATK  SYSTEM  |0  i ME  (.  HANGING  H i AM  jIUi 


COMMON  / Hti  t»/  TENS(ftA*6fcl  . ollAAt.  (iJl&Alt  PIIAsL  1 1 6S  I i HIIAati!  l 6<t  I i 

• CONN  INI  10  I »M2  I 3 | t.svl  ISA  1 . SV  i I OA  I » P AHM  l t)  o I 
COMMON  /AAA/  bPb<6A.6Al.k  I*  JO  < bA  *6,  A I t MK> . , 1 1 • VY  ) i uu-  , • V • 1m  I < 

• A INI  ?«AAI  tl'ALPHl  I 2)  i HAL  Phil  P.  I .1,1  I , 1 , i)  . ALPltl  w , i I . ALP(|2w  , 2 I , 

. BEtli-12).  DIU(?I«  I>2"I2I.  RD1li2).  KD2wi2l»  aHiDlw.21.  aUNliJI 

COMMON  /SINdLS/  F,  PNA,  PNALI-.  PNK*  PNO.  PNat  P.ND  t PNk.  Ha.  hi. 

• h2  * i . /.itiiFt  aN.  i.  aHnAli  Afctnu.  i aiwi  1. 1 U i ei.  Alpha,  mx 

• Vulll.  OMNI.  HI.  rNl.KC.it  AUPIlAC  > Cm,  MUTUAL.  OA.mi-iA  , biO.  C.N> 

• t J TK  J . NH  I.  POOI*  0AM,  At  II  2 > i.l  . PuCub  < Aa/.PCK.jl.iC  OK  1 . ■ C.OK2  I 

• n > R 1 6 JMX  » A2  . R IMaH  A • a 3 . API  N.  lA.  IUMN.  DKAHfcA.  ,I.NoMa. 

• EX.  PHlMX.  f.PSMX  . ERRMA  . OGMXt  111.  MOIMAX.  VlCRMt  PH  1 MX  X dUN'/S  • 

« PI.  I MAX » JMAX.  NX.  MY.  HAD.  NX2 1 NY2  . NXYi  NXOl'1*  MYOlM.  NPT . 

< I Plot  .Ml  TH<  t NfUiF  . NX" .NVI  il.f'S  »l\Fl.  AC, .1) *1)1  *02  .Cl  ..<2  »*H  I U 1 . SR  I IV  • 

« RSHO  I 2 » XCtN  » T L AST  • M)P  f A *P‘|I  •.’  lOCONO  .OCON  *i*l,  I .OCX  I Cj  . i iC  A / O * 1 1C  A 1 N . 

• HC/  2N.HC2  12  t ALPH1 . AlPh/ *l>l.  T l .OON1  tCuN2  |H2*.  .H/N.t  Xy.I.AN.,.  I ImW 

ll  * ^ »***♦#  Pil,:  Cl 

COMPLEX  AlhA.AM.  AA 
1 1>I>1  *H2«KI> }<■  I NS. 

21)1)2*  il2»ttP2i'(Nf»l 
nlOINSInOl 

R D 1 0 ( N S I » 1 • > / 1 > l HNS) 

o;»n;>o in!) i* ( ( l *ai.phj  i mmh/mi?  i « « p»ei*l  p« , -.m.' k>  i 

fiNTI.'2»  MIRT  1 •)?  1 

N?p|MM«r>2 

ROiUMNMaI  . /n?  ■!  >| «.  I 

R:.R0I2»1.  /I  l.WT01«AKtn?  ( 

C 

C A I, I'll  1 . AL  Pll?  AVI.  Hill  AI<I  oli./vli.,,  I,  C'lHKI.Ct  | i,  I..,,.,  IMIMH. 

C r.Al.'V,.  IAN  l-r/HVAPI  !.  It,  Till  i'.H'.i.'-II.I  t>  1. 1 1 I-  ,<«L  f | o'l  •'  ‘I  ; 

c Tnri."Ai.  t<i  or” * 1 1 in 

ALPHI-AlPHIiIMiI*.  , b“  I I ..'♦A.'">„I.III,||||,V..  I «Al  I Hl.ilN.-i  I 

Al  Pit/ • Al  l>ll?i  I NS  I +■  I | ..  ♦(..,  "|)I  HI, 2 ' (0  I I * Ai  I N/.  1 | I 1.1/ 

■ I T | »IU  T I (M,l  I.SRt.  I / '.P  t|>]  I'M 

xfM»xrr*-  M”, ) .(•,  | H fi i 'i i ! 1 1 1 i * ■ i . • '<  i i hi  » i 

.‘.'■Til'  I'M '•":T"| 

41-I.VI  IVI.  !)A|.P|||.  I’Al.l'H/  III.,  I.  . I III.  Ml  L.'iO"  I I I |- ••■,..  I,,u 
( "'IF  IN  A(  I In"  A.'  I Ml  '•  ,.i  i I . . v In. 

■ » ' 

A X * .i  • i* 

\ • X ■ .»*!' 

U“  « > 

I ‘ X • • l' 

A Y ■ • • • 0 

r 

C a f > f ■» u v Oi  Af.r  ( MAV(,' 

I'n  •>  1*1  .NX 

I'll  | | •!  . I I.f  / IN'Mil  I I I *!'( /,.’  ■:•••  I l l-|.|  • . M.II-i.mMh  . 11.1)4 

. r,>«i  xn»i'  pm  1 . i ) 41 1 

A(I.|I"AII«1I  • C '-'  I ' L X ILHJ-  I HI.  I 1 I .Mil  I'll  I | I I 
S riiNTINUI 

DO  1 ' J*2t'.Y 

|.||  I 1 MI.'.IH'I  » (lie  / |f;*i.|  I 1 I +in  /.'fi*'  ,'  I Jl -i,  ■ ( Il-|>  I l-  ll.JI* 

. w r ,’  H xmh  i". 1 1 . .in  i 

A I 1 t.J  I *A  I 1 • J)  *C  I'l  X II  « II  III  1 I I 1 I I 

If)  (.'A!  I Nip 
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DO  14  J>liNYM 
Y-Y2&KO+IFLOATI J)-C.4I*HY 
DO  1 1 I * 1 iNXM 

PH|»u.4«Pl»(HC21N»Gl I 1 + 1 I +hCZ2N*(t2  I J+ 1 ) “liCONO* 

• <KT1«EX0*EPSQ(  1 + 1 « J + l ) +V<  I 2+LXN+LHbl  1 + 1 » J + l I I) 

At 1 + 1 «J+1 1»AI I + 1 • J+ 1 ) +CMPI.X ( COB  I PH  I I ib 1M PHI  I | 
X*XZEH0+ (FLOAT ( II-liSI *MX 
AA«0.4MA(  I + 1 ♦ J 1 +A  ( I » J ) ) 

AA?»CONJlitAAl  »AA 
SU,MX»SlJMX+AA2 
AX«AX+AA2+X*X 
OX.OX+AA2+X 

F AC  TOR*  A lt'ACi  ( CON  JD  I A A 1 + ( A I 1+liJl-nt  [ i J ) I 1 
C X + CX+FACTOS 
OX»DX+FACTOH«X 
A A»0 i 4* ( A ( | iJ+l|+A( I i J | | 

A A 2 »CQN  Jli  I AA  1 » AA 
AY«AY+AA2*Y»Y 

DY«l)Y  + A|N‘A<>tC<>NJ(il  AA)»(  A(  1 ij  + l (-At  | iJI  1)«Y 
14  CONTINUE 

.4UMX  • S'l’-'X •DAMt A 
AX«AX  + |)AMrA 

HX»PX(in/\prA 
CA-CX  + IIY 
l)X»Dx«HY 
AY»AY#l)Ahl  A 
l.'Y»DY«IIX 

KDEIiOM*  1 • 7 ( Ax# ‘lil.  iX«llX«l!X  I 
DAUPlil  I M,|  «,  , !>*  ( I'X  + MI.  X-UX-'CA  1 f .-1  -t.i h..  J 

dal  PI  'if  <Mi)"  • b K 1 1 Y / A Y 

cult  l (MAt»(  Ax*cx-ux«i'X)  tipd-iM.:-; 

FT<V<h  VAU'I-'R  Hi!  U'Tc.K  Ut.TPilT 

x*ii>lot 

AOUTI 1 iK W 
ACPI  I 2 «K.  | ■!) 

ACWT  I •lit:  I ■1.1 

Aljut  I <i  ,K  1 «■"  2 
AOUTI  iiX  I ■ALl-Hl 
AOUT  I (i  iX  I •ALPH2 
AOUTI 7 |X | *KI  T 1 
ADI. T I U iK  )»|'AU'H1  (MI,  | 

AOUT  I Mr  I ■llALPH?  (N>,  I 
AOUT  I 10  iX  I ainiLt  1 I Ns  I 
AOUT  I 1 1 IK  I ■ X C t. M 

if  i had  i>;i.  m im  ru  ?>: 

DALPul  ( + £>>■  • 

0AUPH2  (Mb)  ■ 

.HU.  T 1 1 Fix)*'  . 

2 0 CONTI  NUT 

c 

C I iMI.'A  1 1 TmI  UllitA.v  »*;ii  Vfl 1 *■  i .\PT  I 0 rii..:.i  lii  ... 

r 

ALPlilOl  MS.)  •ALMM1  + I>AU‘M  IN!.  I 
ALI*li?>.(F!M  ■ALP)i?+I>AuI,I!?  IN). ) 
tlTTl,  INM.Ol  T 1 * 01!  f T 1 INF  I 
Ml  K'lHSI  >XCt  I 
C 

C r.O'PUTI  iJUAAT  I T ll-.S  In  Tnl!  IM  111  M i.'.U  !■  >"l  I 

C AliAIMH)  Ct’OKI'INATl  S. 

C 

IF  INH.AU  ilII,  II  1,0  TO  ■)>■ 

Pi  *lil) 

P 2 > , . 0 
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IT  (NMS  .10,  I 02*  '.u 
M/l*COI‘ll*H* 

►f*2«CON?*HZ 

2DU*m/i«hi:Jim?I 

2D12*»V2««l.|..(2! 

202  1 *h/  l • 1/ 1 > 2 < ■ I ?.  I 

20?2»H/2»'»I)2i  (2) 

1)1 1*1)1  u I 21*1  (1.  ,+2..J»AL^HU.  (21  *20111**  2 + /0  11*2  I'll) 

l)  12*0  l.it2l*«  1 1«"*2.«.  »AL^t'U>(  'l  t « /!'  1 2 1 » *2 +2u 12 *201  2) 

021*0  2. '(?l*(  ( t#.i*2«..*ALP»«o(  2)  *2021  I »»2 ♦2.D2  l«ZU?  1 1 

1)22*02" ( 21*1  1 1 *'i *-2*. ■«ALIJH2i  121*2  0221  **2*/.iW*/022  I 

111)1  1 ■ 1 .()/nl  1 

ill)  12*  1 *0/012 

111)2 1 * 1 » 0202  1 

N .)  /’  2 • 1 • i."  2 1.  ? ? 

lliiWl)!i  1 » r.Oil  t ( IU>1 1*102  1 1 

llSKl)!.2»!Ui|(T  I ilolJ}*ui>?2l 

h 1 1*M2*  ( MC l 12-H2  l»lir>H('f.l-H2?«lir<KOfi21 

WT2*HC21Z«WT1 

NHAU*1 

flKTOHN 

K NO 
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5UHR0UT  INI  OUTPUT  ( l.P  I 

THIS  SUltROIlTlNt  STORES  OUTPUT  l)  A T A OliKIMti  THi;  CAUCUL  AT  lut»S 
ANO  PRINTS,  PLOTS  AMO  PUNC'ht-s  VARIOUS  OUTPUT  DATA  '.ilil.li  Tilt 
f ALCULAT I ONS  HA Vh  CONCLUDt H . 

u *»#*««»• 

COMMON  ZAAA  / tPS<  hA,6SI  * fcPSUt 64 *NS I t AOU  t I ] 1 * S»‘>  > .POUT  I V t 1 0 I i 
. A INI  2 *64  I .OALPHlt  2)  *t)ALPH2  I 2 I *1)10  T 1 ( 2 I •AlPhlOU  I ,ALFH2CI2  ) , 

. Ht  T 1" ( 2 1 t l-lu(2l*  |.2i.(2l  » NO  lu I 2 ) » i<u2;(2)t  ji<fl>U(2l*  XLI.Nj(  2 I 
COMMON  /*»  9V/  A 1 I 64  *04  ) * A,!  I NS  * A*.  I . I T l.NSI  N4  *64  ) 

COMMON  /I'.Hh/  Tt.NSI  hS  *64  I • Lill6*i)t  vt2l641*  PnA.'il  1 1 faA  ) * Pi  ii\6l.2  ( tiA  ) * 

• CONMlNI  b)  t M2  I 3 I tSVl  164  1 t SV2 I 64 ) *PAR-H  b.) 

COMMON  /StNCiLS/  Ft  PNA  t PI, ALP  • P IK  * PNC),  PNSt  PNU  * PNZ. » IX*  MY* 

. HZ*  2*  Z.ZtX.ZK*  /MM,  ZFIMAL*  XZ'.ROt  Y2l.R0*  Wlt»Th.  .'LPil,,*  ,<N» 

• v Oft  omdt*  hi*  r i-ii- .R(* v * alppiac*  c.i;*  im>.<ac*  oav -ia*  i.tj*  its* 

• LJTNJ*  RUT  t POUT*  DAI,  I A * -2*  T6*  TPliLSL  t A..2  *PCR  *Sl  *TCOR  I * I CO  (2  * 

. Zl*  RIAAMX*  22*  RIMX'.'Xi  23*  Apv.N*  Z4i  H/I-Nt  Us  A HI./.  • lINMUt 

• LX*  PHlOX*  IPSMX*  LRRMX*  (XiMX*  kl*  HO  I MAX  * VUHM,  Pll  I MX  X * llZNvS  t 

• P|*  IMAX*  JMAX  * NX*  NY*  MAO*  NX2  * NY  2 • NXY,  HXjII'i*  f I Y 0 1 M * NPT  » 

• I PLOT  *M  T|  ,1  *.M’,U*  *.',/..  i *I,Y.'  */,  * ',P|  Au  tii  * } 1 * l >2  * u 1 ti  2 *:•  k T z 1 * SR  1 1)2  * 

• RSRu  1 2 1 XCI.II  * T L AST  * SOR  I B * Pol'  1 , uCwNi'  *,'CU-N  • i'D  1 *m(.Z.  1 0*  / iC2  20  *nC2  IN  * 

« i iCZ 2N  *H(  Z 1 2 *«L  PUT  * i\LPi'2  *■  ■;  T 1 • C.w.'-i  1 *L  ,/>N2  *,i2»  * , u. N * ,. a,.  * i.x, , * 1 1 *,.  T 2 

CU  IMUR  ZOI.TSZ  ,il  1 1 SClP  aC  *Ii/*!  *10  •>  * ,(.  ••  * '.L  A I I • , |i  l. I *, ; P.J,  il*  • 

CO.  IMOM/PUN/HIO 

* u ti  )>  )i  n , i,  „ i,  , ,i  „ i.  ■:  :.  ..  ...  > ..  „ v,  i.  ,,  i.  , i,  i. 

C uw(»l  h x Alt  A2 
LOli  I ( A | Li' 

OlMtNMO.'!  /il  I I'll  * ..PRlIifll,  I*  I ) l,  I * I t ( |:  I * I L i , I ) .' ) 1 1 i-L  u I ( /,  i.  I * 

1 Hit  Ul)  *LlM.I.AH.  ( 21  tSUl.AI  li.  I 

0 I Ml.  *11*  I ON  A V I ( H.  ■ > t AA1M8  I.  Xl-'IUC.)*  XMxIHtl*  X-  Y « H(  > • 

• HZ, fill'  I*  t i", "AX  IS  It  /C,K|U  I*  PH  | *iAX  ( H,  | 

1 1 A T A ( x INCHI  sxit  I * ( y | ; in  I 

f i A t A ( I.  I Ml  |.AI'L»?  I In  II 

(. 

1 Pu.'t  •*  I PLOT*  1 
|PlOT*M|M.,  ( l PL  (IT  . rt  | 

P7  RAX,S»PPllT*f  X 

I T ii,\Nf-»rNF  R(iY»l  X 
S(  P AC  1 ■ V I 0 T Htt  SRT I ' I 
SC  FAC/«Wl  hTM*SI(TI'? 

SC.AKL  A.  sr  I AC  ] »SC  FAc  2 
R I PP  Ak«  1 ZTI  HI, ‘-IX 

c. 

C tiif  aki  /•  » u.x  fw-M  / 1 ;«i  .>  •,  in  i.j  !/*;.:■  (.«»•.  r i-u«\  i.:  vi  i. 

M IN'  i - 

HO  4 6 1-1*1/. 

AN  MIMA  I I I r ARI  A II  I ■ 

HO  s?  J« 1 *NY 
HO  *>,!  I*1»NX 
sin«',ii'i+ii  Nsi  i ,.ii 
F RAC  1 1 C NS  I I *J I HR | P|  AK. 

I'O  M M.«l.lu 

II  I HI  A c ,I.T*  ( <P:M  I h ( II  Ul.  T I : 

DO  AH  L « S'  * 1 1 ■ 

SUM  A ( 1.  I ultlll'A  I L I + TLN.M  l * J I 
A ,i  A HI.  A I L I ■ Alt  t.  A I L I + i jAPl  h 

HO  TO  S? 

>/.  CUNT  I NUI 
CHNTiNUt 
C 

C rOMKrCT  UNITS  ANu  NIIk"AU/l  API/'.  ,M;:  l|UXI' 


.Sul'*  '.i.'HhAI'l  A 
00  '-A  1*1.1/' 
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APRINTI  I ) «ARt  A I 1 )#'SCAKEA 
AF ( ] I “AREA I I ) /AREA ( 9 ) 

FF  ( I ) “SUMAI  1 )*l)AREA 
FLUX!  I ) “FF  ( I )*PTRA|,S 
R i ( I ) -FLUX! 1 l/APR INTI  I I 
SA  CONTINUE 

INTERPOLATE  TO  FIND  Tut  AREA  CONT  A I N I NO  0*63  OF  TOTAL  POWER 
DO  661  1*1*10 

I F ( FFI I) .LF.0.63)  GO  TO  661 
1 LO-  1-1 
I H I * I 
GO  TO  662 
661  CONTINUE 

IF  0,63  EXCEEDS  RANG-  01  FI'*  US*  LAST  TWO  FF  VALUES  To  L X T KAFOL  A 1 E 
ILO.'J 
I H I » 10 

6 6?  AP«APR I NT ( I LO > ♦ I ( ' I *63~F E ( I LO  I I / ( I I t I H I I -FF ( I LO I) I * 

1 I APRINTI  [HI  ) -APR  I NT  ( U.Ul  I 

C 

C FINE  AREA  RELATIVE  TO  INFINITE  ('All.'  S I AM  |)E  AM 
C 

A P 1 * P I » W 2 * I > 

A RE  E * AP /API 

R I b3 ■ P * 6 3*F‘ TRAMS/AI* 

H I ,M X » T F NSHX*PTRANS/ SCARE  A 

C 

C COMPUTE  OTHER  OUANT  I T lf.s 
C 

/ C OR  I I PLOT  ) * 2 

AVI  ( IPl.Of ) ».  *6  3*  PT  SANS/ AP 

A6TI 1PL0TI-AP 

XM  I I I PLOT  | »IM’IX 

XMX  I I PLO  r ) ■ I I I'AX-M X2  MSCFnl  HhX 

XMYI I PLOT  I ■ ( JMAX-NY2l*StFACZ*HV 

Nil  2 I 1 F’L 0 T I “HZ 

El’SMAX  ( IPLOT  I »t.PSHX 

PHIKAXC IPLOT ) “PH | MX 

PARHt  IPLOT  1*1. o/VTI  MM 

II  (AVI  I IPLOT)  < (it  a KI63I1X)  2. 1*2 

M I 6TMX *AMA X 1 t AVI ( IPLOT ) * K 1 6 3 MX ) 

IF  I XM  1 ( IPLOT)  .GT.  W|“X»X)  i'2“2 
R|MXMX»AMAX1  I XI ' I ( I PLOT  I i HP'X-'X) 

IF  (A63I  IPLOT)  .IT.  AP'.'N  ) /./»/ 

APMN*AM|M1  ( A63  ( I PLOT  ) * API'*'  I 
IE  1/  .01.  /FINAL  I GO  TO  2*1 
IF  III/  ,LT.  HZ '-'.M  2*  */ 

II  / MN“  A"  I N 1 I HZ  t H/MM) 

r 

2*1  IF  I • NOT • L P ) 00  TO  11 
C 

c 

C PRINT  RE  MIL  T S 
C 

IF  INI  XlT.NI  *•')  PRINT  0 V 
1 nd  FORMAT  ( // »#  PMI  ' A T'.JKI  I Xl  I ''•»//) 

KLl  I “K  1 MX  * AP  1 / PTRANTi 
Tl  NSMX “ TfNSMX  “HAKE  A 

PRINT  1 (a '!  * 

I AM  FORM  A II///  AfaXl,ill«“l  KIM'IlN  •:»»//) 

PMlNI  l'i  * / f ‘.'I'  *N  11  I K ,1'EluI  *P  I MU, MAP  *AI.  L «!•  111.)  *■•  I X • E - I I 
i'H  FORMAT  I ? A X 2 HH  I III  >.  A|  ( UL  A 1 I Or, i»l  All. I l>  =1  J (','•*  i X ..  t - 1 XM  )/ / 

X 1 / X 3 6 1 > I II  l SHI  OVIR  NOIU  AE  1 /I  l:  lUHT.lIt  •'  l!.!i// 
x /'1x;mhthE  mimiek  in  /-Ml  pn  “I**// 
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X 12XA0HAVERAf>E  POWER  ( KW  I EMITTED  AT  APERTURE  -F10.3.// 

X 20X32HAVEHA&E  TRANSMITTED  POWEH  (KW)  -F1P.3,// 

X 1 AX  3EHAREA  ( S(JCM I CUNTA IN  I NO  0.63  OF  POWER  -F10.3.// 

X 17X35HA  RED  (RELATIVE  TO  INF.  0AUSS1AN)  »F10.3.// 

X 10XA2HAVERA()t  INTENSITY  (Ku/SUOl  IN  Tn|S  AREA  «F10.3l// 

X 26X26HPFAK  INTENSITY  (KW/SOCM)  kFIO.3,// 

X 12XAUHI  HtL  (RELATIVE  TO  INF.  UAUSSIAN  I'lAKI  ■FIO.S) 

C 

C PRINT  NUMERICAL  DATA 
C 

PRINT  AID. 

A 30  FORMAT  ( 1H1  Si-X22H*»*  R 1 C A L.  DATA  **«//) 

PRINT  160. 

1 60  FORMAT!/ 

. 9X  FRANtiE  ( KM  I * (SX  **AS2  ( Ck*2  ) » 6X  <»AS?/A20«  6X  »TP(SECI«  7X 

. #PCH ( W I * AX  * T SAT  l SE C ) * 2 X » 1 SAT  ( ",/CM2  I * SX  •TUKbCOHl*  SX 

. *TUK«fOR2»l 
DO  17  ..  I ■ 1 1 1 

170  PRINT  100.  ( HOUTIKi I I *K«1 ,41 
1 HO  FORMAT ISX .3'  1 3 . S * 6E  1 II . A ) 

PRINT  A20. 

A 2 1 ' FORMAT  ( /lXAi  IINOR'iAL  UH;  A'lPLlU.llf  SAI  P|.l  S AT  APi.RTURL/  I 
I P X » N X / 2 
IPY-NY/2 

PRINT  ADO.  IPX ,NY  » t A I N ( 1 * K ) » K * 1 » NY  > 

ADD  FORMAT  ( SXSMAT  X » I 3 » /X AMY • I 2X2HTO  I A * I I t. ( 32F  A • 1 • / )) 

PRINT  AID.  IPY.NXi ( A|NI 2 .X ) .K»l.hX I 
Al'.  FORMAT  ( SX5HAT  Y»I3,?XAhX*  12X2MTO I A » //2 ( 32FA  * 1 » / ) ) 

C 

PRINT  1?.' 

121  FORMAT  ( //SiXlHZ.  /X2MII/2X  1H(>6X  £111  > 1 1.X  ;’HU2SX0HAlPUA  U Xtil  IALPUA2 

1 3X  fMIlE  T A 1 2X6Hi.'AlPH1  2 XPHiiAUPn;.’  IXtMi.l't  UAXAMXCt  MAXSMI.PS.'iX 

2 3 XftHPM IMX AXAHPAHM/ ) t 
JPU)T*lPLOr-l 

DO  1 3v/  K.l.JPLOT 

130  PRINT  1A0.  AOUTt 1 »K) >K22 IK) * ( AOUT ( I ,K I , I «?.U I * 

1 EPSMAX I K I .PHlMAX(K) • PARN ( K ) 

1 AO  FORMAT  ( l XF  A • 3 • I fi  . I * 3F  7 . A * M H . A • I.F  . 1 1 F i« . A *Et)  . 1 ) 

C 

C PRINT  OUTPUT  DATA 
C 

PRINT  2?n, 

22l  FORMAT  ( lHl  SlX1‘'M»*»  UUTPUl  U/vTA  .k*//) 

PRINT  62 

62  FORMAT  ( / * 37X  » ARK  A*HX»I  I UX  *6X  « AHl-  A « AX » FLUX  « ti  X . I HR  AD  I a\<  L* 

• /t/6X  «LEVEL*  AX 

1 »C>w  CM  I «7X«  (KWI*AX«FRAC  I I>ii(»2X«H</'t  T |i'N»3X«  Il.v:/S..«C  ')") 

DO  66  1-1.1  ' 

PR  1 NT  6 A t CON-' PM  I I .APR  I 111  l I I .1  LUX  ( l ) »AI  I I I .1  ' ( I I « h 1 (I  ! 

6 A FORMAT  I 2 SX.F  6.  A,  r.  12.3,1  l / • I »F  / . A ,1  1 (’ . A , I 1 <,M 

66  CONTlNur 

PRINT  7 , R I h'J.MX  . /.  I . R I X ■-  X . /.  2 . AP  '-1 . /. .'  »M/.--  -I  * 2 A 

70  FORMAT (/AUX*  MAXIMUM  AVO  I ■# 

7 E 10. 3 • 2X  **AT  2-A.ElU.3t/ i AOX . * MAXIMUM  PtAK  I ■ ».fc 10. 3 »2X «*AT  /•*, 
a E10.3,/ .AUX .*  MINIMUM  AREA  t E 10 . 3 . 2 X .» A T l « * . LI  0. 3 . / . AUX  . 

9 * MINIMUM  HZ  «*,Elu,3,2X.»AT  Z»«.Rlu.3) 

PRINT  190 

190  FORMAT  I / /31 X 1HZ9XAH I AVt  7X 3HA636XAH I MAX6X SHXPl. AK 3X5UYPI.  AK  1 
DO  200  K-l.IPLOt 

200  PRINT  21D»  2COK ( K ) . AV I ( K ) . A63 ( K ) . XM I ( K ) . XMX ( K I . XM Y ( K ) 

210  FORMAT (29X6F10.3I 
C 

C PUNCHED  OUTPUT 

C IF  INPUNCH, EO. Cl  00  TO  230 

PWlO-WlDTH 
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1QCW-2HMP 

IFlNCW.EO.il  I 0CW*2hCW  ii*.  v , _ . , 

PALFS-0.  [(  Pf!.A  8 /'  ' " ' !' 

PVO»VOOT/HT*»Ol  <?  . , , 

PALF«ALPHA#l.Ei 

PE»ENtRl>Y/lU00. 

OM-OMOT/HT 

PF»F/1.E5 

PUNCH  169  »N  1 1),  lOCW  »N8M  ,NAl)  »PN  A .PNALF  »PhX  .I'NO.P'IS  *PND  iHFii  i 
1 PF  .WN.OM . PE  *HT  »PVU  i CALF  .HALF  5 

169  FORMAT  I A6.1X.A2.8H  NSHAP  t » . 1 1 .6X  . AHNAD-  , I 1 , Ft)  .6  ,F  7 .<>.  2F  / . 

1 78.3»F7»<t»/,A6,Fb.2»P'7.3»F'IUj3il'8«4»F10>3,P9,{i» 

2 ,F7.4I 

PUNCH  73,  NllJ.HX,HY,H2»l'hlMXX.MX.NY»NAC'.lHL0T 
73  F0KMAT1A6.4I;  h:.3»4l!)l 

PUNCH  7?  . M 1 D .PNOtl  tCiCUMO  .OCliN  »h,>  1 I’M  AC  • UAH.i'A  , PUU  f 

PUNCH  72  , NtU«21.Rl(i3MX,2P,«lMXMX,?.  l.APMN.ZA  .H/.’IN 
PUNCH  72,  NID.TCORl  .TC0R2.IU  *P2,D»XCl  N 
7?  F0RMAT1A6.8F. 9.21 
DO  71  K»  1 » I PLOT 

71  PUNCH  72  » N 1 0 ♦ 7.C0W  1 K. ) tAVlUl  ,A63U  I .Xf.ltKI  .Xi-.XIK,)  tX-IYUl 


: 7 , 3 , 
i ii  » F 7 . A 


PLOT  FLUX  FRACTION  AND  AkL A FRACTION  V£.  I KRAO 1 ANCl. 


IF  (NPLOT.LT. 31  GO  TO  7 A 
HO  76  |»l«y 
FPLUT1  ll.AFI  1 I 
FPl.OTl  1 + 9 I ■ F‘  F 1 1 I 

coNTinur 

X t f 4 c : H E fi « R • • > 

Y I NChtb* 8 • . 

CALL  PLOT  1 1 » * • t>  *“3  I 
CALL  LAPTLC  .1.6, 2. n| 

CALL  iikil  1 LUliXIM  , I-  PLOT  ,9,2  , 1 , 1 « V »X  1 , iCl  IL.'i , Y INC  HI  I,  , 111  , 

X 1 31  iP  LUX  F WALT  ION ,10, 

1 L 1 Ntl. API ,“  1 , 1 , XI . 1 U , )X  , Y.’  l'.,i,Y,YT  ICK1 

IIY  Ah  !U  TUhNt  P uY  CAI  L t,i  Mil  .-..nl*  1 M-.  ' I ,Hia  I :>  it.  w 


UNIT  OF  Y,  AMP  YTIU.  lb  Pl.'.UMi  1 1:  ll'. II  !.  ii.T'.I.iM  11(1 
PYIRAC,  f 1(.  . I :■!  LU".  I i»  I Vi  I IsC  nl  . U I ••  ; i , ..Xlb  , I..';,,  t A I K i; 
DYi  KAL-YT  U.K/liY 

PYF  Ll  IX • IT  YF l<  AC. x h PI  l-.'liY  M.l.T  “I  .1 II  J 
PY  A III  A«|.YlT!AC:»APP  I r.r  I '■  I 
Tfc.-'PX  ■ 3 » P. 

(ALL  SY.* '.POL  I T I ■•■l,X  , - . 69  , * ■ , 1 1 6,2  1 1 1 PPM ' 1 AMP  1. 1 r.  .•  / 6 t i •),«■, 2. 
r 1 1'ipx  > lu<iF  ( i<  i I y > / xc.IM  i »ux 
T r.'TI’Y *rPLOT  (9l»|)Y 

CALL  SY'IHOLITL  PX,Tr:.-.PY-.2ii.:,.A*>‘.',7lji  ni..  i . i\  . ■ H.L.ll 
CALL  !, Yl  itJOL  1 TEMPX-.?.  , T t-'.PY- . i 6 , . 1 , 1HA  ,;  ,1) 

Ti  MPy«FPi.PT  I 18)  #PY 

CALL  SYMMOL  1 TITil'X  , T|  ;..l'Y*«2l|-  . . A bb  , 7'.*;  L l . . . . IL  .)  n , 1 Ml . i 

CAH.  SYMHl'L  I tL  IPX.  .2v- , 1 1 M-Y-.i  6,.l  ('!’  , lliF  ,'i  , I I 

CALL  AX  I S 1 -1  • , ••  (III  IF  LUX  ( l.v.-l  ,b  , Y l - ,t.i  II  : ■ * 9 1 • * , Y I U ,. , ' , ,.Yi  LL> 

CALL  AX  I A lx  I MCI  U *.  • . . , 13HA,<I  A I KAi  7 | i.h  ,-13, 

1 Y|l.flil  •,  ,91).  ,Y|  Iv..;  , , .ilYlKACi 

CALL  AX  I S I X I HC  HI  r.+  1 • »■  • , 1 ImAWi  A(  i,  ( I | ,-  |1  ,Y  I *IU  i.  '<0.  >>■  I 
1 PVARl  A.I.nl  H .2  I 
CONT  l l.l  1 1 


PLI'T  AVLWAGF 


Pi  ak  ir-Ti  k;.  lilt;. 


.)  I LiPL  “ 1 1 »i.‘ 

IF'  INPLiJTil  T .2 1 III'  TO  H. 
DO  7 6 1-1,1  PLOT 

i plot  i i i - a v 1 1 1 1 
f plot  i i + nu.ni  i -xc  1 1 1 1 

COST  INUI 
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call  PLOT  t D I SL»L  » 0.5*  -3) 

CALL  UNI  PLdTUCUK.  FPLOT  . IPL"T  * 2*1.  V*  iPLlll.  XlNCHtftl  YlNCIU 
. * In  • 1 »«H|  (KW/CM2  ) >V.L  INLLAnl  i-l  » 1 *XM  I N *I)X  . YU  I N n)Y  • VT  KM 

TtMPX«o.5*X  INCHES-C.O 
CALL  LARrLHJ.t'b.B.n) 

CALL  iiYMUOLl  TLoPX*  — v..Kb«  bHZ ( KM  I . O.bl 
TCMPX*  I ZCONM  I PLOT  ) -XM I N I »PX 
TLMPY«  1 FPLOT  < I PLOT  )-YIJNI*l)Y 

CALL  SYMIAOU  TFt'PX-.  .A  i T FIT  Y-- . 3 « ).10b.  A'UAVti  b*  A I 

TbMPYo  ( FIHC'TI  ?.»  1 PLOT  I -YMN  I *UY 

CALL  SYA'IUlUTEOPX-i  .A»  U.VPY+P.3,  0.10*).  AHlPAX*  L't  Al 
CALL  PLOT  ( 2 » b * 0»..»  -'4  1 
CALL  ,'iYVHOL  K . « H . f , < 1 A i 3' 1/ 

CALI.  MUMP  I'  ll  I • 36  * 0 • * • 1A \t.  * t>.  t Alih  H • b ) 

CALL  !.Y*riOL  I 1 « 3?  |H  »L  I » I A .All  IT  itl.tAI 
C 

C PLOT  r J NA|.  CONTOIIftS 

V 

HU  IF  IMPLOT.LT. II  (»0  TO  II 
CALL  PLOTth.-  . .4  I 

CALL  LApPLIt  .ll.A.I  I 
call  plot  (2.b*>.*-3i 
XCKVTLW-S. 

CALL  SY'IHOl,  ( XCPN  TIT  » b » * • 1 A * 3 . ' i — 1 I 
XX/I'IH)»X2I  HOP’iCFAl  1 + XC.l  M 
YYZl.  HO«Y/t.Nl'»FCPAC? 

IF  INIK.NL-  ■ ) fin  To  bh 
X P'OTt 

Y INC  Htr»»  1'  •• 

|iX«MXlL'>CFACl 
PY  AlY'iC.I/itV 

I X»-A 

If-  I vp;. • fi . . I Pi1  to  l i 

aljiot  i | :a{  (.!•:, roui,  plot  In  oil  u;  iru-v  ax  1 1 •• 
contii  nr 

,'  CLFAl  • | ”, 

X I'  AC  ■*  L , 

! >■  ■■  ■ J»l  i '/ 

i-r  * a 

II  I f L L'l  I I |JI  .1.1  i.u  n • i.  •...  . 

YY<  . iO-i.  YY,  ■’  l.l  -II  I.  0.1  ( .11  ■ I.  I i/i  1.  ii  I I ' Y ) 

n.  I - ' 'I 
i II  T 
. I A|  i --'.(  I I A C 
X .(■  A(  I = YY.'I  ii./  A A|  ' 

Y'  C A|  i -iYV.'I  HO/:,(  A | i 
Ii  I Y i .IT  . 1 . ) (.  i T 1 1 |- 

II  I 1. 1 AC  .(-1  |I.‘>(  I XI  ■ C * / . 

M 1.  F <■«»  « I.f  L • in  “ Xi  A( 

no  In  '/'/ 

y II. ( iu  'hi  .■  «x'.c  Ai.l 

Y INC  III..,.  1 . Y'»C  M.  I 
XXf.l  P'li-i  ( I I AC 
YY/I-  !<■•■-■.(  l.l  A C 
I X " -i  1 
f ON  T I NO! 

CALI  T-ipom  a.  i ( li  i".  .-ixi'  |;  >i  I .'  1 1 1 . ■ I 

1 I I t 1 1 1 . v > . 1 i-1: 1 * I X . A l II  (, . 1 i A 1 1 X ( l h O > .’I  l.i-f  y ! I :..  * I . 'i  t i v . I 'I  I 
|.  I Tl  VA. 
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5UBH0UUNE  LAHhUX.Yl 

COMMON  /MNGLS/  F*  PNA*  MNAlF*  H':k*  Hf.O*  HN!>*  PNU*  MNZ*  MX*  MV* 

HZ*  it  ZZ./-ZF.  ZNM*  ZFINAL*  XZhKO.  Y2l.NO*  /.IDTm*  ALPHA*  aiN* 

Vl/DT  * 0*11/ T « HT  * LNKRCiY.  ALPHAS:.  CS  • I<1  M<AC  * GAMMA  t ITJ*  CTX. 

KJTKJ*  RMT  • MOOT*  OAKtA,  *2  t T.S,  TPuLoL*  AiZ  *MCH  . S i . TCON 1 • KOR2  * 

Zl»  HI63MX*  22.  NIMXM.t*  2 3 * AI'M.N  * lUt  HZI'.N*  DKANtA,  TtNSI'X* 

EX*  PH ! MX  * EPSMX  * ENRMX*  DGMXi  lit*  HOIMAX*  VTLUM*  PH  I MXX  .MZN'/a  t 1A 

Ml*  IMAX*  JMAX  • MX*  NY*  MAD.  NX  2 * MYZ  * NX  Y t NXUlM*  NYDIU*  NMT  * 

I PLOT  *Nl  TEN  .NHMF  *NX,"  *NYM*,V1S  ,MF  L At,  *|>  ,0 1 .02  *M  1 * M2  ♦ tRTl!  1 . SKTD2  * 
WSKU12*XCCN*TLAST  •!>0HT8*HMIX' t&CUMP  *OCuN*HOI  *MCZl0*MtZ20*m.ZlN* 

HCZ2N  * HC  Z 1 2 * AL  PH  1 * AI.PM2  * l>L  T 1 « CUN  l*  COM2  *M2U  *M2li  *L  XL-.tXN  *1/ U *i'.T2 
COMMON  /OUTS/  NHvi,SCLFAC*NH*.*NM'I,iic*'*NEXIT*NMC,01  * NMUi'iCM 
X 1 *XYO. 7J 

IF  ( NUM  . L T , 0 1 CiO  TO  1- 

CAUL  SYMHDl  IXlY*i',l  , SIINALF  « it:  ♦!'**/ ) 

CALI.  NU'lnCH  (XI  * Y * • 1 *M,  lALF  *u*  ••  * Al.l'  7 • 3 I 

Y-Y-.2 

CALL.  SYV'IIOL  ( X * Y * 1 1*1)11  NX*  *•'.'■*':  I 
CALL  NIIMIU  Tit  XJ  *V,'.i*  t *M  > i'*i  ,AMI'7.  J) 

Y*Y-,2 

CALL  r.V'  HOL  I X ♦ V * ' • 1 *bM  NO*  *u . ..  * S I > 

CALL  NO’ '!'l  HI  X 1 * Y * i' ♦ t *Mi«r *;/ *i  t A' ii  7*  1 1 *" 

Y*Y-.2  f , 

CALL  ;.Y.*'IU>1  IX*Y,  • 1 *1/11  N,-i»  *•■*"  *'■  I 

CALL  Nl/I  "U.  M I X I * Y , . • l *l’N'.  * • • *Aiil/.-l 

Y«Y-,2 

CALL  ‘.y  MM.  ( X*v,  •*  1 ,MI  'll  * *r.  t " I 
C A|.l.  Ml'im  MIX  1 *Y,'  , 1 *M  •■'*•  i,  7.3  I 

I/I  T|i!/‘, 

( '"IT  If, Ml 

n » ! / 1 ■ ■ . 

< A l 1.  .villi  ( X * Y,  , 1 , /m  I •*■,.,/  I 

( A|  l .,|,IV!  • (XI  |V|  • | <1  - i * * A i ‘ I I 

V .i  V - , . ■ 

f * I I.  • I l I ( X t > , * 1 * / 1 I 1 ' • * * • . * / I 

call  m.-"  -i  .i  I x 1 iY  • < 1 1 l ; 1 1 • • ■ */.i  i /«,■  i 

V«Y-,,' 

AlI  -■  A l »•!  IA  • 1 • 1 ' 

CALL  :.Y'i"l'L  I X * Y*  >.  J » 7 1 1 .LI’m/i**  • *fl 
CALL  f.Ul'ii  Id  XI  *Y*  *l«/.|.l  i * * A i 'I  '..'I 


Y * V - . 

CALI.  '.Yi’m'i  ( X * Y * * ] > / 1 1 ■,-•*  • i .'  I 

CALL  'll  i‘V  c M 1 X 1 * Y , .1*1*  ■ I'M"  '*11 


Y » Y - . 

V *V  I I /lll/l  • 

< .‘.LI.  ’ Yl’l  01.  I X l Y i • 1 * /I  i . « > • l / I 

(ALI  fill""  I < I X I i Y * . I i , ■ I'M  li  /.  II 


V » Y - , 

l,'«l'l|lT/HT 

CALI.  1 Y.'T  l l.  IV*  fi  i | * Ml  l •:  l i ‘ * • • * 7 ! 
C.  AL  I Nil,  "M  I X 1 • Y * •*  1 *1  i.  > ■ I 
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c 


c 


CALL^SYMBOL  <X«Y*u.l *7H  0T»*3.0*7 I 

CALL  NUMpEK  ( X 1 • Y »()  • 1 »HT  n.  •(.'•'♦HI'  7 »5l 


call'  SYMBOL'  I X » Y*''«  1 * 7HF  VtRl'Y*  t C • n »7  I 
CALL  NUMBKH  ( X 1 1 Y *''«  1 *t  Nt  HOY  »0  • C lAHEH.ll 
HE  TURN 
END 
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SUBROUTINE  ToPOGRAFlF  ifiXOIM*NVC>lv.»N**NY«ri  »DtLF  .NC  *X  INUItS  * 

X Y INCHES. IMAGE, 

1 XM I N *DX  » XF  URMAT  .XLAbEL  *NC.X  . Yl. IN»l>Y . YMJRHAT  . VLAotL tNC. Y ) 

TOPOGHAF  DRAWS  A TOPOGRAPHICAL  PLOT  Op  lilt  VALUES  IN  THE  ARRAY  FINX.NYI. 
USING  5IJPROUT  INK  CONTOUR.  CONTuilK  LINES  WILL  lit  DRAWN  KUU  NL  VALUES  OF  K. 
AT  FI.  F1-0ELE  » F 1- I NC- 1 1 OtLF  . 

OR.  IF  DELF • THF  POUT  INI.  WILL  CALCULATE  THE  MAXIMUM  AND  MNjI’Uf'  VALUES 
OF  F AND  DRAW  NC  CONTOURS  BF.TWEEN  THEM. 

X I NCHt'S  ■ C.HAPH  LLNGTH  IN  INCHLS*  Y I FiCHEiS  ■ GRAPH  liLIUhl  IN  INOED. 

IMAGE  « NX«NY  STORAGE  LOCATIONS  FOR  tJSF  PY  CONTOUR. 

XM I N • 1ST  X-VALUE.  DX  ■ X-VALUF  INCREMENT.  XFORMAT  » FORMAT  FUR  X-VALUES. 
XLAHEL  ■ NCX  HOLLERITH  CIimKALTLRS  To  I.At'LL  TnE  X-AXIS. 

YMIN.  DY.  YFuRMAT  . YLAOEL.  .'ICY  PROVIDE  CORRESPONDING  VALUES  FOR  Y-AxlS. 


COMMON  /PPM/  TENSI6A.6AI.  I>1  IGA  I . UJMftMl*  oiiASi-  l ( GA  I . Pl.A.NF  / I NR)  » 
• CONMlFlt  It  1 * "ZI'Jli  SV1K.A  I.  f.  VU I GA  ) . PARM  ( 20)  » D I f TOR  1 20 . { I 
DIMENSION  F INXMM  .MYP 1" ) 
tr  INCX.GT.nl  GO  TO  It. 

DX  T *DX»  I NX- 1 1 / It  ,U 
DY  T »0 Y*  I NY- 11/1.  ,1. 

GO  TO  Hi- 
lo DXT ■ Au.’jl  XwlNI  /{>• 

DY  T *AU5  I Ylil  Nl  FS» 

NCX. -A 

PC  CONTINUE 

IF  IULLF  .Nl  ..i  I UU  TO  A 

F " I N * K * ' A X ■ F I 1 ,1  | 
no  s l ■ l . n x 
DO  f;  J t.1.  MY 

IF  IF  II  .J)  - F'llNI  H.S.-I 
? r • ■ i \*  » f 1 1 .j) 

GO  TO  !i 

I II  IF  I | .01  - I "AX  I • 

It  F.MAX  " F I I . Jl 

S tOF.TII.Ui 


I'l  . 1 1 ‘ > A Y. 

r l i vi  l ■ 
tiO  TO  7 


I MIN)  /I'll  - II 
I "AX 


%■ 


v'V 

A? 

'H& 

s y-A 

\j  >* . 

t. 


**• 

✓ 


F L I 01  l.  " 11 
111  « 01  LI 


/ NCl.NC-l 

CALL  a X 1 rt  I . . . XI  A|t|  |.  . .(.  A « 1 . » . , . - . 1 . . A.  1 1.  .11 A I . Al  . .1  • I.  I I 
(.ALL  Axis  It.  .YI.Ai.l  I. » ,C  V ,!■  |,  .Y  ilN.uvt  , H t t ||) 

F L VL  *1  1.1  VI  l 
I'  * U 

IF  I mil, I X l‘K  HI  1 . • I.ET.  .-II  '■(  T . ‘ 

X * 1 

r all  i-i  oti s. »s. .-M 

*.  t-  (t.uriFii'i 

In;  1 I « l . NC  1 
I II  vi  E»i < ■'!'  1 1 '.I  I I '‘I  L Vi. 

IF  IR. El, '.'I  CALL  ( ONTi.nR  1 1 • xo.  . -■■»  " 1 .07.1  I'l.  Vi  E. 

IX  I NC  > It  !, « Y I NC  III  IVAt.l  1 

I I (.■,10.11  (All.  I O’.  1 1 - 1 1 .NX  I - 1.  .’  Y I ; ,.  ,x  .'  t . ■ i i vt  I. . 

IX  |f, (.HI  . V I NC  > .1  '.»  IP.. I.  I 
I (C'NTlNUI 

Yl-'"'  , 

UU  l-'l. '-T  IIP.  *YI>,-  < ) 

PI  TllHf 
I VI , 
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SUBROUTINE  KONTUR  I F »NXD  lM.NYDlll.NX  »NY  .FLt'VFL  »X  INCHES  *V»  *•  ,‘HES. 

1 I MAGE  ) 

DIMENSION  F I NXD I M »NYD I M ) » IN.AGL  tNXDIM.NYDIM  | 

CONTOUR  DOES  INVERSE  LUHlHLt  INTERPOLATION  ON  A 2«D I MENS tGHAL  ARRAY.  FIX»Y| 
WHEN  CALLED  WITH  A GIVEN  FLLVEL  VALUE « IT  RFTuRNS  AFTER  HAVING  PLOTTED 
A SET  OF  CONTOUR  LINES*  WHEKt  F « FLLVEL*  ON  A GRAPH  XlNCHES  LONG* 

AM)  Y INCHES  HIGH. 

XF  AC  TOR  « XlN(  HI  S/INX-l  I 
Y FAC  TOR  « YINCHES/IMY-II 
X0«XlNCHES/2. 

YO«Y INCHES/ 2* 

LOAD  IMAGE  ApRAY 

DO  2 1 Y ■ 1 *NY 
DO  2 I X **  1 .NX 

IF  (FI  IX.IYI.Gt.ELI  VLL I GO  TO  1 
I M AGk I IX* I Yl  ■ -1 

GO  TO  2 

1 IMAGE < 1 X * I Y ) ■ 1 
? CONTINUE 

SCAN  IMAGE  EUR  THE  1ST  POINT  "f  A REGION 
IYSTANT  . 1 
1 no  i,  |Y«|YRTART.NY 
DO  a I X«  1 « fix 

IE  ( I VAGI'  I IX.  I Y|  ,i  0,1  I 00  To 
'>  CONTINUE 
RE  TURN 

LIFT  PEN  AF.p  OR  | N 0 T"  START  I ' 0 P,  I NT'*  AND  1ST  Tnl.  KH,U  ' H.UNI'i 

*•  I YSTART  « IY 

IT  ( |Y,,  i!.l  I f'O  TO  (» 

II  I IMAGE  I IX.|v-l)  i),)|  I.,..  TO  H 

CYO  • I IY-1-IM  IX.  I v)-l  1.1  Vt  I I / (■  NX.  I v | _*  ( | » , | Y-|  | | |»Yt  „(  T t.i\ 
o>’  H-  / 
fi  OY'.i  » P 

7 CXU  • I IX- 1 I*  XT  ACTOR 
C.XpP“LXi.-Xr) 

C Y.'P'CY ,i-YD 

cal  i.  i'i.pt irx',  p *(  v.'P,;* i 
INPUT  « ? 

TO  if- 

Ti  T a it  I ati  |f;T!TI<  p.nuNi  Ary 
Y INOUT  » I 

CX->  ■ (lx-?  I«XFA<  TOR 

i i y»?+  i r 1 1 x-i . i y- 1 1 -T  il  pi  1. 1 / 1 1 1 1 x- 1 « i v-  j i -i  1 1 x - 1 * i v 1 1 1 * yt  r.  i. 
t xcp.’cx  -xn 
OYPP.CYi  -YD 
CALL  PLOT  If  Xi  P.CY  P , .1 1 

GO  TO 

NLlKT  D I Ml  C T I o':  Is  .'.li-i’Y',  Ui.O'  I Si  I'R  't.  l XT|  i *.A|  ■ io'.i..,’.  Y * 

AND  COUNT!  R-ri  OCX  W I M’h  i\’  I ; 1 1 1 I " A i ;<•  >i  I -A., Y 

I THE  INMI'I  OE  Till  HI  (-loN  Is  Tw  IN  i|.,nl  N tm  ■-*.  1 1 I |i.|.lh  D 

"US  I T I VI  X-(  ROSSINI, 


C X * I Ix-n  vxl  nCTOK 
II  ( I Y.l  ().F|Y|  l.'l  TO  1 | 

CY  • I 1 V-  1 * (I  I |X,  1Y|-I  u VI  L I / II  ( IX.  1 Yl 
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GO  TO  12 

11  CY  • (NY-1 l*YFACTOH 
CONT 1NUK 
CXP-CX-XO 
CYP*CY-YO 

CALL  PLOT  ( CXP »C YP  .2  I 

I F ICX.NF.CXO)  GO  TO  16 
YC.AP  » ABSFICY  - CYPI 
IF  ( YGAP,LT..Uwl)  (.0  TO  3 

16  IF  ( IX.E0.NX1  GO  10  AO 

IF  (K'llX  + l.lYI  - F'LLVLLI  AO.13.13 
13  IF  (IY.EO.NYI  GO  TO  1A 

IK  (F  I IX  + 1 1 1 Y+l  1 - FLLVLL  I lA.lb.16 
1A  IX  ■ IX  + 1 
GO  TO  in 
16  IX  ■ ix+l 
IY  ■ I Y+ 1 
GO  TO  2“ 

positive  y-cros-iing 

?.■  IK  I IX. I 0.1  1 GO  TO  ,’1 

CX  ■ I lx-l-tu  IX.  lYI  - KLK V-Ll/O  I |X»|Y|  - K'llX-l.lYIIMXIAL  Tim 
GO  TO  ?.! 

?\  CX  . I. 

2?  CY  « I lY-l IHYFACTOU 
c.xp»cx-xn 
CYPiCY-YO 

CALL  PLOT(rxP.CYI',2) 
i.v  ?(,  I » i x » n k 

IK  I M'Al'K  ll.lYI.LT.il  I'll  |(,  ;OJ 
26  I I I « I Y I « 

I K ( IY.K../.YI  1.0  TO  1 

II  1 1 I I X.  IY-1  1 - K‘L  ■ VI  LI  1>  .T  ( ' 

:> . li  i I x . -.11  i.f>  ro  ?'• 

II  t U I Y-l  » I Y-*  1 I i LI  vi  l.  I -'A  . ».  ■ 

?A  I Y . | Y 4 1 

(V-  T'l  7 
? A lx  . I x - 1 

| V i | Y + 1 

00  TJ 

Ki‘  (tA  T I Vi  X-<  K’O.V,  100 

».i  CX  • ( I X - 1 I + X K AC  TOO 
IK  I I Y . i 1 I I.o  TO  31 

CY  • I I Y-l-l  K I I X , lY  1 - loV  I.  I t I M I ' • I Y ) - II  I > » 1 V-  1 ) I | * Y,  ,■<  I,.  ■ 

r.<x  re  '^7 

1 I ( Y * 

7 Co  NT  I fill 
( xi'  + CX-VD 
f yl'.OY  -Y,5 

( Al  I.  I I ,r>TI<  Xl>«(  Yv.2  I 
1 1 ( cx.ip  .r  x I i.o  to  n 

II  (<  Y.'  -)„l  Y | 1-0  TO  ) 

■ii  1 1 I I x . I 1 . 1 ) i..  ■ To 

IT  ( l|  | V-|  , | Y 1 - I I I VI  I I .’A  . *'. 

w,  n ll<'.  >'o,  II  i.n  TO  t'l 

II  (I  I I V - l . IY-1  I - ' LI  VI  L I *1 
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00  TO  >0 
36  IX  • 1X-1 
IY  ■ 1Y-1 
GO  TO  60 

C NEGATIVE  Y-CROSSINO 


60 


61 


CY  • ( IY-1I*YFACT0R 
IF  ( IX. EG. NX)  GO  TO  61 
CX  • < 1X-1+IFI  IX* I Y ) - 
GO  TO  62 

CX  ■ <NX“1 ) *XFACTOR 


FLEVtL  I / ( F ( 1 X * I Y I 


62  CONTINUE 
CXP»CX-XD 
CYP*C V-YD 

CAUL  PL0T(CXP*CYP.2I 


6) 
46 
6 5 


( IY.E0.1 I GO  TO  10 

<FI1X»1Y-1>  - FUEVEU  30*61*43 

( 1X.F0.NX)  GO  TO  44 

(Fit  X+l  * I Y-  1 1 - PLKVl'.t  I 44  1 4 6 *4!i 


IF 
IF 
IF 
IF 

IY  « 1Y-1 
GO  TO  40 
I X + 1 
1Y-1 


IX 

IY 


00  TO  1 ' 


FNO 


F 1 1 X+l * l Y I > ) *XF  ACTOR 
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SUBROllT  t Nfc  CONTOUR  ( K *NXD  If.  ,D  YD  I'.' « NX  »NY  »FL  EVI.L  , X I NO  It  S » Y 1 NOIt  $ » 

1 IMAGE  I » 

DIMENSION  F ( NXII I M ,NVD  IN  I * I MAGH  NXI)  I M *NVIJ  I M ) 

C CONTOUR  OOFS  INVERSE  DOUBLE  IhTTROOLAT  ION  CN  A 2-D  I '•'t  NS  I UNAL  ARRAY,  F I X * V > 

C WHEN  CAU.ni  I'U  TH  A GIVEN  FLLVTL  VALtlki  IT  RETURNS  AFTER  HAVING  PLOTTED 

C A GET  OF  CONTOUR  LINES*  WtltKL  F » F LEVEL  * ON  a OKAPI  i XIUCnFS  LUND, 

C AND  YlNCuCS  H|(iH, 

XFACTOR  « XlNCHCS/INX-n 

yf actor  « yinchfs/ihy- J i 

C LOAD  ll'ACti:  ARRAY 

00  2 |Y«1*NY 
DO  ? I X«1 (NX 

IF  IFI  | Xi  lYl  *(.r*FLI  VI:  L I 00  10  1 
IMAf.f  I IXilYI  * “1 
GO  TO  2 

I IMAGE (IX, IV)  • 1 
? CONTINUE 

c.  scan  image  for  thi.  irt  point  of  a region 

1 YiiTART  • 1 
'<  DO  A |Y»|YSTAXT»NY 

DO  A I X • | , NX 

tr  I IMAM  1 IX,  I Yl. I 0.  I I DO  TO  >1 
A CONTI  NOT 
Ml.  TURN 

C LIFT  I'l'1  ADD  I MINT  T'  RTAOTIli  D.'ItT,  AM'  |.M  Im  M|d|.i“  K.'U'.O, 

‘>  lYSTARI  u |Y 

IF  I I V , I (,  , | I DO  TO  I, 

IF  I U'ADCI  IX* IY-1  I , '-'U  . I »D  tu  I' 

f.Y'  » I I Y- I -I  M I X . I V | -I  LI  VI  l I/O  l IX  , I > I -•  I l>  , I V- I I II  '■  Yl  At  T OK 

do  t ' : 

/,  f’Y-1  <■  •' 

/ ( x • > ••  I I X - I l*x»  t,(  to.: 

(■AIL  • I.  T l<  X ilY  i'l 

| Ml H 1 1 * ? 

«...  f , / 

r SFAvr  A *,  | *iM  Im'H /*;■ • A *< Y 

h |\.*»T  " 1 

( x ' * (hi-; ) * xr  i < tin* 

( Y‘  » I I W + IF  I I X-l  , I Y-  1 I -M  I VI  I I / I I I I > - 1 , t Y - | I “I  I I I • I Y , t|  ~YI  /.l  T <( 

I 'M  I I’l/'T  O X , V Y ','<1 

r,c  rr 

c VIM  D|I„CT|0N  |S  i'I.U  AY.'i  CI'UI.  IM  I OR  At,  IXTl.ft.-V.  M'M  OYi 

r A’:"  r ■ I IDT.-  '(•  ( I lift.:  I V |m<  A.  I \ T I - ’ . A I . : -i’l  Vi! r.  Y 

!'  ITM  1 1;  ■'  1 1 ' ‘ O*-  Till.  Ril.li'N  It.  .'.I'./V,  Tr  1 1 it  I < H I vt  111  ' I h I - ■ I : li  I I I 

C POMTIVI  X-(  INI!...  INI. 

10  f X • I I X«  ) I * XI  ACTOR 

It  I IV. I 'J.DYI  DU  TO  l l 

CV  1 llY-l*IIIIX.|Y|-l|IVIll/OI|X,l'l  * I I I M I V.  Ill  I vi.hK,, 

I*''  TO  1? 

lo  tau  i i.i -t  if  x.c  y,;m 

II  rr  » inV-i  i hyfactoi, 

II  ((  X..DI  ,1  X-  I (iO  Tu  K> 

YrAD  u Alls  I li  v - i.Yi  | 

IF  I Yi  Al>,|  T • • II  ( o To  i 
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16 

13 

1* 

IS 


! 

70 

\ 

V 

71 

r 

77 

;* 

?p 

73 

76 

2*. 


17 


WHITNEY,  MADER,  AND  ULRICH 


IF  ( I X * Kq. NX  I 00  TO  60 
IF  (FTIX+l.lYI  - FUFVtU  60  * 1 3 i 1 3 


IF  tlY.FO.KYI  00 
IF  (FI  1 X*1 » I Y t J | 
IX  » tX  + 1 
00  TO  10 
IX  » 1 X ♦ 1 
IY  « I Y+ 1 
00  TO  2v 


POSITIVF  Y-CR05S1M0 


TO  16 
- FimU 


16.U.1* 


if  i ix. rev. 1>  oo  TO  21  • 

CX  » I |X-1-<F(  IX.IYI  - F'LHVI;L1/IH  lX.IYI  « FI  IX-l.lYl)  l*XF ACTOH 
00  TO  27 
CX  « 0 

CY  ■ I I Y-ll * YFAC  TOR 
C Al  L PL0TICX.eY.7l 


on  26  l » i x . n x 

IF  I IMAt.FI  I . I YI.LT.l  I 00  10  2P 
I MAOR II  . I Y I ■ . 

IF  (|Y.U).hY|  t.O  TO  10 

IF  (FI  I X . I Y F 1 1 - FU'VI  L)  10.23. 2.) 

IF  I IX.r.O.  I I 00  TO  26 

IF  imX-l.lY  + ll  * FLt.VIU  26.2'i.2‘> 

I Y « I Y + 1 

r,n  jo  ?i 

IX  « lx-l 

1Y  « 1 Y * 1 

00  TO  3' 


> , , A T | V 1 »■■(  '<l)AA|MO 
I.X  . I |X-n»XFAl.Tl>K 
1 1 | |.  0.1  1 0(1  TO  7 1 

< > ■ I | Y.I-UII  X . I Y » - I 1. 1 > I I / 1 1 I I * . I Y ! 

on  TO  37 

CY  - 

I All  I’l  "TirX.fY.7l 


r i lx  . I i-l ) 1 )'  3i-  Ai.  i of 


IF  irx.ur.ex  1 on  to  ').< 


If 

1 1'  Y .I’ll 

,r  Y ) 00  To  3 

A A 

K 

( 1 x • i ■ :• 

a | 1 00  TO  ?' 

IF 

i r ( ix- 

1 .lYI  - FLI  VI  1. ) .o-.Jii  ,36 

JO 

II 

1 1 Y.l  (J 

■ II  >.0  TO  3 b 

II 

1 F ( IX- 

I.IY-)  i - n,  vi  i.i  vm  . »• . 

AH 

1* 

on 

■ ; x- 1 

T'i  3-i 

A 0 

IX 
1 Y 

1.0 

. Ix-| 
• IY-1 

TO  A' 

m i ' a r I vi 

v.f|.|l‘,'i  | F.1 1 

h.l 

C.  V 

• l I V- 

■ l 1 » Y 1 AC  T0I3 

1 r 

( l X »!  r 

'.'.X  » On  TO  <i  1 

c \ 
( 

» l 1 X - 1 + < M U • 1 V | - H 1 *.  t U t 1 1 

Tv-  ‘«P 

t*  \ 

( X 

- (NX* 

. 1 1 « XF  Al  Tun 

Ut> 

( Al 

V H '1 

1 1 ( X 1 1.  Y . 7 1 

ii 


I; 


1 


1 


1 


i] 


V 


J 


02 


i:  NRL  REPORT  8074 

l 

II F ( IV.tO.l)  00  TO  It) 

IF  (Ft  I X 1 1 Y“  1 1 - FLFVEU  30*43*43 
43  IF  I IX*EQ*NX ) GO  TO  44 

IF  (Ft  1 X*1 1 1 Y-l  ) -FLKVtUI  44*43*43 
44  IV  • IY-1 
fiO  TO  4U 
43  IX  • 1X4-1 
IY  • IY-1 
00  TO  Id 

FND 


i' 

s 

| 


I 


t 

i 

( 


B 
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SUBROUTINE  GETT  ICK  I DY, TICK  I 

C OETTICK  TAKES  ANY  GIVEN  NUMBER  OF  UNITS  PER  INCHi  DY,  FROM  SCALE  AND 
C GIVES  BACK  TWO  NEW  VALUES.  DY  AND  TICK.  FUN  USE  lN  Axis.  wMEjKt  THE  NEW 

C DY  ■ THE  NUMBER  OF  UNITS  PER  TICK.  AND  TICK  IS  BETWEEN  .«  AND  i INCHES. 

C THE  NEW  DY  ■ l, 2, OR  ! TIMES  ROME  POWER  OF  10. 

0  ■ AHSFIDYI 
CALL  NORMAL  1 0 . 1 EXP  I 
IF  ID-S.O)  2.1. A 

1 TICK  ■ 1 . u 
RETURN 

2 IF  ID-2. SI  3.3.  'I 
’ IF  tn-1.01  7.1.6 

A DY  ■ 1 0 « *1 ' 1 • ** 1 1 XP 
TICK  • lu./lf 
RETURN 

S OY  « >.»lU.»»ltXP 
TICK  ■ S.u/D 
Rl TURN 

A DY  • 2. • |D . ** I KXP 
TICK  • 2.U/I'. 

RETURN 

7 PRINT  l.O,  DY 

ID"  FORMAT  l/3'H  *«*»*MkUk  II,  GlITlLK  --  DY  « I.  1;) . 3/ I 
RETURN 
FND 

SMBRGUT  1 Ml*  NDRVAtj  lAKO.ItXPI 

C NORMAL  TAKES  ANY  NUMBER,  Al<(i,  AMD  NORMALITIES  IT,  II.,  CONVERTS  If 

C TO  THE  FOR",  AKC'.M.  tt#ll.XP,  WHERE  I ,LE . ARC.  ,L  T , 1 0. 

SION  * +1." 

II  XP  ■ 1 
IF  IARC.I  ft.S.l 
A SION  ■ -l.r 
AI?G  . -ARG 

1 IF  I A KB  - U'.SI  2, A. A 
.»  IF  I AMT  - 1.0 1 1,*j* 

1 ARM  ■ ARG* 1 P • P 
itxP  ■ irxp  - i 

no  to  l 

t>  ARC.  • ADCi/ 1 ?:«.■ 

IFXP  . IfXP  + 1 

no  to  l 

A AKG  » SI  OMARI, 

RETUHN 
END 
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SUBROUTINE  OETSCAUE  I Y»N .HE  1 GH I . 1M 1 N .0 i »K  » I ICK  * I FORMA  i ) 

CiE  T SCALE  OBTAINS  SCALING  PARAMETERS  FOR  I HE  N VALUES  IN  ARRAY  Y 
(OF  DIMENSION  N*K) . HEIGHT  IS  THE  GRAPH  HEIGHI  IN  INCHES. 

YM  IN  »DY » T 1 CK » AND  ! FORMAT  ARE  PROVIDED  BY  GE I SCALE  FOR  USE  BY  Axis 
YM IN  * THE  DATA  VALUE  Al  0 INCHES. 

DY  « THE  DATA  INCREMENT  / TICK.  AND 
TICK  ■ DISTANCE  BETWEEN  TICKS  IINCMESI. 
iformat  ■ a format  for  axis  labeling  to  fit  the  data  scaled 

ENTRY  GETSCALZ  INCLUDES  EERO  AMONG  THE  VALUEo  SCALED 
DATA  IN  ARRAY  Y MAY  BE  SCALED  BY  I HE  FOLLOWING  CONVERSION  10  INCHES 
Yl  I ) ■ (Y! I I -YMIN) /DY#T ICK 

ENTRYS  DOSCALE  ANU  UOSCALLA  CIiANLL  IIiu  VALUES  IN  ARRAY  y iO  INCHES 
ENTRY  SCALE  REPLACES  THE  SYSILM  ROuiIhE  SCALE  (U-.Iio  IFur.'IAi  HAKAl  .t  1 El<  I 

DIMENSION  YIN) 

IGATE«1 

YMIN  ■*  YMAX  « Y ( 1 ) 

GO  TO  1 

ENTRY  GFTSCAtZ 
I GATE *1 
YM  1 N ' YMAX  « 

GO  TO  1 
ENTRY  DOSCALE 
I C.ATfc -2 

YMlN*YMAX«Y(  1 I 
GO  TO  1 

ENTRY  DOSCALE?. 

IGATE-Z 
YM1N«YMAX*L 
GO  TO  1 

ENTRY  SCALE 
IGATU‘3 

YMIN*YMAX»YI  1 1 

1 M . N*K 
00  5 I ■ 1 t M.K 
IF  ( Y ( I I - YMIN I 2.5.3 

2 YMIN  * VIM 
GO  TO  3 

3 IF  ( YMAX  - Y ( I ) I 4.5.5 

4 YMAX  » Y(  1 I 

5 CON  MUM 

IF  1YNAX  - Y‘HN|  S.ft.7 

0 SCAl.t  A CONSTANT  ARRAY  if.  INCHIS 

5 DY  = TIC.'.  * 1.0 

IF  ( YMIN.FU*.'.  I GO  TO  1(. 

CALL  NORMAL  (Y\  IN,  ILXUI 
T V I N • 1 N T F ( Y M I M I 

IF  I TXMM.GT.YMIMI  T " 1 N ■Tl  - 1 N“  1 • 

YM  I N«TMI  M«  1 1!."1'  I r XP 

DY  IMCll»|)v.  ( 1'  . »4  I II  Xl’-t  I)  I /HI  I '.lit 

C.A.L  GMT  ICK  (UY.riCKl 

GO  TO  M 

C ROUND  OilwM  Y 'Ii'I  To  UY  SC  A I l 

l hY  = ( YMAX  -YMIN) /MIGHT 
CAI  L M.TTICK  (I'Y.T  1CK  ) 
n«DY 

CALL  NNRVAI.  11-ilHY) 
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call  normal  iymin.ivi 

X*  YM I N* 10 •**  t IY-IDY1 
TX* INTF ( X2D) *0 
tF  (TX.GT.XI  TX«TX-D 
YM1N«TX*K.**1DY 

C HOUND  YMAX  UP  TO  DY  SC  ALL  t .SO  MAX  Y vALob  PALL  IN  LA*i  ilU  oF'ACLl 

CALL  NORMAL  I YMAX . I MX) 

X'YMAX*M.*»(  IMX-IDYI 
Tx«  1 NTF  t XJP ) *0 
IF  HX.LT.X)  TX«TX+D 
YMAX»TX*10***IDY 

C ADJUST  TICK  LENGTH  FOR  WIDEN  KA.NUL  AND  ROUNDOFF 

PY INCH* I YMAX-YMINI /HE IGHI 
T 1 CK*P  Y/PY  I NCH  -,u  -v-OI 
IF  lYMlN.LT.'.l  YiM  I N*  I NTF  ( Y,  Ii,/.,Y)*1JY 

8 CONTINUE 

if  i iijATr.no. 1 1 oo  to  u: 

PO  9 1 * 1 1 !•’ . K 

9 Y|  1 ) "( Y I 1 1 -YMlN) /OYlNCM 
IF  ( lr.ATF.FO. 3 I RETURN 

C SELECT  APPHOPRIATU  FCH'AT  Fi  R THIS  DATA 


10  ENPT  KK»YMIN+l'Y*|NTFItiUuHI/  l 1 C<  1 
H 1 0»MAX  1 F t APS  ( Y.M  | N ) . AUS  ( END  I ILK)  I 
SWALl»DY 

CALL  NORMAL  I SMALL. I EXP  I 
CALL  NORMAL  Il'IG.NEXPI 
IF  UEXP.LT.-3)  00  10  1A 
IF  INEXP.C.f.A)  CiO  TO  1'* 

IF  ( IfcXP)  12*11.11 

C NO  DECIMAL 

11  IOEC-0 

I RANGE  "2  + NT.XP 

00  TO  13 

C WlTrl  DECIMAL 

12  I :m.  c ■- 1 f xp 

1 RANGE*  I PEC*? 

IF  I NtXP.Ml, .;  ) IRANOi.*  IRAf.o^  + Ni  ak 
C CONSTRUCT  FORMAT 

i»  tore i jFr*H*»P 

I RANOE • (RANGE *8**12 
t FORMAT  .4.111  ...  . OR  . I IM. L . O K . 1 H >\,\i  t 

return 

W.  I FORMAT* /.Hi; 8.1 
RETURN 
F NO 


66 


'WK1 


NRL  REPORT  8074 


SUBROUTINE  MULT1PLT  ( X • Y .N  iN Y S »K X *K Y * K T s * ALENG i H . i HE IGH , , wLABEL >N a 
l.YLABEL .NYC tLlNELABL. MARK. NMAHK»XMINvAL»DX.YMJNvAL»Dt  * » I ICK) 

C 

C MULTIPLT  WILL.  PLOI  NYS  LlNEs  OF  N POINTj  EACH.  ON  i HE  jAME  OHAPH. 

C ON  A LINEAR  OR  LOG  SCALE  i DEPENDING  UN  , HE  Eiy.k,  -.ED. 

C X VALUES  ARE  TAKEN  FROM  ARKAi  a Ai  Alii*  A!  1+Ka)  l+.N-i IKa) 

C THE  1ST  Y ARRAY  IS  SI  OREO  IN  ill),  i il+Ki  ).«»•« i i 1> i N-l ) K . ) 

C SUCCESSIVE  Y ARRAYS  BEGIN  KY5  LOCAlIONs  APARi. 

C THE  PLOT  WILL  RE  XLENGTH  INCHES  LONG  Bt  i HE  1 GH t INCHEo  HIGH. 

C EG.  THE  2ND  ARRAY  BEGINS  Ai  nl+Kisl.  ,Ht  LAw , A,  i . 1* . N ■ u-l  1*  , „ I 

C THE  X-AXIS  WILL  BE  LABELED  ..  I , H , HE  n„  CHAaAC  , E r\~  In  ..uABti.'.  Ai.D 

C MAY  BE  HOLLER!  i H Fukn.Aij,  EG.  BHEhEkG  ,,  ub  luCA.Ivn.  ..HEnE  .HE 

C CHARACTER  COOLS  ARE  SIOREO.) 

C L 1NELABL  IS  AN  ARRAY  OF  NYS  HOLLLRllH  wOHDs  WHICH  wILL  Lit  ua£D  iO 
C IDENTIFY  THE  NYS  LINES  WHENEVER  NYS  .Gi.  1 

C MARK  IS  AN  ARRAY  OF  NYS  IMLGLKS  DENOllNG  SPECIAL  jiMBOLs  t EG.  -1-N0IHING. 

C O-SQUARKSi  1 -OCTAGON....  ) 10  HE  DRAkn  un  EACH  lInE  K'vEh'  ii'iAra  Putins. 

C A NEGATIVE  NMARK  SUPPRESSES  I H£  CUNNECiInG  lInE  tlt,„EEi'i  - V.BULo. 

C 

C ENTRY  LOOPLOT  will  give  a LOG-LOG  PLOI 

C FNTRY  SEMILOGX  WILL  GIVE  A SEMILOG  PLOT  win.  X ON  THE  LOG  SCALE. 

C FNTRY  SEMILOGY  GIVES  A SLMII.OG  PLOI  wlIM  Y ON  I HE  LOG  SCALE. 

C FNTRY  LINFPLOT  GIVES  A LINEAR  PLU  i . 

C 

C OUTPUT  PARAMETERS  XMlNVAL.  ija.  YivIrvAli  Ui.  Glvfc  I i.F  uim»iA  . tun  LunLL  im.  1 nG 

C THE  SCALING  OF  THE  PLOT.  SO  THAI  A Pi.  Pis. I I MAi  HE  PLOmED  ON  iHt 

C SAME  GRAPH  BY  USING  THE  FOLLOWING  CONVERSION  10  INCHES— 

C XIINCHESl  ■ LOG  1 PX/XM 1 NVAL I *0X 

C Y I INCHES  I ■ LOG  I PY/ YM I NVAL I *0Y  FOR  LOGPlOI. 

c with  a linear  or  semi  log  ploi.  slhsiMuie  one  or  bo  i h of  iiit  follow  ing— 

C YllNO'FSl  ■ (PY  -YiI1NVAlI«Dt  Fur  lIrI>uui  Ai.D  uEilIluGa 

C XIINCHESl  • IPX  -XM1NVALI»DX  FOR  LlNO'LOl  AND  sO.lLuGt 

C 

C NB  IT  IS  EXPECTED  THAI  IHt  CALLING  PROGRAM  ..ILL  HAyt  AIRE  AI.)  i MADE 

C THE  INITIAL  PLOTS  CALL.  AND  WILL  ALSO  CALL  SIOI'PIUI  |(l  I CRMINA I E PLOllING. 

C ALLPLOTS  INCORPORATES  THE  PROCEDURES  OF  SCALE »L INE « A* I S . AND/OR  LOGPLOIS. 
c without  changing  the  contents  or  mi:  arrays  plumed. 

C ALLPLOTS  MUST  HE  USED  V,  I I M PLAIN  Ilg)  PLOmEH  GAPER.  Aj  1,  PROOgCEo 

c A VARIABLE  TICK  AND/UK  LOGCYCLE  LENtilH.  i SCALING  KUyilNEs  ou  1 i ABLE.  FUR 

C USE  WITH  AVAILAHLi  LINED  PLumING  PAPER  AI.  E LOGO  A PER  AND  -CALEiEN.I 

C ADDRESS  ANY  UUESTIONS  10  JEANNE  uLKlCM.  L«.  232P. 

C 

DIMENSION  Klllt  Y(l|,  L INELAI'L  ( I I 

DATA  If  XPSI2E-...7)  . I TENS  UI.  ■ . 1 ) . I I /.  tl.O«D  I . IM«.0?l 
C 

ENTRY  EOGf’lOT 
t GATE* 1 
GO  TO  1 

c 

FNTRY  SEMILOGY 
I GATE* 7 
(iO  TO  IV 
C 

f-NTRY  SM'lLOGX 
IOATE *2 
GO  To  1 

ENTRY  1.  INI.'PLOT 
ir.ATt  */. 
fiO  TO  19 
C 

C lOG-S<:AL!N<i  of  X 

C IS  I FIND  XMIN  AND  XNAX  VAl.lll 

c 

1 XMIN  * XIII 
X 'AX  a mil 
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M » N*KX 

00  5 I«2*M*XX 

IF  IXU)  -XMIN)  2*5*3 

2 XMIN  - XI t) 

GO  TO  5 

3 IF  ! XMAX  - XIII)  4.3*5 

4 XMAX  ■ XI M 

5 CONTINUE 

IF  (XMIN)  11*11*12 

11  I GATE*  I CiATE  + 2 
GO  TO  19 

12  CONTINUE 
C 

C DETERMINE  APPKQPR I AT  t SCALING 

C VALUES  ON  THE  GRAPH  RANGE  FROM  XMINVAL  Ui  u INCHEjI  iO  Iu*»MAaEaP 

C <AT  XLENGTH  INCHES). 

C NSCAI.ES  « NUMBER  OF  SCALES  (POWERS  Oh  ■ fc.NI  F OR  iMla  UA  i A jt  i . 

C FIPS  • iNCHtS  / SCALE 
C 

CALL  NORMAL IXNlN* I UASEXP I 
XMINVAL  • ll'.(.»*IUASEXP 
XMINLOG  u LOGFt XMINVAL I 
ARG  » XMAX 

CALL  NORMAL  I ARG *MAXEXP I 
IF  (ARG.EO.l.U)  GO  TO  6 
MAXEXP  ■ MAXEXP  + 1 

6 NSCALES  *■  MAXEXP  - I U AS K X P 
FIPS  • XLENGTH/NSCALES 

IF  IF  IPS. GT. 0.5  I GO  TO  7 
XLENGTH  « NSCALES 
FIPS  * 1 

7 DX  * FIPS/LOGFI lu.o) 

C 

C TO  scale  XARkAY  to  inches 

C XMINVAL  CORRESPONDS  TO  V INCHES*  X'-  lNLOG  • Lotil  aI.INvAL ) 

C DX  « INCHES  / DATA- LOG-UNIT!.  * Ii  • DA  Ci/nivl-u  i lvOuhIis  io  IkCHEj 
C LOG!  101  * ONE  SCALE  IN  DATA-UOG-UN I I S 

C LOG!  10)  »DX  " FIPS  ■ ONE  SCALE  IN  INCH!..* 

C 

C DRAW  AND  LAHt. L X-AXIS 

C 

YDN  • -.OH 
YUDN  ■ -.2 
EXPLOCY  ■ »»3>- 
TKNl.OCY  ■ -.4,: 

CALL  PLOT  (o*u*1l 
IFXP  • IMASf  XP 
DO  10  In  | ZERO* NSC. ALES 
X LOCI!  « XLOC  « I » F I o.'i 
CALI  PLOTIXLOC*.’*2) 

CALL  F L1 >T I XL  DC  * YDDN  *2  I 
TENLOC  « XLOC-.15 

CALL  SYMBOL  ( TENLOC  * TLf.LUCY*  1 1 NS  ILL  »/) 

CALL  NU'TUt  R ( XLUC  * L XPLGC.Y  * I XPSl/.l  * Ii  AP*o.,Hil  (I 

CALL  PLOTIXLOC.  .31 

IF  ( I .LU.NSCALI SI  GO  To  U 

ILXP  « IEXP  + 1 

no  9 J»2.9 

XLOC  « XLOCi.  ♦ LCWU  LpATI-  I Jl  l«HX 
CALL  PLOT  ( XI.OC  *•■  *2  ) 

CALL  PLOT ( XI  DC  lYliM  *2  I 
IF  (FIPS.LT.  1 • *■  I (>(>  To  9 
CALL  NUMPI  R I XLOC  * VPl.fi  1 1 XPS  IV i J » . »P»(  I 1 I 
9 CAUL  PLOT  ( XLOC  *..  *1  I 
10  CONTINUE 

XLAHLOC  • Xl.LNGTH  + ;.5 
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CAUL  SVMBOl.  I XLAbLOC  «>-  * TtNS  UE  • XLAbEl,  »0  .NX  ) 

GO  TO  I20.JBI.  IGATE 
C 

C FOR  NON-LOG  SCALING  OF  X 

19  CALL  GETSCALE  t X ,N  .XLENGTh.XI'"  IN  »0X  |KX  .X  ! ICK  I XFL/KM  ) 

CALL  AXIS  ICL.XLAliEL.-NX.XLtNGTH.OlXI  1CK  .xv lN*DX .XKLmm) 
GO  TO  (20.2G.2C3HI  . IGATE 

LOG-SCALING  OF  Y«S 

20  YtENGTM*  YHK  I fil  If 
IF  IYrtfclGHT.GT.lt  . I YLINOT>i«lo. 

21  YMlN  • Ylll 
YMAX  « Ylll 
f,YS.‘IYL»KYS 
GO  21  I 1 ■ 1 .MY ft  . KY  S 
M. | 1 .MKY-1 

no  ?<>  i-ii.m.xy 

|t‘  (Ylll  -YMU1I  22.24*21 

22  YM  I N * Ylll 
GO  TO  24 

23  IF  I YMAX  - Y(II)  24.25.24 

24  YMAX  ■ Ylll 
24  C0NT1NUF 

IF  (YMINI  14.14.14 

14  I GATE ■ I GATE* 1 
GO  TO  3 a 

15  CONTlNlir 

OflfcHMINF  APPHOPN  I ATi  SCALING 

CALL  NORMAL  I YM|N»  IHASI  XP  I 
YMlNYAL  » 1 o»u«*  I llAbL  AP 
VM 1 NLt-HJ  * LOGF  I YM I MVAL  I 
ARG  ■ YMAX 

CALL  NORMAL  I AWh.'IAXFXM  I 
IF  | AkLi.FL.  1»'  I GW  |u  2t> 

MAXES!'  • KAXF  XP  . 1 
j I,  NSC  ALES  » MAX  fc  XP  - I It  A;, I XP 
F I P S * YLl  NGTH/VRCAU 
2 7 nv  ■ FI  Pf./LOGF  II  .f  I 

c liRAp  ANO  l.AI'l  L Y-AXIS 

c 

XOUT  * -.!  « 

XOUTT  « -.?l‘ 

XOUTN  » -.IS 
TENLOCX  « -.4o 
EXPLOCX  ■ -.30 
CALL  PLOT  1 1.  . 3 I 

DP  30  1 « 1 2.E  HO  .NBC  A L T It 
YLOCO  . YLOC  « l*F IPS 
I F XP  • I'tASI  XP  + 1 
CALI.  PLOT  ID.YL0C.2l 
CALL  PLOT  | XOUTT « YLOC  .2  I 

CALL  SYMBOL  (TIMLOCX.  YLOC.  TEri.'  IFI  » 2"1  * ' • > 

EXPLOCY  . YLOC  + •'  7 , . 

CALL  MJM.PfcH  ILXPLOCx.  I.XPl.tlY,  iXP'.lAi..  lt>P.  > • • " 1 1 

CALI.  PLOT  It  . YLOC  i 3 1 
IF  I I .It),  NSC  AL I SI  GO  TO  )■ 

DO  29  J ■?»•< 

YLOC  1 YLOC  t'  + Lt!GF  1 1 L'-  AT  M J ) I 1 ' ' Y 
CALL  PLOT  It.  Yl.OC.  21 
CALL  PLOT  IXUlT.  YLt  C . 21 

u ir  ipg.lt. 2.  ■)  l>o  10  <"1 

CALL  MOOfTH  IXIHITN.  YLOC.  IXP:.I/I»  J.  ■'  * 2011) 
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29  CAUL  PLOT  (0,  YLOC . 31 

30  CONTINUE 

YLABLOC  • YL£NGTH*«S  - TENS  1 2L#NY* ( 3 » /7» I 
CALL  SYMBOL  <-•  9 » YLABLOC  iTENSlZt  tYLAtik  LlYC.  iNYC) 

00  TO  40 

C FOR  NON-LOO  SCALING  OF  Y 

3fl  YLENGTH*  YHF  I GUT 

IF  I YhEIOHTtOT.lrt.  ) YLf.MGTH»10« 

C 39  CALL  GETSCALE  I Y iN*NY5  * YLENGTH,  Yl  ,IN  *DY  «K  Y , r l I lx » tFUKi*,) 

39  CALL  6ETSCAL2  t Y*N#NYb  »YLtNGilt»  T i*.  1 ih > U i im  i i i ILm  i Fwititil 

CALL  AXIS  I 0 tc  tYLAUEL  iNyC  hLENGi  lli9„,  * i i ILK  * i .'-UN  «0  , • , FORM) 

DRAW  THE  PLOT  TO  SCALE  t LEAVING  THt  CUM  I tN  IS  OF  ARRAYS  X ANI.1  Y UNCHANGED 

40  CONTINUE 

43  GO  TO  <9i>»60»7,i»Bi;l  t IGATE 

C LOG-LOG  SCALING 

30  DO  32  NY*1 iNYS 

1 1 * ( NY- 1 )*KYS+1 
X IN- 1 LOOP  < X ( 1 ) l-XP I NLOG ) *DX 

Y tN*  (LOOP  1 Y ( It)  ) -YM  INLOG  I *DY 
CALL  PLOT  t X I !% » Y I N * 3 ) 

DO  91  l*2»N 
i x- 1 t-n«Kx*i 
IY«I I-l l*KY+l l 

X IN*  I LOIiFI X I IXI  l-XM INLOO )*0X 

Y IN*  ( LOGF ( Y I I Y I I-YMI NLOfi I *0Y 

31  CALL  PLOT  ( X I N t Y I N * 2 ) 

32  IF  (NYS.GT.t)  CALL  SYllDOL  I X In  i Y I.,  »H*L  INEL  AHL  t NY  I »i;  .#  I 
RETURN 

C SEMILOGX  SCALING 

6U  DY*YT  ICK/DY 

DO  62  NY*liNYS 
I 1 • t NV- 1 l»KYS*l 
X IN*  ( LOGF  I X I 1 I ) -XT’ t NLOG)  *UX 
Y IN* ( Yl  I 1 I-YM1NI »DY 
CALI  PLOT  IXIN.Y1N.3I 
DO  41  I *2  tN 
I X “ ( 1-1  l*NX+l 
I Y * ( 1 - 1 ) * NY  + I 1 

x in*{logf(xi  ix)  i-xvinldc.iM'X 

Y I N * ( Y I lYl-Yf’lM  *I’Y 
M CALI.  PLOT  <XlN»Y|N*2> 

62  IK  I N Y S « Ci  T • 1 ) CALL  SYIm«>L  I X I « Y I *M  tL  1 Ml  L AhL  ( NV  ) ( c,  • U I 

wf turn 

C 5LMILUGY  SCALING 

70  PX*X  T ICX/DX 

DO  7?  NY* 1 iNYS 
I 1 • l M V-  1 ) * K Y S ♦ 1 
X IN*  I X (1  I-XMINI  *I)X 

Y IN*  I LOGF  ( Y I I 1 I ) -Y," 1 INLm.  ) « I ■ Y 
CALI  PLOT  (XlNiYIM.Y) 

DO  71  1*2,'! 

IX*( 1-1 l*KX+l 
I Y* I 1-1  | »NY+  1 1 
X IN*  I X I I X )-XN'lNI»DX 

Y IN*  I LOGF  ( Y ( I Y I l-YMM.uli  1*1.9 

71  (All.  PLOT  Ix1NiY|N«2I 

7?  I)  INYG.C.T.ll  CALL  NY.-iluL  ( X 1 1. 1 r I >.  »D  »L  I Hi  I.  Apt  l .( t I * I H : 

we  turn 
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UM'a.x  r.rAU>(, 
nr  nx-xT  ICY /r>x 
[iv»yi  icx.yiY 
!>')  8?  NY«l»fiYS 
I \«iny-i  )*h'.y;,+i 
X IN*  (XII  i-XVP:)*nx 
Y IN*  t Y I n l-YI-.INI  »I'V 
(.ALL  I'LilT  1X11  • YIN  *11 
110  <U  l»2.H 
|X«I 1-1  )»KX+1 
lY«l  1-1  HKY+lt 

x I n ■ I x l lxl»xvtni*nx 

Y 1 « I Y I I Y ) -Y  'l  J>! ) *i)¥ 

PI  CALL  I’LOT  <XIA:»YII,«?) 

82  IF  (NYS.C.T.ll  CALL  SYM;f>L  I X 1 li  »Y  1 " »t(»L  I NLL  Ai-L  I NY  1 »o  t b ) 
mi  turn 
l-NT) 
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SUBROUTINE  FASTFOUH<A,M.INV»SilFS*lFERR>  EASTF  1 

FASTFOuR  ...  DISCRETE  FOURIER  TRANSFORM  ...  FORIRAN  63 
INPUT  PARAMETERS  TO  BE  SET  BY  USER  BEFORE  EN  I EHlN(i  FAST  FOUR 

A IS  A 3-UIMENSIONAL'  AHHAY  OF  COMPLEX  COEFFICIENTS.  OF 
DIMENSION  lN(l)tNI2)iN(3)l« 

THE  A'S  ARE  STORED  WITtl  REAL  PART  OF  Atll.12.13)  IN  THE 
LOCATION  WITH  INDEX  2*<I3*N1«N2  ♦ I2*N1  * ll)  + l ANU  THE 
IMAOINARY  PART  IN  THE  LOCATION  IMMEDIATELY  FOLLOWING. 

IF  THE  FOURIER  SERIES  IS  REQUESTED.  ARRAY  A IS  REPLACED 

BY  .... 

XI J1  ,J2.J3)«SUM  A(Kl«X2.K3>«Wl*»IKl»JI|ttW2«»(K7*J2)*W3»*U3*J3l 
SUMMFD  OVER  Ki*2*  Ni-i*  K<t-.n  M.-.*  X j 0.  Nj-. 

WHERE  W!  » N I -TH  ROOT  OF  UNITY. 

Nl »Z**MI I I IS  THE  NO.  OF  POINTS  IN  THE  l*TH  DIMENSION. 

THE  DIMENSION  OF  A IN  Tht  CALLING  PROUHAM  SHOULD  HE  TWICE 
THE  NUMBER  OF  COMPLEX  ELEMENTS  IN  THE  LARGEST  A ARRAY  TO 
BE  PROCESSED. 

THE  COMPLEX  X'S  AKL  STOKED  IN  THE  SAML  MANNER  AS  A. 

IF  THE  FOURIER  TRANSFORM  IS  REQUESTED.  The  ARGUMENT  A 
IS  TAKEN  TO  BE  X AND  IS  REPLACED  BY  THE  ARRAY  A SATISFY- 
ING THE  FOURIER  SERIES. 

LFT  MT-MAXIK1, M2. l.<3)-2,  NT  ■2*<MT  * WITH  M l-  K I NO  Tub 
M GIVEN  WHEN  THE  TABLLS  ARE  SET. 

SI J I «S INI J»PI /2  »NTI  .J«l.  NT-5 

lNVIJ  + ll»WORi)  CONTAINING  BITS  OE  J IN  INVERTED  ORDER  IN  ITS 
RIGHTMOST  NT  BIT  POSITIONS.  FOR  J»0,  NT-1 

TO  SET  UP  SIN  AND  I NV  TABLES  ... 

CAUL  SE  TUP  I A »M » I NV  »S  *L . I E LRK ) 

ONE  NEED  NOT  REPEAT  THL  CALI  To  SETUP  IF  ONE  Out!.  TUT 
CHANGE  the:  MAXIMUM  M« 


DIMENSION  Al  l ) »'  1 31  . INVI  1 1 »M  i 1 <31  *KP(  3 1 

t*  A S T h 

2 

EQUIVALENCE  1 N1  *N  1 1 1 1 » 1 N2  *i.  1 2 I i » 1 .3 *ul  3 m 

lh 

3 

c. 

CALL  WITH  1ES  • +1  1 OR  lUURli.R 

IN  1 1. 

c 

CALL  WITH  IE'S  ■ -1  FOR  Kwok 2 i. K 

T .<»>  N ‘ih  \':i<  • 

c 

f 

CALL  RLT^MI  A lfiViii • II  l «N  1 TO  M.  T T*\out.’>  Lf* 

i i N ANt  > l K V 

t. 

c 

IFEPR  • i,  IvEII  N AkUIIi-.I  .,1 S ,,  ,’m| 

UK 

c 

c 

IFFRR  • 1 WHEN  THERE  IS  Ap,  | ivu 

1 1 ; C Ai.  1. 1 1 

DATA  (SlNAS  » • 7v  7 1 . bTlt  1 2 ) . 1 P 1 * 3 . 1*.  1 

v fc  •>  1 h } 

I A.',  II 

u 

1 u 

1 F l IF;,)  12,13 

► AMI 

3 

1 2 

E'TT  • XMAXoF  !.•'  1 1 1 .,'U2  1 ,1  1 31  1 -2 

1 Ah  II 

b 

IFIMTT  .Lt.  MT)  U.IJ 

y a:,  i k 

7 

1 3 

lEtRR  « 1 » GU  TO  9VVV 

» a*  ;•  1 1 

It 

1 A 

IFE.RR  • IT 

I AMI 

V 

Ml  ■ ,M(  l 1 i M2  ■ .( 2)  J 

■ 3 ■ . : t j ) 

» Ah  71 

1 i 

Ml  ■ 2»»M1  S N2  * ?*»M?  T- 

\i  5 ■ ? * m f 1 » 

I- AM  I 

1 1 

IFIIFSI  lf>,9VV'»,2,j 

I-  A:,  1 1- 

12 

C 

IALCULATI  TKAN'.FUi  ■ ...  Blil.AU 

A ti  V ( < /<•  J1  • ( A ) / • . 

1 6 

N To T « N1*N2*N3  i Hi  « 1 . /I  1 E 1 . 1 1 v II  . i , 

1 1 A\  1 1 

1 ' 

00  1H  I • 2.NTOT2.2 

» Mil  1 I- 

1<* 
1 ') 

Al  1-1  1 • A ( 1 - 1 | M'.'l 

1 »A  , ] f 

1 II 

All)  « -Al  1 l»EN 

I AMI 

Kl 

2 n 

NP(  1 1 * 111  +N  1 

I-  AMI 

1 7 

NIM  2 | ■ HIM  1 )»N? 

I A ' , 1 1 

IP 
1 u 

NIU  3)  ■ NP  1 2 1 *N3 

1 Ah  11 

HERE  BEGINS  I HI  LoiiH  mi  Til  Ti-  . I i.|. 
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DO  230  ID  ■ 1*3 

IL  ■ NPI3>  - NPI ID)  * 111  ■ 

Ml  ■ Ml  1 0 > 

IF  I .NOT.  MI  I 230*30 
10 1 F ■ KBIT  • NPI 101 
IF  (Ml  .EG.  2»IMl/2l ) 60*40 

Ml  IS  ODD.  DO  L ■ 1 CASE 
KBIT  » KBIT/2  * KL  ■ 

DO  SO  I • 1.IL1.IDIF 
KLAST  ■ KL  + l 
00  SO  K • I *KLAST  *2 
KD  * K * KBIT 

ONE  STtP  WITH  L* 1 * J-0 
AIK)  ■ AIK)  * AIKD)  * AIKO 
HEAL  PART 

T » AIKD  ) S AIKO  I ■ AIK 
IMAG  PART 

T ■ AIKD+1)  » AlKD+ll  » AlK+ll-T 
CONTINUE 

IF  (MI  .EO.  1)  2S‘i * 52 
LFIRST  « 3 * JLAST  • 1 

DEF  ...  JLAS1  ■ 2*»<L-2)  -1 
LFIRST  ■ 2 S JLAST  • 0 

DO  24u  L ■ LUHST  . Ml*  2 


IL1  * IL+1 


KL  ■ KBIT  - 2 


AIKD)  • AIK) 


A I K * 1 1 


GO  TO  70 


Kb  IT/* 


FASTF  20 
KASTF  21 
FASTF  22 
FASTF  23 
FASTF  24 
FASTF  23 

FASTF  26 
FASTF  27 
FASTF  28 
FASTF  29 
FASTF  30 


• AIK  ) + T PASTE  31 

» vlK+1  ) * T FASTF  32 
FASTF  33 
FASTF  34 
FASTF  33 

FASTF  36 
8 AS  TP  3T 
Kl.  • KH1T-2FASTF  311 


00  FOK  jmi 

DO  60  1 • 1.IL1.1DIF 

KLAST  « I+KL 

DO  80  K ■ I *KLA5T*2 

Kl  ■ K+KBIT  * K 2 * Kl+KUlT 


K3  » K2+KUI t 


FASTF  39 
FASTF  40 
FASTF  41 
FASTF  42 


DO  TWO  STEPS  WITH  J-0 

AIK)  ■ AIK)  ♦ AIK2I  « AIK21  « AIK)  - /•  I K.  2 ) 

AIKD  ■ AIKl)  + A I K3 ) * AIK3I  * AIKD-  AIK3I 
AIK  I ■ AIK  I ♦ A(Kl)  . AIKl)  . AIK  ) - AIKD 
AIK2I  ■ A I K2 ) ♦ A ( K 3 1 * I « MK3)  * AIK2)  - A|k.))*1 


FIRST  STEP  RtAU  PART 
T ■ A I K2 ) S A I K2 1 « A ( K l-T  * 

T ■ A I K 3 ) 1 AK3I  * «IK  1 l-T  V 

FIRST  STF.P  IF  AM  part 
T • AIK2+1I  4 AIK2+1I  ■ AU  + D-T 

T • A ( K3+1I  t AIKJ+ll  * AlKl+D-T 

SECOND  STEP  REAL  I'ArT 
T.AIKll  * AIK  1)  » « I K )-T  1- 

T ■ AIK3I  1 MKT)  ■ A I K2  I + A I K.  1+  1 ) 
SECOND  STEP  iMAii  PART 

A I K 3 + 1 ) ■ A IK. 2+  D -T  » 

T « AIKI-D  t Alsl-.ll  ■ AD.  + II-I 
CONTINUE 

IFI.FIOT.  Jl  AST  I 233*62 
JJ  • JJDIF  + 1 
DO  FUR  J*1 
I LA  ST  * Il+JJ 
Kl  H5  I . J J*  HAST  * I OIF 
KLAST  • r.l  ♦ I 
IK)  83  K • I .KLAST  *2 

Kl  • K ♦ KBIT  7 K2  « El  + Mi  I T 


AIK  ) * A(K  l+T 
AIKl)  " A | K D ♦ T 


AIK  ♦ 1 ) 
A I K 1 * 1 I 


AIK  +1I+T 
A I k l * 1 I * I 


AIK  ) • A(k  I + 1 

AIK.)  - '•  IK2l-AIKJ.il 

A I K2+  1 I ■ « 1 1-  2*1  l + T 
•*  A I K ♦ 1 I * A I K ♦ 11  + T 


« K l + Kli  I T 


LETT  I Nil  W * F XI’  I 1 *P  1 /4 '-  » w2  « W**2»  V'  • """l 

AIM  • AIK)  + A I K 2 ) . I * A I * 2 I • Mil  - M K 2 I * I 

AIKD  ■ A(Kl)*W  + A(KJ)»vT,  a (1*1  « A I K D * M - MMI‘|-,I 

AIK)  ■ AIKl  + AIKD  * AIK1I  » AIK)  - ADII 

A ( K 2 I ■ A I K2  ) + A(K.n*l*  A I k .» I . AIK.)  - M K.)  ) » I 


FASTF  AJ 
FASTF  AS 

FASTF  AS 
FA*.,  TF  AA 

MvSlI  4 7 
I-  A S I F AM 

I AS  It  AS 
)AS1F  31' 
FA.STF  31 
FASTF  32 
FASTI  3 + 

I AS  IF  SA 
I ns  1 1 3 3 

F A S If  flA 
FASTF  ! 
F a ; 1 1 1-  3 o 
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R • AIK2+1I  % T ■ A IK2 ) 

FASTF 

39 

A ( K2 ) * AIK)  * R * AIK)  • AIK)  - R 

FASTF 

60 

f • 

AIK2+1I  ■ A I K*1  1 - T » AIK+1)  ■ AIK+ll  ♦ T 

FASTF 

61 

V. 

AWR  ■ AIK1I  - AIKl+ll  * AWI  ■ AIKl+l)  +AIK1) 

FASTr 

62 

R • AIKS)  + AIK3+1)  4 T • AIK3I  - AU3+1) 

FASTF 

63 

AIKS)  « (AWR  + RI»SIN4S  t AIK'  .■  IAWR  - RI+51N4S 

FASTF 

64 

p 

AIK3+1I  • 1 AW I - T1+SIN4S  $ Alll.  J)  » (AWI  + T)»SIN4S 

FASTF 

63 

i. 

T • A ( K1 ) 4 A (K  1 1 ■ AIK  1 -T  $ AIK)  ■ A 1 K ) ♦ T 

FASTF 

66 

r 

T « AIKl+l)  V AU1  + D ■ AIK+ll-T  » A(K+1)  • AtK+ll+T 

FASTF 

6T 

V 

R » AIK3+1)  t T ■ AIKS) 

FASTF 

68 

A 1 K3 ) • AIK2I+R  4 AIK2)  « AIK2I-R 

FASTF 

69 

AIK3+H  ■ AIK2+1I-T  4 AIK2+1)  ■ AIK2+1I+T 

FASTF 

70 

8 3 

CONTINUE 

FASTF 

71 

i 

IFIJLAST  *IE*  1)  23S»90 

FASTF 

72 

ri 

90 

JJ«JJ+JJOlF 

FASTF 

73 

C 

NOW  no  THE  REMAINING  J'S 

f; 

DO  230  J ■ ?»JLAST 

FASTF 

74 

i 

c 

KETCH  W'S 

l 

c 

DEF  Wl»tXP(  l*Pl/4)+*tNV( JI/MT  * W2“W1**2  • w3«Wl**'J 

i 

98 

1 * lNVtJ+1) 

FASTF 

73 

f> 

98 

Wlt»8(NT-ll  '*  K 1 2 * S I 1 ) 

FASTF 

76 

y 

12  » I + 1 3 12C  • NT-12 

FASTF 

77 

i 

IFII2CI  1 2o » 1 lu»  luO 

FASTF 

78 

c 

C 

12  IS  IN  MKST  QUADRANT 

f. 

100 

W21  « 3I12CI  4 W22  « S(  1 2 1 4 00  TO  130 

FASTF 

79 

k- 

i„ 

C 

12  IS  PI/2 

(■■ 

s 

no 

W21  • 0.  4 W22  • 1.  A 00  TO  130 

FASTF 

80 

!'■ 

S' 

c 

12  IS  IN  SECOND  QUADRANT 

| 

no 

I2CC  ■ 12C  + NT  » 1 2C  ■ w 1 2C 

FASTF 

HI 

£' 

W21  ■ -S( I2C)  * W22  “ SI  I2CCI 

FASTF 

H2 

f 

no 

IS  « M2  4 ISC  ■ NT- 1 3 

FASTF 

H3 

li 

IMI3CI  160  * ISO  * 14" 

FASTI- 

B4 

i 

c 

13  IN  HIRST  OUAl’RANT 

no 

W 3 1 ■ 5II3CI  t W32  » SI  131  4 (.0  TO  200 

FASTF 

83 

i 

c 

13  • PI/2 

ISO 

W31  ■ 0,  % W 32  ■ 1*  4 00  Tu  200 

FASTF 

86 

c 

11  IN  2ND  OH  3RD  QUADRANT 

no 

I ICC  ■ ISC  + NT 

FASTF 

8 7 

IFII30CI  190  * 1 Hu  * 1 7 1 

FASTF 

HH 

r 

11  IN  2ND  QUADRANT 

1 7u 

130  * “I3C 

FASTI- 

89 

W31  « -5II3C)  * W32  * SI  1 ICC  ! 4 1,0  TO  2C0 

FASTF 

90 

c 

11  « PI 

1 Ho 

W31  » -1*  4 W 12  ■ i)«  1.  C>0  TO  2u0 

FASTF 

91 

c 

13  IN  3R0  QUADRANT 

rt 

1 90 

note  • nt* net  * net  * -net 

FASTF 

92 

W11  « -SlliCCCI  i W32  * -SII3CCI 

FASTF 

9 3 

2 uu 

1 L A s T ■ IL+JJ 

FASTF 

94 

DO  220  1 » JO*  HAST i IUIF 

F A A 1 f 

93 

KLAST  • KL+I 

FASTF 

96 

no  220  K ■ I *KL AST  * ? 

FASTF 

47 

K 1 » K + KMlT  4 K ? • K 1 -*  K It  | T « )">  • K2  + >,I!|I 

F ASTI 

9H 

c 


C DO  TWO  STEPS  WITH  J NUT  u 

C 

C A ( K 1 ■ A ( K I + A I K2 ) *W2  * AIK?I  « AIK)  - A(k2)«w? 

C A ( K 1 ) » A(K1>»W1  + AIKH»,V)i-  Al(1|  » AIk1)»'.\1  - A(K3l«w1 

C 


c 

AUI  • AIK)  ♦ A 1 K 1 ) • 

i*  U 1 ) 

■ AUI 

- AIK)  1 

c 

r 

A 1 K 2 • AIK2)  + A 1 K 3 ) * 1 

• AIK.3) 

• A 1 K 2 ) 

- AIKS)*) 

V 

K ■ A I K2  1 *W2 1 - AIK2  + 1 1 »W22 

f-  A r*  T r 

T « AIK/I *W22  + A ( K 2+ 1 I »W2 1 

h A.srr  ln-1 
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AIK2)  » AIKI-R 
AIK2+1 I • A (K*l )“T 


AIK)  ■ AIKHH 
AlKHI  ■ A t K + 1 ) »T 


R • AIK3 )*W31  - A I K 3*1 1 *W32 

T - A(K3)#W32  + A I K3*l  I *W3 1 

AMR  « A I <1 ) *Wl 1 - AIK1+1)«W12 

Awl  » A ( K 1 1 *Wl 2 * AIKl+1 l*Wll 

A ( K, 3 1 • AWR-R  * AIK1I  • AWR*R 

AIK3*U  > AWl-T  » A I K 1* 1 1 « AWi*T 

T » A(<1)  * A ( K 1 1 ■ AIKI-T  * 

T • AIKl+l)  * A I K 1+ 1 ) • AIK+ll-T 
R ■ —A  l K. 3+ 1 I S T ■ AIKS) 

AIKS  I ■ AIK2I-R  $ A ( K2 ) ■ AtK?)*R 


AIK) 

% 


» A I K )*  T 
AlK  + l)  < 


AIK+ll+T 


A I K2* 1 1 • A(K2*1)  + T 


A(K3*1I  • AIK2+1I-T  % 

220  CONTINUE 

f END  OP  I AND  K LOOPS 

230  JJ  ■ JJ0IF'*JJ 
C END  OF  J LOOP 

234  Jl AST  « 4*JLAST  *3 
2 A U CONTINUE 
C END  OF  L LOOP 

240  CONTINUE 
C END  OF  10  LOOP 
C 

C WE  NOW  HAVE  THE  COMPLEX  FOUR  It R SUMS  PUT  THEiN  ADDRESSES  ARt 

C bIT-KEVKHSFID.  THE  FOLLOWING  ROUTINE  PUTS  THEM  IN  UNDER*..** 


FASTF101 

FASTF10? 

FASTF103 
FASTFTO* 
FASTF10J 
FASTF106 
FASTF107 
FASTF10B 
FASTF109 
I AS  TFT 10 
FASTF1 1 1 
FASTF112 
FAS  THU  J 
K A S Y F 1 1 A 

F'ASTFUb 

FAS  TFT  U 
FAS  TFT  17 

F AS  TF 1 18 


NTSQ  ■ NT ANT 

FAS  TFT  IV 

3 50 

IFIM3  *LT.  MTI  370*360 

FASTF120 

360 

N3VNT  ■ N3/NT  i MINNS 

■ NT 

F ASTF12 1 

1003  • U 1 GO  IU  380 

1 As TF 122 

3 70 

N3VNT  « 1 i NTVNS  ■ NT/NS 

* 

MINNS  « 

N3 

F AS  TF1 23 

1003  • 1 

FA6TF12A 

3 80 

JJD3  « NTS0/N3 

F ASTF125 

IFIM2  *LT*  MTI  A70.A6U 

FAS  TFT  26 

A 60 

N2VNT  « N2/NT  * MNM2 

■ NT 

FAST  F'l  27 

1002  « o 1 00  TO  A80 

F A S T F 1 2 8 

A 70 

NlVNT  ■ 1 1 NTVNT  ■ NT/N2 

MINN 2 ■ 

N2 

FASTF12V 

I 002  - 1 

FAST)  1 SI) 

A8C 

JJD2  » NTSO/NZ 

I-  AJ, TFT  31 

IFIMl  *LT»  MT)  STlXDbvi 

FAS  TIT  ii 

36U 

NlVNT  '■  N 1 /NT  E . t NM 

» NT 

F A5TF 1 SJ 

1001  ■ u t CiO  TO  SbO 

FAST)  ) 34 

5 70 

NlVNT  » 1 A MTVN1  ■ F.T/N1 

* 

MlNNl  * 

.N1 

F A S T F 1 3 3 

1001  » 1 

F AS  TFT  36 

981. 

JJO 1 ■ NThU/Nl 

F AS  1 F 117 

600 

J J S ■ J ■ 1 

FAS  TFT  38 

(70  880  JPP3  ■ liHUVNT 

) A.STFlSV 

IPP.3  ■ INVIJJ3I 

F ASTI- 1 AO 

DO  87u  JP 3 ■ l.'UNNS 

) A),  THAI 

IF  1 1003  1 621). 61u 

F A S T F 1 A < 

610 

1 P 3 • INVI  JIM  )»l,3VNT  » 

Go  TU 

63  U 

) AS  TFT  AS 

62u 

IPS  ■>  1 MV  I JP  1) /NTVNS 

f AS  t FT  A 4 

610 

13  ■ I t PP’3+  1 \>'S  I *N2 

F A.',  T IT  A 4 

7 0 0 

JJ2  • 1 

F a S T F 1 A 6 

DO  870  JPP2  • 1.N2VNT 

F ASTI  167 

IPP2  ■ INVIJJ21+13 

1 AS  T )T  A 6 

DO  860  JP2  • l.MIMNr’ 

r As  DU1/ 

IMI002I  7?i;*71  = > 

1 As  til  *>(' 

710 

1 P 2 ■ INVI JP?)*N2VNT  > 

Go  Tu 

V }\> 

F AS T F 1 VI 

7 20 

1 P 2 « INVI JP2I/N7VN2 

1 AfTI  1 4,' 

710 

12  ■ I 1 PP2+ 1 P2 ) «N 1 

1-  A s 1 1 1 6 ) 

BOO 

JJ1  • 1 

1 AS  T 1 1 V A 

DO  86o  JPPI  ■ 1 .NlVNT 

i <> :,  1 1 l v ,i 

IPP1  ■ INVI JJ1 1 + 1? 

1-  AsTI  1 Sli 

DO  H50  JPl  ■ 1 * M t N M 1 

1 As  Tl  1 a / 

IE  I |GOt ) 820*810 

I-  as  TIT  6 A 
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a iu 
8 20 
830 

840 

8 4 5 
830 

8 60 

C 

8 TO 
C 

8 80 
C 

c 

c. 

882 
8 84 

C 

c 

voo 

9U4 

003 

806 

C 

C 

c 

c 

c 


4 20 

9i.fl 

9 50 

C 

c 

c 

c 

9 M' 
C 


191  ■ INVl JP1 )*N1VNT  * 00  TO  830 

IP1  ■ INVl JP1 l/NTVNl 
I * 2* I 1 PP1* 1 PI ) ♦ 1 
tP(J  .GE.  I)  843*840 

f “ All!  3 All)  » A I J I i At  J)  * T 
T ■ A 1 L+l  I ' S A(  1 + 1 ) » A ( J + l I » Alj+1)  * T 
J ■ J*2 
CONTINUE 

jji  ■ jji  + join 

ENO  OF  JPP1  AND  JP2  LOOPS 
JJ2  ■ JJ2  + JJP2 

END  OK  JPP2  AND  JP3  LOOPS 
JJ3  * JJ3  ♦ J JUS 

END  OF  JPP3  LOOP 

WAS  THIS  A TRANSFORM,... 

IF(-IFS)  9999  *4994, lid? 

YES,  REPLACE  A nY  COMJl'ilA), 

DO  884  I ■ 2.NTOT2.2 
All)  ■ -All) 

GO  TO  9V99 

THAT  IS  THE  END  ....... 

ENTRY  SETUP 

THIS  PROGRAM  COMPUTES  iHt  AIN  AM)  1NV  TAnLES 
MT  « XMAXOFI J * | XMAXOKIMI 1 ) ,W( 2 I iKO ) )-2 ) I 
IFIMT  ,LE.  11)  906,905 
IKERR  • 1 1 Go  TO  9949 

1PERR  • 0 

NT  « 2**IIT  * NTV2  • NT/2 
SET  UP  SIM  TAIILti 

THETA  * PI/2»»IL+1)  ...  FUR  L“  1 ... 

theta  « ,2s»pi 

JSTEP  • 2**1  MT“L* 1 ) ...  Fui<  L«1  ... 

JSUP  • NT 

JO  1 F ■ 2** IMT-L I ...  FOR  L • 1 ... 

JO  I F • NT  V2 


S( JD1H  ■ SINK  I THETA  I 
DO  9$u  L ■ 2 * M T 
THETA  • , 5#TME T A 

J S T E P 2 • JSTtP  1 JSH'P  ■ JI'IF 
SIJD1KI  • SlNF(THFTA) 

JC1  • NT-JOIF 

51  JU  ) « COSF  I Thl;  TA) 

JLAsf  * NT-JSTEP2 
I FIJI, AST  ,LT.  jr.TEM)  91)0,92. 

DO  94. i J • JliThP*  JLAST  tJSTLI- 

JONT-J  J JD-J+JulF 

SUDI  » S(JI*SIJC1>  + M JhlF  ) #5!  jr  ) 

CONTINUE 


)■  J.'ll  • J.STFP/2 


SLT  UP  INv(J)  TAHLI 

M T L E X P ■ 2**(KT-LI  ...  Fi.<  L»1  ... 
MTLEXP  ■ MTV? 

LM1EXP  * 2 11  * I L - 1 I ...  K.K  Ul  ... 

L ‘ 1 1 L X P ■ 1 
INVII)  ■■) 

UO  9«a  L»l»r,T 

INV  I L'-.ll'XP*!  1 ■ ITLFXP 

DO  9 Tu  J • /.E-llEXP 

JJ«J+LHlrxP 


FASTF159 
T AS  TF 1 60 
FASTF161 
FfSTF162 
FASTF163 
FAS TF1 64 
FASTF165 
F ASTE1 66 
FASTF167 


FASTI  1 68 
FAS  TF  169 


FASTF17D 

FASTF1T1 
FAS  TF 1 72 
FASTFl  73 

FAST  FI  V4 

FASTF1  75 
FAS  TF 1 76 
FAST F 177 
FAST F 178 
FASTFJ7V 


F ASTF1 HO 

FASTF1H1 

fasti-i«2 

F AS  TF  1 U> 
FASTF184 
F AS  TF | 85 
FA5TF186 
PASTE! 8 t 
I A,’.  PIN 
F AS  IF J 89 
F AS  Tl-  1 90 
F AS tl  1 41 
F uSTI  l 92 
FAS  II  14.! 
FASTI- 144 
F.iSTI  195 


i Asti-  l 4t, 

fas  tl-  I VI 
I ASTI  1 4 . 'I 

i As  tf  1 y-i 

I "S  Tl  2D0 
F-  .Tl  /('l 
FAsTI-,  :V 
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97; 

9 80 
9 999 


invi  jji  ■ I n v c ji+'^tu  xt- 

s'  T L E XI1  ■ KTLl  XI'/2  l.'11-.xi-  * LUUi'f; 

CONT  t NUf 

HRTllMN 


FASTI  ?0J 
FASTFilOA 
F AF.  TF/MVj 

FA5TFZ06 


END 


FA.9TF207 
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